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Abstract 

The STIRAP-based unitary decelerating and accelerating processes have 
been proposed to realize the time- and space-compressing processes in the quan- 
tum control process to simulate the reversible and unitary state-insensitive halt- 
' ing protocol (Arxiv: quant-ph/0607144). A standard three-state STIRAP pulse 

sequence may act as a basic unitary decelerating sequence or a basic unitary 
accelerating sequence. A STIRAP-based unitary decelerating (accelerating) pro- 
cess then consists of a train of these basic STIRAP unitary decelerating (acceler- 
^ ' ating) sequences. The present work is focused on investigating analytically and 

O^. quantitatively how the momentum distribution of a momentum superposition 

state of a pure-state quantum system such as a momentum Gaussian wave- 
|- [ packet state of a single freely moving atom affects the STIRAP state transfer 

J> ■ in these decelerating and accelerating processes. The complete STIRAP state 

, transfer and the unitarity of these processes are stressed highly in the investi- 

gation. It has been shown that the momentum distribution has an important 
influence upon the STIRAP state-transfer efficiency. In the ideal adiabatic con- 
dition these unitary decelerating and accelerating processes for a freely moving 
atom are studied in detail, and it is shown that they can be used to manipulate 
and control in time and space the center-of-mass position and momentum of 
\ a Gaussian wave-packet motional state of a free atom. Two general adiabatic 

conditions for the basic STIRAP decelerating and accelerating processes are 
derived analytically. They are strict and accurate. They can be used to set 
■ up a conventional STIRAP state-transfer experiment and also the basic STI- 

^ I RAP decelerating and accelerating processes. With the help of the STIRAP 

theory and the unitary quantum dynamics it confirms theoretically that the 
time- and space-compressing processes of the quantum control process (Arxiv: 
[quant-ph/0607144) can be realized almost perfectly by the STIRAP-based uni- 
tary decelerating and accelerating processes in the ideal or nearly ideal adiabatic 
condition. 
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1. Introduction 

The stimulated Raman adiabatic passage (STIRAP) processes are very im- 
portant coherent double-photon processes [1, 2]. The STIRAP method has been 
extensively applied to the complete population transfer in energy levels of atoms 



1 



and molecules in the laser spectroscopy [3, 4], the laser cooling in neutral atomic 
ensembles [5, 6, 7, 8, 9] (here may include the conventional Raman adiabatic 
processes), and the atomic quantum interference experiments [10, 11, 12, 13, 14] 
as well as other science research fields. A standard three-state STIRAP pulse 
sequence [3, 4, 15, 17, 18] consists of a pair of Raman laser light beams which 
are selectively applied to three chosen energy levels of an atomic or molecular 
system. The largest advantage of the STIRAP method is that the complete 
population transfer or state transfer in an atomic or molecular system may be 
achieved by the STIRAP pulse sequence with delayed and overlapping Raman 
laser light beams [15, 16, 17], and the STIRAP method is tolerant to the ex- 
perimental imperfections. The basic theory of the STIRAP method has been 
well set up [16, 17, 18]. The experimental confirmation for the STIRAP method 
has been carried out first in the atomic and molecular laser spectroscopy [3, 15] 
and then in diverse other science research fields. The basic STIRAP theoret- 
ical and experimental methods have been first extended to study and design 
the quantum interference experiments in cold atomic ensembles [10, 11, 12, 13, 
14], which involve the atomic motional momentum transfer. The theory also 
has been developed to study the Raman laser cooling processes in a neutral 
atom ensemble [20, 21] by combining the velocity-selective coherent population 
trapping [5, 19]. In these two types of experiments [5, 10, 11, 12, 13, 14, 19, 
20] the Raman laser pulse sequence such as the STIRAP pulse sequence usually 
consists of a pair of counterpropagating Raman laser light beams. 

It has been shown that the dynamical state-locking pulse field plays a key 
role in constructing the reversible and unitary state-insensitive halting protocol 
[22] and solving efficiently the quantum search problem [22, 36]. A general 
state-locking pulse field [22] consists of a sequence of the time- and space- 
dependent electromagnetic pulse fields and could also contain the time- and 
space-dependent potential fields which could be generated by the external elec- 
tric and / or magnetic field. A unitary decelerating (or accelerating) laser light 
pulse sequence that is used to decelerate (or accelerate) a moving free atom 
could be thought of as the component of a dynamical state-locking pulse field. 
The standard three-state STIRAP pulse sequence may act as either a basic uni- 
tary decelerating sequence or a basic unitary accelerating sequence, which is 
dependent upon the parameter settings of the two counterpropagating Raman 
adiabatic laser light beams of the STIRAP pulse sequence. The STIRAP-based 
decelerating (accelerating) process then consists of a train of the basic STIRAP 
unitary decelerating (accelerating) sequences which are applied to a moving 
atom consecutively. It has been proposed that the STIRAP-based unitary decel- 
erating and / or accelerating processes may be used to coherently manipulate the 
halting-qubit atom in the quantum control process to simulate the reversible and 
unitary state- insensitive halting protocol [22] . Thus, the STIRAP-based unitary 
decelerating and accelerating processes are the important building blocks of the 
reversible and unitary state-insensitive halting protocol and the efficient quan- 
tum search process based on the unitary quantum dynamics in time and space. 
A unitary decelerating process is much like the conventional laser cooling pro- 
cess in an atomic ensemble. The essential difference between the two processes 



2 



is that the unitary decelerating process is reversible and unitary, while the laser 
cooling process usually is irreversible. A unitary accelerating process could be 
more like the momentum transfer process in the quantum interference experi- 
ments of an atomic ensemble. The atomic momentum transfer process in these 
quantum interference experiments is generally transverse with respect to the 
initial atomic moving direction, while here the unitary accelerating process is 
longitudinal. The unitary decelerating and accelerating processes stress their 
unitarity and complete state transfer in the quantum control process. 

A quantum computation tends to avoid using a space-dependent unitary op- 
eration as its basic building block, since such a unitary operation usually is more 
complicated and inconvenient to manipulate and control in a quantum system 
with respect to the conventional quantum-gate operations. It is necessary to 
manipulate and control at will the internal motion of an atom (or atomic ion) 
in quantum computation when the specific internal states of the atom such as 
the hyperfine ground electronic states or the nuclear spin states are taken as a 
quantum bit, but at the same time the atomic centcr-of-mass motion tends to be 
kept unchanged simply (or to be constrained simply) so that it does not affect 
these quantum-gate operations and hence in quantum computational theory it 
is usually not considered explicitly. However, as pointed out in the previous 
paper [22], it is of crucial importance to unitarily manipulate and control at 
will in time and space the center-of-mass motion, the internal motion, and the 
mutual cooperation and coupling of the two motions of the halting-qubit atom 
in order to realize the quantum control process to simulate the reversible and 
unitary state-insensitive halting protocol. The unitary manipulation and con- 
trol in time and space for the atom is also a key step toward the realization 
to solve efficiently the quantum search problem. An electromagnetic wave field 
such as a laser light field can manipulate and control in time and space not 
only the center-of-mass motion and the internal electronic (or spin) motion of 
an atom separately, but also it can create and control the coupling between 
the center-of-mass motion and the internal motion of an atomic system. This 
is the theoretical fundament for the laser cooling of a neutral atomic ensemble 
[21, 23, 24] and the unitary decelerating and accelerating of a free atom [22]. 
One large advantage to use a laser light field to manipulate an atom is that 
the space-selective and / or the internal-state-selective unitary operations of 
the atom can be realized easily. The STIRAP-based unitary decelerating and 
accelerating processes could be a very useful double-photon method to coher- 
ently manipulate and control in time and space the center-of-mass motion, the 
internal motion, and the coupling of both the motions of a moving free atom. It 
usually uses a pair of Raman adiabatic laser light beams to couple the atomic 
center-of-mass motional state and internal electronic states (or spin polarization 
states) to realize the coupling between the two motions. The STIRAP-based 
unitary decelerating and accelerating processes have been proposed to realize 
the unitary time- and space-compressing processes which are necessary compo- 
nents of the quantum control process to simulate efficiently the state-insensitive 
reversible and unitary halting protocol [22]. On the other hand, as far as the 
Gaussian wave-packet state of a free atom is concerned, the STIRAP-based 
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unitary decelerating and accelerating processes could generate a type of time- 
and space-dependent unitary propagators which can manipulate and control 
the center-of-mass position and momentum of the atomic Gaussian wave-packet 
state. This important result will be shown in the paper. 

The basic STIRAP theory generally does not consider explicitly the effect of 
the center-of-mass motional momentum distribution of an atomic or molecular 
system on the STIRAP population transfer in the laser spectroscopy [16, 17, 
18], although the Doppler effect in these atomic and molecular systems has also 
been considered suitably in the STIRAP experiments [3, 15]. This could be 
due to that the center-of-mass motion of an atom or molecule is generally much 
slower than the internal electronic motion of the atom or molecule and could be 
neglected in the basic STIRAP theory, and the STIRAP experimental settings 
[3, 15] are also favorable for the experiments to minimize the center-of-mass 
motional effect, for example, the STIRAP experiments may use a pair of co- 
propagating Raman laser light beams. However, a pure-state quantum system 
such as a single freely moving atom may be in a superposition of the center-of- 
mass motional momentum states which may has a broad momentum distribu- 
tion. It is generally hard to realize the complete state transfer in a quantum 
system with a broad momentum distribution by the standard STIRAP method. 
One therefore must consider the effect of the momentum distribution of a su- 
perposition state of the quantum system on the STIRAP population or state 
transfer. In concept the momentum distribution of a momentum superposition 
state of a pure-state quantum system is essentially different from the conven- 
tional momentum distribution of a quantum ensemble, the latter is a statistical 
distribution of momentum of the particles which form the quantum ensemble. 
But the effect of the momentum distribution on the STIRAP state transfer is 
similar for the two types of quantum systems according to the unitary quantum 
dynamical principle that both a closed quantum system and its ensemble obey 
the same unitary quantum dynamics [22] (it seems that a quantum system in 
the presence of an external electromagnetic field is not a closed quantum system, 
but in theory such a quantum system may be treated conditionally as if it is a 
closed quantum system, as shown in the section 11 in Ref. [40]). It has been 
investigated in theory how the momentum distribution of an atomic ensemble 
affects the population or state transfer in the atomic laser cooling based on the 
velocity-selective coherent population trapping [5, 19, 20] which also uses the 
Raman laser light beams, the STIRAP-based momentum transfer in the cold 
atomic interference experiments [10, 12], and the conventional Raman-laser- 
light-based cold atomic interference experiments [11], but these investigations 
are usually either qualitative or numerical. It is important to investigate ana- 
lytically and quantitatively how a superposition of the momentum states affects 
the STIRAP state transfer in a pure-state quantum system when the STIRAP 
method is used to perform the state transfer or the unitary operation in quan- 
tum computation. This paper is devoted to such a theoretic investigation: how 
a superposition of the momentum states affects the STIRAP state transfer in 
a single freely moving atom. According to quantum mechanics a freely moving 
atom in the presence of the STIRAP pulse sequence may be described by the 
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complete set of the product states of both the center-of-mass motional states 
and the internal states (the electronic states or spin polarization states) of the 
atom. The atomic center-of-mass motional state may be a superposition of the 
atomic momentum cigcnstatcs. When the atom is transferred from one internal 
state to another by a STIRAP pulse sequence, the transfer efficiency is gen- 
erally dependent upon the atomic center-of-mass motional state. The purpose 
to investigate this dependence is to understand how the momentum distribu- 
tion of the center-of-mass motional state affects the transfer efficiency and then 
design a better STIRAP pulse sequence so as to achieve a complete STIRAP 
state transfer over the whole effective momentum distribution of the atomic 
center-of-mass motional state. 

In order to manipulate and control the halting-qubit atom in time and space 
in the quantum control process [22] it is necessary to investigate the time evo- 
lution process of the halting-qubit atom in the STIRAP-based unitary decel- 
erating and accelerating processes. In order to investigate quantitatively how 
a superposition of the momentum states affects the STIRAP state transfer for 
a single freely moving atom in the STIRAP-based decelerating or accelerating 
process it is also necessary to calculate the time evolution process of the atom 
in the presence of the Raman laser light field. For example, it needs to solve 
the unitary dynamical equation to set up a general adiabatic condition for the 
basic STIRAP decelerating and accelerating processes for the atom. It seems 
that a single freely moving atom in the presence of an external electromagnetic 
field such as the Raman laser light field is a simple quantum system, but the 
time evolution process of the atomic system is not so simple, partly because the 
interaction between the atom and the external electromagnetic field is usually 
time-dependent, and on the other hand, because the atomic internal motion 
and center-of-mass motion as well as the coupling of the two motions induced 
by the external electromagnetic field need to be considered explicitly and si- 
multaneously in a theoretical treatment. The time evolution process of a free 
atom in the presence of a laser light field becomes so complex that it is generally 
difficult to solve exactly the unitary dynamical equation of the atomic system. 
It is also quite inconvenient even to use an approximation method to solve the 
unitary dynamical equation with a high accuracy. The STIRAP-based unitary 
decelerating and accelerating processes for a free atom are relatively simple, be- 
cause these decelerating and accelerating processes are adiabatic and the Raman 
laser light beams of the STIRAP decelerating and accelerating pulse sequences 
affect only the three chosen atomic internal states. These special points may 
simplify greatly the investigation of the time evolution process of the atom in 
the STIRAP decelerating and accelerating processes. Because the time evolu- 
tion process of the atom involves only the three chosen internal states of the 
atom, it may be investigated conveniently in the atomic three-internal-state 
subspace. On the other hand, it is well known that a unitary process of a free 
atom absorbing or emitting a photon has to obey the energy, momentum, and 
angular momentum conservation laws. The energy, momentum, and angular 
momentum conservation laws put a restriction on the time evolution process 
for the free atom during the STIRAP-based unitary decelerating and accelerat- 
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ing processes. These laws lead to that the atomic motional momentum cannot 
be changed arbitrarily by these Raman laser light beams in the STIRAP uni- 
tary decelerating and accelerating processes, but rather it can be changed only 
within the Raman-lascr-light-induccd momentum state subspace [5, 19]. Then 
the time evolution process of a free atom may be investigated conveniently in the 
Raman-laser-light-induced momentum state subspace in the STIRAP unitary 
decelerating and accelerating processes. Here the unitarity of these processes is 
emphasized again. This greatly simplifies the evaluation for the time evolution 
process. 

In this paper the basic STIRAP theory [15, 16, 17, 18] has been developed 
to study and construct the STIRAP-based unitary decelerating and accelerat- 
ing processes for a free atom by combining the quantum superposition principle 
[25] and the energy, momentum, and angular momentum conservation laws for 
the atomic photon absorption and emission processes [5, 19]. This research is 
focused on investigating analytically and quantitatively the effect of the mo- 
mentum distribution of a superposition of momentum states of the atom on the 
STIRAP state transfer in these processes. Both the ideal and the real adia- 
batic condition for the basic STIRAP decelerating and accelerating processes 
are derived analytically. One important result of the paper is to confirm theoret- 
ically the time- and space-compressing processes of the quantum control process 
[22], which are realized by the STIRAP unitary decelerating and accelerating 
processes, with the help of the STIRAP state-transfer theory and the unitary 
quantum dynamics. 



2 . The Hamiltonian for a single atom in the presence of the Raman 
laser light beams 

The STIRAP-based laser cooling processes and the unitary decelerating and 
accelerating processes are generally involved in the center-of-mass motion, the 
internal electronic motion, and the interaction between the two motions in an 
atomic system. A complete theoretical description for these processes need 
consider the atomic center-of-mass motion, the internal electronic motion, and 
the coupling of both the motions of the atom. When an atom is irradiated by 
an externally applied electromagnetic field such as the Raman laser light beams 
in the STIRAP experiments, the total Hamiltonian for the physical system 
consisting of the atom and the electromagnetic field may be generally written 
as [25] 

H = Ha + Hrad + Hint (1) 

where Ha is the atomic Hamiltonian in the absence of the externally applied 
electromagnetic field, 

k 

Hrad the Hamiltonian for the electromagnetic field which may be written as 

Hrad = ^ I d^x{E.E* +B.B*), 
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and Hint the interaction between the atom and the externally applied electro- 
magnetic field, 




Airk,tf. 



The atomic Hamiltonian Ha describes both the atomic center-of-mass motion 
and the internal electronic (or spin) motion of the atom, while the interac- 
tion Hint induced by the external electromagnetic field creates the coupling 
between the atomic center-of-mass motion and the internal electronic motion. 
The external electromagnetic field usually is weak in the conventional STIRAP 
experiments and the interaction HinL between the atom and the electromagnetic 
field could be considered as a perturbation. 

It is usually inconvenient to use directly the total Hamiltonian of Eq. (1) to 
describe the STIRAP-based decelerating and accelerating processes, although 
it is the most exact to use the total Hamiltonian to treat these processes the- 
oretically. Without losing generality one may use a simpler form of the total 
Hamiltonian of Eq. (1) to describe clearly the unitary decelerating and accelerat- 
ing processes. It is well known that the semiclassical theory of electromagnetic 
radiation has been extensively used to describe the laser spectroscopy in the 
atomic and molecular systems [1, 16, 17, 18], the atomic coherent laser-cooling 
processes [20, 21, 23], and the atomic quantum interference experiments [10, 11, 
12]. If the semiclassical theory is also reasonable for the unitary decelerating 
and accelerating processes, then the Hamiltonian Hrad of the electromagnetic 
field itself may be omitted from the total Hamiltonian of Eq. (1). Though 
the semiclassical theory of electromagnetic radiation generally can not describe 
exactly the spontaneous emission in an atomic system [1, 25] and especially the 
spontaneous emission in the long-time atomic laser cooling process [26] , it may 
be suited to describe the STIRAP-based unitary decelerating and accelerating 
processes as these unitary STIRAP processes can avoid the atomic spontaneous 
emission by setting the suitable experimental parameters. Furthermore, if the 
electric dipole approximation for the atomic system is also reasonable, then one 
may use conveniently the electric dipole interaction H^ to replace the interac- 
tion Hint to describe the unitary decelerating and accelerating processes. The 
electric dipole interaction H^ may be written as 



where D is the atomic electric dipole moment and E{x, t) is the electric field of 
the externally applied electromagnetic field, and the coordinate x in the electric 
field E{x, t) is the center-of-mass position of the atom. The electric dipole ap- 
proximation is reasonable when the wave length of the external electromagnetic 
field is much larger than the atomic dimension under study so that the atom 
may be considered as a point particle — a point electric dipole — in the electro- 
magnetic field [1]. It has been shown that in the electric dipole approximation 
the electric-dipole Hamiltonian Hd is really equivalent to the interaction Hint 
up to a gauge transformation which is also a unitary transformation [1, 27]. 



Hd = -D.E{x,t) 



(3) 
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The electric dipole approximation is very popular in the theoretical description 
of a variety of laser light and matter interactions. For example, one generally 
uses the electric-dipole approximation to deal with theoretically the dynamical 
process of atomic laser cooling in an atomic ensemble [1, 4, 5, 10, 19, 20, 26]. 
Now in the semiclassical theory of electromagnetic radiation and in the elec- 
tric dipole approximation the total Hamiltonian (1) of the atom in the external 
electromagnetic field is reduced to the form 

H = PV(2M) + V{x) + H{r) + (4) 

where the sum of the first three terms is just the atomic Hamiltonian Ha of 
Eq. (2). The center-of-mass Hamiltonian Hem = P'^/{2M) + V(x) describes the 
atomic center-of-mass motion, while the internal Hamiltonian H(r) describes 
the internal electronic (or spin) motion of the atom. In the center-of-mass 
Hamiltonian Hem the term Hk = P'^/{2M) with the atomic mass M and mo- 
mentum P is the atomic kinetic energy and the term V{x) the atomic potential 
energy in an external potential field. In the unitary decelerating and acceler- 
ating processes of a free atom the external potential energy V{x) of the atom 
is zero, i.e., V{x) = 0. If there is an external electric (or magnetic) field during 
these processes, then the external potential energy V{x) could not be zero. For 
example, the potential energy V{x) ^ when an atom in a harmonic potential 
well is applied by an external electromagnetic field. Generally, both the center- 
of-mass kinetic energy operator P^/(2M) and the internal Hamiltonian H{r) 
are not commutable with the electric dipole interaction Hd. The Hamiltonian 
H of Eq. (4) still needs to be further simplified so that it can describe con- 
veniently the standard three-state STIRAP experiment. This simplification is 
based on the facts that the externally applied electromagnetic field in the STI- 
RAP experiments can only affect some specific internal electronic energy levels 
of the atom and the internal electronic states are discrete. In the following only 
one-dimensional center-of-mass motion of the atom is considered. 

Because there are the center-of-mass motion, the internal electronic motion, 
and even the coupling between the two motions of an atom in the STIRAP 
processes, one must use simultaneously both the atomic center-of-mass motional 
states and the internal electronic states of the atom to describe exactly the 
STIRAP processes. Suppose that the wave functions \il)p{x)) and \ipj{r)) are the 
eigenstates of the center-of-mass Hamiltonian Hem and the internal Hamiltonian 
H{r), respectively, 

Hem\i^p{x))=Ep\i;p{x)), i?(r)|V,-(r)) (r)), 

where Ep and Ej are the eigen-energy of the center-of-mass motional state 
\tpp{x)) and the internal state \tpj{r)), respectively. For the STIRAP processes 
of a free atom the eigen-energy Ep is equal to the atomic center-of-mass ki- 
netic energy P'^/2M since the potential energy V{x) = and Hcmlijjpix)) = 
HKli^pix)) = P^ / {2M)\ip p{x)) . According to quantum mechanics [25] the com- 
plete set {\'>Pp{x))\tpj{r))} of the product states of the center-of-mass motional 
states and the internal states of the atom can be used to describe completely 
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both the ccntcr-of-mass motion and the internal electronic motion as well as the 
coupling of the two motions of the atom. Now the total wave function \ 'i/{x,r,t)) 
of the atom may be generally expressed as a linear combination of these product 
states [25], 

\^{x, r, t)) = ^ a{j, P, t)IVp(x)) |V,(r)) 
i,P 

where a{j, P, t) is an expansion coefficient. In the decelerating and accelerating 

processes an atom may be in a given internal state, while its ccnter-of-mass 
motional state is a superposition state. Then in this case the total wave function 
of the atom may be expanded as 

\^{x,r,t)) =^a{j,p,t)\%{x,t))\i)j{r)) (5) 

where the wave function \iljp{x, t)) isa superposition of the atomic center-of-mass 
motional states. For example, the center-of-mass motional state \^pp{x,t)) may 
be a wave-packet state of the halting-qubit atom in the quantum control process, 
which was denoted as \CM,R) in the previous paper [22]. In the decelerating 
and accelerating processes of a free atom any center-of-mass motional state 
\'tpp{x,t)) of the atom may be expanded in terms of the complete eigenstate set 
{|^p(x))} [25], 

=^a(p,P)exp(-^^i)|Vp(a;)). (6) 

p 

Note that the state \tjjp{x)) is also an eigenstate of the atomic center-of-mass 
momentum operator P and in one-dimensional case it may be written as 

|^p(x)) = |P) = -^exp(iPa;/fi). (7) 

The momentum wave function |P) represents that the atom moves along the 
direction +x with the center-of-mass motional velocity P/M. 

In what follows it is supposed that the halting-qubit atom has a three-level A 
configuration for the three-state STIRAP experiments. In the STIRAP experi- 
ment the external electromagnetic field can have a real effect only on the specific 
three-state subspace {|V'o(''))> IV'i(^))) IV'2(''))} of the internal electronic states 
of the atom. This three-state subspace will further simplify the Hamiltonian 
of Eq. (4) as the time evolution process of the internal states of the atom in 
the STIRAP experiment is confined in the three-state subspace. Those internal 
states of the atom outside the three-state subspace will not be affected by the 
external electromagnetic field and are not considered in the STIRAP experi- 
ment. Thus, it is sufficient to use the internal Hamiltonian {H{r)) projection 
onto the three-state subspace to describe the three-state STIRAP experiment. 
Since the internal states {\ijjj{r))} are the cigcnstatcs of the internal Hamilto- 
nian H{r), this projection Hamiltonian onto the three-state subspace may be 
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given by 



H{r) = Eo\Mr)){Mr)\ + E,\i,,{r)){i,,{r)\ + E2\i>2{r)){i,^{r)\. (8) 

On the other hand, the electric dipole interaction Hd of Eq. (3) can also be 

simplified further in the three-state subspace. In the STlRAP-based decelerat- 
ing and accelerating processes the total external electromagnetic field generally 
consists of a pair of counterpropagating electromagnetic fields, which may be 
amplitude- and phase-modulation Raman laser light beams. The total electric 
field for the two counterpropagating linearly polarized Raman laser light beams 
may be expressed as 

= ^^Lo{t) exp[i{-kLo-x - wlo*)] 

+ ^ELi{t) exp[i{kLi.x - iVLit)] + C.C. (9) 

where C.C. stands for the complex (or hermite) conjugate of the first two terms. 
Notice that the three internal states {IV'fe have quite different energy eigen- 
values. The first Raman laser light beam with the electric field i?Lo(0 couples 
only the two internal states iV'ol^)) ^^'^ IV'2('')) of the atom, while the second 
with the electric field ELi{t) connects only the two internal states \ipi{r)) and 
1^2 (''))• The frequency difference {ujlq — ujliI should be near the resonance 
frequency of the two atomic internal energy levels Itpoir)) and \ipi(r)). The first 
Raman laser light beam usually is named the pumping laser pulse and the sec- 
ond the Stokes laser pulse in the laser spectroscopy [15]. If among the three 
internal states the two states |'0o(^)) {ipiir)) which are usually the ground 
states have the same energy eigenvalues, then one may use a pair of cr+ and 
(7_ circularly polarized laser light beams [19] to replace the present two Raman 
laser light beams, one circularly polarized laser light beam coiipling only the 
two internal states |'0o('')) ^^'^ IV-'2('')) another connecting only the two 
internal states \ipi{r)) and |V-'2(''))- Then in the three-state (internal) subspace 
the electric dipole interaction of Eq. (3) may be written as, in the rotating 
wave approximation [1], 

Hd = nflo2{t) exp{i{-kLo-x - ujLot)}\i'2{r)){ipo{r)\ 

+ nrii^it) cM^kLi.x - LULit)}\4>2{r)){4>i{r)\ + C.C. (10) 
where the Rabi frequencies for the two Raman laser light beams are defined as 

fto2{t) = -l{Mr)\D-ELo{t)\Mr)), 

n,2{t) = -l{Mr)\D.ELim,{r)}. 

The electromagnetic field of the Raman laser light beams is usually weak in 
the conventional STIRAP experiments, so that the effect of the electromagnetic 
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field on the atom could be considered as a perturbation. Hence the rotating 
wave approximation is reasonable. On the other hand, if each Raman laser 
light beams in the STIRAP experiment is replaced with a pair of the laser 
light beams with the orthogonal electric field vectors and the suitable phases 
[38] or one circularly polarized laser light beam [19], then the rotating- wave 
approximation may be eliminated and hence the electric dipole interaction (10) 
may be constructed exactly. The total Hamiltonian of Eq. (4) associated with 
the electric dipole interaction if^ of Eq. (10) and the internal Hamiltonian H{r) 
of Eq. (8) and the product basis set {\4>p{x))\4'j{r))} may be used conveniently 
to describe the STIRAP-based imitary decelerating and accelerating processes 
of a free atom. The transition matrix elements of the electric dipole interaction 
Hd can be calculated in the product basis set, 

W{j',P';j,P) = (Vy(r)K^p,(a;)|ifd|^p(a;))|V,(r)). 

These matrix elements are not zero only when both the internal states \ipj{r)) 
and \ipji(r)) are in the three-state (internal) subspace. They are also subjected 
to the constraint of the energy, momentum, and angular momentum conserva- 
tion laws for the atomic photon absorption and emission process in the STIRAP 
decelerating and accelerating processes. This is an instance of the velocity- 
selective rules which have been used in the atomic laser cooling processes [5, 
19]. Below it is shown how the energy and momentum conservative laws have a 
constraint on the electric dipole transition matrix elements {W{j' , P';j, P)} in 
the three-state STIRAP experiments of the STIRAP-based unitary decelerating 
and accelerating processes. 

In the STIRAP experiments the internal states \ipQ{r)) and |V'i(^)) of the 
three-state subspace usually are taken as the hyperfine ground electronic states 
\go) and \gi) of an atom, while the internal state |V'2('')) be taken as some 
excited state \e) of the atom. For example, the two internal states jgo) and 
l^i) may be the hyperfine ground electronic states 35*1/2 {P = 1) and 3Si/2 
{F = 2) of sodium atom (Na), respectively, while the excited state |e) may be 
the excited electronic state 3P3/2 {F = 2) of the sodium atom. Suppose that at 
the initial time in the STIRAP experiment the atom is in the ground internal 
state l^o) and the center-of-mass momentum state \P') = {V^)~^ exp{iP'x/fi), 
which also means that the atom is in the product state |P')|go) and it travels 
along the direction +x with the velocity v' = P'/M. Now a laser light field 
propagating along the direction —x is applied to the atom. Then the moving 
atom may absorb a photon from the laser light field if the frequency (cu = koc) 
of the laser light field is just equal to the transition frequency of the atom in 
motion between the ground internal state \go) and the excited state |e) after the 
Doppler effect is taken into account. After the atom absorbs a photon from the 
laser light field, the atomic motional momentum becomes P' — Hko according to 
the momentum conservation law and hence the atom is decelerated by Hko/M. 
Since the atomic internal energy levels are discrete, the momentum change of 
the moving atom is also discrete after the atom absorbs a photon. This means 
that when the atom is excited to the internal state |e), its motional momentum 
can not take an arbitrary value but it has to be P' — Hko due to the fact 
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that the atomic optical absorption process obeys the energy and momentum 
conservation laws [19]. After the atom absorbs a photon it jumps to the excited 
state |e) from the ground state l^o) and its initial motional state \P') is changed 
to \P' — hko) and hence the atom is in the product excited state \P' — fifco)|e). 
This is just the atomic decelerating process based on the optical absorption 
mechanism [23, 24]. The atom in the product excited state \P' — fiko)\e) may 
be further decelerated by another laser light field. This laser light field travels 
along the same direction +x as the atom and its frequency {lo = kic, ki ^ ko) 
is equal to the transition frequency of the moving atom between the ground 
internal state \gi) and the excited state |e) after the Dopplcr effect is taken 
into account. Thus, this laser light field may stimulate the atom in the excited 
internal state \e) to jump to the ground internal state \gi). Because the energy 
difference between the two internal states \gi) and |e) is quite different from that 
one between the two internal states \go) and |e), this laser light field will not 
affect the transition between the two internal states j^o) and ]e). Likewise, the 
first laser light field docs not affect the transition between the two internal states 
l^i) and |e). Different from the first decelerating process this atomic decelerating 
process is based on the stimulated optical emission mechanism [23, 24]. An atom 
in an excited state may jump to the ground state when it is stimiilatcd by an 
external laser light field [25]. When the atom in the excited state |e) jumps to the 
ground state \gi), it may emit coherently a photon to the laser light field. Since 
the atomic motion direction is the same as the propagation direction of the laser 
light field, the atom really sends part of its motional momentum to the laser 
light field and hence is decelerated in the stimulated transition process from the 
excited state |e) to the ground state \gi). This part of motional momentum is 
just Hki according to the momentum conservation law and accordingly the atom 
is decelerated by hk\/M. Thus, after the stimulated optical emission process the 
atom is in the ground internal state \gi) and has to be in the motional state 
\P'—hko — hki), that is, the atom is in the product state \P'—hko — hki)\gi). Both 
the atomic optical absorption and emission processes are required to be unitary 
here, as pointed out before [22]. The imitary decelerating process based on the 
three-state STIRAP process just consists of the reversible optical absorption and 
emission processes mentioned above, where the two Raman laser light beams are 
adiabatic and usually counterpropagating. This basic STIRAP-based unitary 
decelerating process may be expressed in an intuitive form 

|P + fifco)|.go) ^ |P)|e) \P-hki)\gi) (11) 

Here for convenience in the later discussion the atomic motion momentum P' is 
denoted as P+ hko. Thus, a conventional three-state STIRAP pulse sequence 
may be really used as a basic unitary decelerating sequence if the two Raman 
laser light beams of the STIRAP pulse sequence are arranged suitably such 
that the moving atom is decelerated consecutively by the two Raman laser light 
beams. Obviously, the inverse process of the STIRAP-based unitary deceler- 
ating process (11) may be used to accelerate the atom in motion. However, it 
is more convenient to use directly the reversible optical absorption and emis- 
sion processes to accelerate an atom in motion and this may be achieved by 



12 



setting suitably the parameters of the two Raman laser light beams of the STI- 
RAP pulse sequence. In the STIRAP-based unitary accelerating process the 
first laser light field that induces the optical absorption process travels along 
the motional direction {+x) of the atom, while the second laser light field that 
stimulates the atomic optical emission process propagates along the opposite 
direction (— x) to the moving atom. When the atom in the initial product state 
|-P')l.9o) is excited to the product state \P' + hkQ)\e) by the first laser light 
field, it absorbs a photon from the laser light field and is accelerated by Hk^/M. 
When the atom in the excited state \P' + fifco)|e) jumps to the ground state 
\P' + Hk'f^ + Hk'i) \gi) under the stimulation of the second laser light field, it emits 
a photon to the laser light field and is accelerated further by tik[/M. Therefore, 
the basic STIRAP-based unitary accelerating process may be expressed in an 
intuitive form {P' = P — ftfeg) 

\P-nk'o)\go) - \P)\e) - \P + fik[)\gi). (11a) 

In what follows only the basic STIRAP-based imitary decelerating process (11) 
is treated explicitly. In an analogous way, one can also deal with the basic 
STIRAP-based unitary accelerating process (11a). 

It follows from the basic decelerating sequence (11) that the time evolution 
process of the atom in the unitary decelerating process (11) is restricted within 
the three-state (product state) subspace {\P+fiko)\go), \P—fiki)\gi)} for 

a given atomic motion momentum P. Then during the STIRAP-based unitary 
decelerating process (11) the total wave function \'^{x,r,t)) of the atom at any 
instant of time t can be expanded in the three-state (product state) subspace [5, 
10, 12, 19, 20], according to the superposition principle in quantum mechanics 
[25], 

I'i'ix, r, t)) = J2 PiPHMP, t)\P + fiko}\9o) 
p 

+ Ai{P,t)\P)\e) + A2{P,t)\P - Hki)\gi)}. (12) 

where the sum over the momentum P is due to the fact that the atom may be 
in a superposition of momentum states, as can be seen in Eq. (6), p{P) is the 
time-independent amplitude which has the physical meaning that |/9(P)P is the 
probability in the superposition to find the atom in the three-state subspace 
{\P + fiko)\go), \P)\e), \P — Kki)\g\)} labelled by the momentum P, and the 
time-dependent amplitudes {Afe(P, t)} satisfies the normalization condition: 

|^o(P,t)|' + \Ai{P,t)\'' + |A2(P,i)|2 = 1. (13) 

The amplitude p{P) is time-independent because the two Raman laser light 
beams induce a change only within the three-state subspace {|-P + hko)\go), 
|P)|e), \P — Kk\)\g\)} for each given momentum P during the unitary decelerat- 
ing process [19]. If at the initial time to the atom is in the ground internal state 
\go) and in a wave-packet motional state, then the initial wave packet state of 
the atom may be expanded as 

\^{x,rM))=Y.P^P')\P')\9,)- (12a) 
P' 
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Obviously, here the coefRcients Ao{P,to) = 1 and Ai{P,to) = A2{P,to) = 0, 
which can be deduced from Eqs. (12) and (12a) with P' = P+Hko, while |p(P')p 
is the probability to find the atom in the three-state subspace {|-P + ^^0)150)) 
|P)|e), \P — fiki)\gi)}. As an example, the initial state \'^{x,r,to)) may be 
taken as the Gaussian wave-packet motional state of the halting-qubit atom in 
the right-hand potential well of the double-well potential field in the quantum 
control process [22] . The three-state (product state) subspace and the basic de- 
celerating process (11) show that only special dipole transition matrix elements 
{W{j',P';j,P)} can take nonzero values, that is, for a given momentum P 
there are only four matrix elements to take nonzero value: W{e, P; go, P + hka), 
W{e, P;gi,P-Hki), W{go,P + hko; e, P), and W{gi, P - hki; e, P). For the ba- 
sic STIRAP decelerating process (11) the total electric field of the two Raman 
laser light beams can be explicitly obtained from equation (9) by setting the 
parameter sets: {ELo{t),kLQ,oJLo) = {Eoi{t),ko,LOo) and {ELi{t),kLi,LOLi) = 
{Ei2{t), ki,cui), where the first Raman laser light beam (£'oi(f), fco, wo) couples 
the two internal states |go) and |e) and its propagating direction is opposite to 
the motional direction of the atom, while the the second beam {Ei2{t), ki,LUi) 
connects the two internal states \gi) and |e) and it travels along the motional 
direction {+x) of the atom. Here suppose that the energy difference (measured 
in frequency unit) between the two ground internal states l^o) and \gi) is much 
larger than the detunings of the two Raman laser light beams. Now these four 
nonzero electric-dipole-transition matrix elements for the basic decelerating pro- 
cess (11) can be obtained with the help of the total electric field of Eq. (9) with 
these parameter settings, the electric dipole interaction Hd of Eq. (10), and the 
momentum eigenstates of Eq. (7) as well as their orthogonalizations, 

Wo2(i) = W^oit) = {e\{P\Hd\P + nko)\go) = ^02(1) eM'i'^ot), 

W2o{t) = {go\{P + fiko\Ha\P)\e), 
Wi2{t) = W^^it) = {e\{P\Hd\P - nki)\gi) = ?ini2(i) exp(-iwit), 

W2i{t) = {g^\{P - fiki\Hd\P)\e). 

Here the star * stands for the complex conjugate and the Rabi frequencies f2o2(0 
andf2i2(i) are defined in Eq. (10) with the states: |V'o('')) = Ido)^ IV'i(^)) = l^i)) 
and |V'2(0) = |e)- 

Now the time evolution process of the atom in the presence of the Raman 
laser light beams in the basic decelerating process (11) is described by the time- 
dependent Schrodinger equation: 

ih^^^ix,r,t) = H{t)-^{x,r,t). (14) 

Here the total Hamiltonian H{t) is given by Eq. (4), in which V{x) = and 
H{r) and Hd are given by Eq. (8) and (10), respectively, while the wave func- 
tion |\l/(a;, r, t)) is given by Eq. (12). By using the four nonzero electric-dipole- 
transition matrix elements and the orthonormalization of the momentum eigen- 
states of Eq. (7) the Schrodinger equation (14) can be reduced to a three-state 
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Schrodingcr equation for a given momentum P, which may be written in the 
matrix form 



ih 



d_ 
dt 



AaiP,t) \ / AoiRt) 

A,{P,t) =H{P,t) A,{P,t) 
A2{P,t) / V A2{P,t) 



(14a) 



where the three-state vector {AQ{P,t), Ai{P,t), A2{P,t))'^ (here T stands for 
the vector transpose) satisfies the normalization of Eq. (13) and the reduced 
Hamiltonian H(P,t) is a 3 x 3— dimensional Hermitian matrix, 



H{P,t) 



2M ' 





So 



p2 

2M +^2 







(P-hki 
2M 



El 



Here the three basis vectors |1) = (1, 0, 0)^, |2) = (0, 1, 0)^, and |3) = (0, 0, 1)"^ 
of the three-state vector space {{AQ{P,t), Ai{P,t), A2{P,t))'^} stand for the 
three basis product states \P + tiko)\go), \P)\e), and \P — tiki)\gi) of the original 
three-state subspacc {\P + Hko)\go), \P)\e), \P — fiki)\gi)}, respectively. This 
reduced three-state Schrodinger equation can describe completely the thrcc- 
statc ST I RAP experiments just like the original Schrodingcr equation (14). 
By making a unitary transformation on the three-state vector in Eq. (14a) [20, 
12, 18, 25]: 

.1 .{P + nkaf 



Ao{P,t) 



exp[-(- 



2M 

32 



Eo)t]AoiP,t), 



Ai(P,i) = eM^i^ + E2)t]Ai{P,t), 



A2(P,t) = cxp -(-^^ — 



Ei)t]A2{P,t), 



the Schrodinger equation (14a) is further reduced to the form 



d 



ih- AiiP,t) 

''^ ^ A2{P,t) 



H{P,t) 



MP,t) 
Ai{p,t) 



(15a) 
(15b) 

(15c) 

(16) 



Now the Hamiltonian H(P, t) is a traceless Hermitian matrix, 



H{P,t) = 



WS2{P,t) 

Wo2{P,t) Wl2{P,t) 

wt2{P,t) 



where the time- and momentum-dependent complex parameters W02 {P, t) and 
Wi2{P,t) are given respectively by 

Wo2{P,t) = mo2(t)exp{-i[ - (a;o2 - u;o)]t}, 



Wi2{P,t) =Hni2{t)exp{- 



2M 

-2Pki + hkl 
2M 



(Wi2 - OJl)]t}, 
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and the transition frequencies for the atomic; internal states arc defined by 
fia;o2 = E2 — Eq and Huii2 = E2 — Ei. Suppose that the Raman laser light 
beams are amplitude- and phase-modulating. Then the Rabi frequencies for the 
two Raman laser light beams (the pumping pulse {flp{t)) and the Stokes pulse 

are written as 

r2o2(i) = r2p(t)exp[-i^o(i)], ili2{t) = ^^(i) exp[-i(/)i(t)], 
and the parameters in the Hamiltonian H{P, t) therefore are given by 

Wo2{P, t) = Mlpit) exp[-iap(P, t)], W^i2(f , t) = Mls{t) exp[-ias(P, i)], 
where the phases ap{P,t) and as(P,t) are dependent upon the momentum P, 

as{P,t) = [ ~'^^\'^^^' - {^12 - uj,)]t + Mt). 
Now the Hamiltonian H{P, t) can be rewritten in the explicit form 



H{P, t) = n 



Qp(t)e'«f(^'*) 

1^,(^)6*"' 



(17) 



This type of Hamiltonians often have been met in the three-state STIRAP ex- 
periments in the laser spectroscopy [4, 15, 17, 18] and in the atomic interference 
experiments [12]. The three-state Schrodinger equation (16) and the traceless 
Hamiltonian (17) are the theoretical basis to design the Raman adiabatic pulses 
of the STIRAP-based decelerating and accelerating processes. There is a special 
point in the STIRAP-based unitary decelerating and accelerating processes that 
the Hamiltonian (17) is dependent upon the center-of-mass momentum of the 
atom besides the frequency offsets {uJk^ — w^} (fc = and 1). This is similar to 
the situations of the STIRAP-based atomic laser cooling [20] and quantum inter- 
ference experiments [10, 12]. The effect of the frequency offsets on the STIRAP 
population transfer has been examined in detail in the laser spectroscopy [18b]. 
Though here considers only the pure-state quantum system of a single atom 
instead of an atomic ensemble, the atomic momentum distribution could have 
a great effect on the population transfer of the atom from an internal state to 
another in the STIRAP decelerating and accelerating processes. This is because 
the atom may be in a superposition of the momentum eigenstates and hence 
has a momentum distribution. For example, though a freely moving atom is in 
a Gaussian wave-packet state in coordinate space, it is also in a superposition of 
the momentum eigenstates of the atom in momentum space. Here the position 
of the momentum P in the unitary decelerating and accelerating processes is 
similar to that one of the frequency offsets in the STIRAP experiments of the 
conventional laser spectroscopy [18]. The frequency offsets (with respect to the 
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transition frequencies of given atomic internal energy levels) could be set at will 
for a single atom, since they are the parameters of the Raman laser light beams, 
while the superposition of momentum eigenstates of the atom (i.e. the atomic 
momentum distribution) is an inherent property of the atomic motional state. 
Whether or not the Raman adiabatic pulses obtained from the Schrodinger 
equation (16) and the Hamiltonian (17) are suitable for the decelerating and 
accelerating processes are dependent upon the atomic momentum distribution. 
The excitation bandwidth for a good STIRAP pulse sequence must be much 
larger than the effective spreading of the atomic momentum distribution. 

In the quantum control process [22] to simulate the reversible and unitary 
halting protocol and the quantum search process the halting-qubit atom needs 
to be decelerated and accelerated by the STIRAP-based unitary decelerating 
and accelerating processes, respectively. As rcqiiired by the quantum control 
process, if the halting-qubit atom is completely in the ground internal state 
l^o) (151)) at the initial time in the decelerating or accelerating process, then 
it must be converted completely into another ground internal state (|.go)) 
at the end of the process. Because the conversion efficiency from the initial 
state to the end state in these processes has a great effect on the performance 
of the quantum control process, it is of crucial importance to achieve a high 
enough conversion efficiency in these processes. This is the first guidance to 
design the STIRAP pulse sequences for the unitary decelerating and accelerating 
processes. On the other hand, in order that a high conversion efficiency is 
achieved in the STIRAP experiments one must also consider the decoherence 
effect due to the atomic spontaneous emission when the atom is in the excited 
internal state. The atomic spontaneous emission becomes an important factor 
to cause the decoherence eff'ect when the atom is in a short-lifetime excited 
internal state in the STIRAP experiments. Therefore, the atom should avoid 
being in the excited internal state during the STIRAP experiments. This may 
be realized by setting the favorable detunings for the two Raman laser light 
beams. Since the three product states \P + ^0)150)) and \P — fiki)\gi) 

are the eigenstates of the total atomic Hamiltonian Ha of Eq. (2) and have 
eigenenergy: {P + Hko)^ / {2M) + Eq, P^/{2M) + E2, and {P - hkif / {2M) + Ei, 
respectively. Then the energy difference between the ground state \P + hkQ)\go) 
and the excited state \P) |e) is given by {E2 - Eq) - Phko/M - n^kl/{2M) and 
the energy difference between the ground state \P — hki)\gi) and the excited 
state |P)|e) is {E2 - Ei) + PHki/M - H'^kf/{2M). Since the Raman laser light 
beam with the carrier frequency luq couples the ground state \P + tiko)\go) and 
the excited state |^')|e), the detuning Ap of the beam is given by 

A,(P) = (.o.-a;o)-(^)P-g, 

while the detuning Ag of another Raman laser light beam with the carrier 
frequency Wi is 

A.(P) = + (|)P-§^ 

By using the two formulae one may set conveniently the detunings Aj, and A^ 
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for the STIRAP-bascd decelerating and accelerating processes. Both the detun- 
ings are dependent upon the motional momentum of the atom and also their 
frequency offsets, respectively. If the atom is in a wave-packet state or generally 
in a superposition of momentum states, then there is a distribution of momen- 
tum and the detuning settings should consider the momentum distribution. In 
next sections the parameters of the Raman laser light beams will be determined 
so as to arrive at the main goal that the conversion efficiency from the initial 
state to the end state in the decelerating and accelerating processes must be 
as close to 100% as possible. There the atomic momentum distribution must 
be taken into account. Thus, this is involved in the adiabatic condition for the 
three-state STIRAP state transfer of the decelerating and accelerating processes. 

3. The basic differential equations for the STIRAP-based decel- 
erating and accelerating processes 

The three-state STIRAP experiments have been studied extensively both in 
theory and experiment in the laser spectroscopy [15, 16, 17, 18]. These studies 
tell ones which conditions the complete population transfer can be achieved by 
the STIRAP method among the energy levels of atomic and molecular systems. 
From the point of view of theory an important feature of the three-state STI- 
RAP experiments is that there is the special eigenstate of the Hamiltonian (17) 
that corresponds to the atomic trapping state [28]. This special eigenstate is 
independent of the intermediate state [17], which is usually taken as the excited 
internal state of the atom in the STIRAP experiments. Under the adiabatic con- 
dition the complete state or population transfer through the special eigenstate 
can occur from the initial state directly to the final state, while the intermedi- 
ate state is bypassed in the transfer process. Such an adiabatic state-transfer 
process is particularly favorable for the STIRAP-based decelerating and accel- 
erating processes. This is because the atomic spontaneous emission could be 
avoided in the decelerating and accelerating processes when the excited state of 
the atom is bypassed in the STIRAP state-transfer processes. This also shows 
that the semiclassical theory of electromagnetic radiation is suited to treat the 
STIRAP experiments. On the other hand, from the experimental point of view 
the three-state STIRAP experiment needs to set suitably the experimental pa- 
rameters for the Raman laser light beams. The important thing in experiment is 
the settings for the Rabi frequencies of the two Raman laser light beams and for 
the pulse delay between the two Raman laser light beams [16] . The atomic sys- 
tem should first interact with the Stokes pulse and then with the pumping pulse, 
and an appropriate overlapping between the two Raman laser light beams is also 
required in experiment [3, 4, 15, 16, 17, 18]. These requirements are generally 
related to the adiabatic condition of the three-state STIRAP experiment. 

The special point for the STIRAP decelerating and accelerating processes 
is that one must consider the atomic momentum distribution when a general 
adiabatic condition is set up for these processes. According to the adiabatic 
theorem and the adiabatic approximation method in quantum mechanics [29, 
30] one should first calculate the three eigenvectors and eigenvalues of the in- 
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stantaneous Hamiltonian ofEq. (17), 



H{P,t)\g{P,t)) = E{P,t)\g{P,t)). (18) 

The three eigenvectors {\g{P,t))} and their eigenvalues {E{P,t)] of the Hamil- 
tonian i?(P, t) are usually named the adiabatic eigenvectors and eigenvalues, 
respectively. Notice that the three basis vectors of the three-state vector space 
{{A^{P,t),Ar{P,t),A2{P,t)Y} are |1) = (1,0,0)^, |2) = (0,1,0)^, and |3) = 
(0,0, 1)"^, respectively. The three basis vectors really stand for the three basis 
product states |P-|-fifco)|(?o), l-P)|e); and \P — hki)\gi) of the original three-state 
subspace for each given momentum P, respectively. Any one of the three eigen- 
vectors of the Hamiltonian iJ(P, t) may be expanded in terms of the three basis 
vectors. By the explicit Hamiltonian of Eq. (17) one can obtain one of the three 
eigenvectors [15, 17, 18], 

\g\P,t)) = exp[-i7(t)]{cos0(t)|l) - sin0(t) exp[-z(ap(P, t) - a,(P, i))]|3)} 

(19a) 

where the phase ^{t) is a global phase, the mixing angle 9{t) and the Rabi 
frequency Q.{t) are defined respectively by 

cos^(t) = sin^(i) = and n{t) = y^y^j^TiW- 

The eigenvalue corresponding to the eigenvector \g^{P,t)) is = fko^ = 0. 
The eigenvector \g^(P,t)) is special in that it does not contain the intermediate 
eigenvector |2) which contains the excited internal state |e) of the atom. It is 
the so-called atomic trapping state [28]. The adiabatic population transfer is 
achieved through the eigenvector |g''(P, t)) in all the three-state STIRAP ex- 
periments [15, 17]. Therefore, when the initial state |1) is transferred to the 
final state |3) through the eigenvector \g°{P,t)) in the STIRAP state-transfer 
process, the intermediate state |2) is not involved and hence is bypassed. Since 
the phase difference ap{P,t) — as{P,t)) in the eigenvector \g^{P,t)) is depen- 
dent upon the atomic motional momentum P the adiabatic state transfer in the 
STIRAP experiments is really affected by the atomic momentum distribution. 
In the laser spectroscopy the similar trapping state is obtained [4, 15, 17, 18] 
but that state is not dependent upon the momentum P. Thus, there is not any 
theoretical problem involved in the effect of the momentum distribution on the 
adiabatic state transfer in the laser spectroscopy. Of course, in the laser spec- 
troscopy the Doppler effect usually is considered in the STIRAP experiments 
of those quantum systems such as an atomic or molecular beam [15]. However, 
in the atomic laser cooling [5, 19, 20, 21], quantum interference experiments 
[10, 11, 12, 13, 14], and the atomic decelerating and accelerating processes the 
atomic momentum distribution generally needs to be considered explicitly. The 
other two adiabatic eigenstates of the Hamiltonian of Eq. (17) are given by [15, 
17, 18] 

\g^{P,t)) = -^exp[-z,5(t)]{sin0(t)|l) Texp[-zap(P,t)]|2) 
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+ COS e{t)exp[-i{ap{P,t) - as{P,t))]\3)} (19b) 

and their corresponding eigenvahies are = fujj^ = ^H^l{t). Here the phase 
5{t) is also a global phase. The phase difference ap{P,t) — as{P,t) in the adia- 
batic eigenvectors is given by 

ap(P,t)-a,(P,i) =0oW-<^iW 

+ {-[a;o2-a;o] + [a;i2-a;i]+ ^ °^ + ' }t, (20) 

and there is a relation between the phase difference and the detunings: 

ap{PJ) - as{P,t) = [Mt) - Mt)] + [^s{P) - Ap{P)]t. 

The phase difference is dependent upon both the frequency offsets {u)02 — ijJq) 
and {u)i2 — wi) and also the momentum P. It could be considered that the term 
{ko + ki)P/M in the phase difference is generated by the Doppler effect. A large 
Doppler effect is not favorable for the decelerating and accelerating processes. 
Assume that the atomic motional state is a wave-packet state with a finite wave- 
packet spreading. Then in order to minimize the effect of the Doppler-effcct term 
{ko + ki)P/M on the STIRAP state transfer in the decelerating and acceler- 
ating processes one must choose suitably the frequency offsets {u)q2 — '^o) and 
{uJi2 — LOi) for the two Raman laser light beams. Now Po and /^Pm are denoted 
as the central position and the effective momentum bandwidth of the atomic 
momentum wave-packet state, respectively, and AP = P— Pq as the deviation of 
the position P of the momentum wave-packet state from the central position Po . 
For simplicity, hereafter assume that the momentum Pq is always much greater 
than APm/2 in the unitary decelerating process, this means that the atom al- 
ways moves along the same direction -\-x in the decelerating process. Suppose 
that the absolute amplitude at the position P of the momentum wave-packet 
state decays exponentially as the deviation AP. A typical example of such 
wave-packet states is Gaussian wave-packet states [25, 31]. Then the amplitude 
of the position P outside the effective momentum region [P] = [Pq — APm/2, 
Po + APm/2] in the momentum wave-packet state is almost zero and the prob- 
ability to find that the atom is not in the effective momentum region [P] is so 
small that it can be negligible. This is just the definition of the effective mo- 
mentum bandwidth APm of the atomic momentum wave-packet state. Then 
it is sufficient to consider only the effective momentum region [P] of the mo- 
mentum wave-packet state when the adiabatic condition is investigated for the 
three-state STIRAP experiments of the decelerating and accelerating processes. 
Here the carrier frequencies ujq and ldi oi the two Raman laser light beams are 
not determined until the adiabatic condition for the STIRAP experiments is 
examined later. 

The unitary decelerating and accelerating processes of the quantum control 
process [22] require that the atom in the initial internal state |(7o) (Isi)) be com- 
pletely transferred to another internal state \gi) (l^o)) and at the same time the 
initial atomic wave-packet motional state be completely converted into another 
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wave-packet motional state. It is well known that in an ideal condition the STI- 
RAP population transfer process can achieve 100% transfer efficiency from one 
atomic internal state to another [3, 4, 15, 17, 18]. The adiabatic state-transfer 
channel for the three-state STIRAP experiments is formed generally through the 
special adiabatic eigenstate \g^{P, t)) of the Hamiltonian H{P, t) [17]. Therefore, 
if the unitary decelerating and accelerating processes want to achieve their main 
goal, they had better make full use of this adiabatic state-transfer channel. If 
now the atom is prepared in the adiabatic eigenstate \g^{P,tQ)) of the Hamilto- 
nian H{P, to) at the initial time to, then according as the adiabatic theorem [25, 
30] the atom will be in the adiabatic eigenstate |g°(P, f)) of the Hamiltonian 
H{P, t) at any instant of time t in the adiabatic process for io <t<tQ + T if the 
adiabatic condition is met, that is, if the time period T is infinitely large or more 
generally if the eigenvectors of the Hamiltonian H{P, t) vary infinitely slowly 
[30]. However, in practice the time period T can not be infinitely large and ro- 
tating of the eigenvectors of the Hamiltonian are not infinitely slow. Thus, the 
ideal adiabatic condition can not met perfectly in practice. In fact, the quan- 
tum control process does not allow the time period T of the adiabatic process 
to take an infinitely large value. Obviously, if the time interval T is taken as 
a finite value but large enough or the eigenvectors of the Hamiltonian ff(P, t) 
rotate sufficiently slowly, then the adiabatic theorem still holds approximately 
and the real adiabatic state at the end time to + ^ of the adiabatic process will 
be very close to the ideal adiabatic eigenstate | g"(P, to + T)). Since the real 
adiabatic state at the end of the adiabatic process is close to the ideal adiabatic 
eigenstate \g^{P,to + T)), one may calculate the real wave-packet state of the 
atom at the end of the STIRAP-bascd unitary decelerating or accelerating pro- 
cess with the help of the ideal adiabatic eigenstate \g^{P,to + T)). Then this 
simplifies greatly the calculation for the real wave-packet state of the atom at 
the end of the decelerating or accelerating process, although such a calculation 
could generate an error for the real wave-packet state of the atom. Similarly, if 
the real transfer efficiency for the STIRAP process is calculated with the help 
of the ideal adiabatic eigenstate \g^{P, to + T)) instead of the real adiabatic state 
at the final time of the decelerating or accelerating process, then there exists 
certainly an error for the real transfer efficiency. However, these errors may be 
estimated. In fact, if one finds the real adiabatic condition for the STIRAP 
process, one may estimate these errors, as shown below. Therefore, for the STI- 
RAP process satisfying the real adiabatic condition there is a simple scheme to 
calculate the real transfer efficiency and the real wave-packet state of the atom 
at the end of the STIRAP process: one may use the ideal adiabatic eigenstate 
\g°{P,to + T)) of the Hamiltonian II{P,to + T) to simplify the calculation of 
the real transfer efficiency and the real wave-packet state of the atom, then 
evaluate the generated errors and control these errors to be within the desired 
upper bound by setting suitably the experimental parameters for the STIRAP 
process. This is really the procedure to design the STIRAP pulse sequences for 
the unitary decelerating and accelerating processes. Thus, the problem to be 
solved below is how to design the STIRAP pulse sequence such that in the real 
adiabatic condition the final adiabatic state for the real adiabatic process is still 
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very close to the ideal adiabatic eigenstate \g^{P, to+T)) even if the time period 
T takes a finite value. 

In the STIRAP-based decelerating and accelerating processes the adiabatic 
evolution process of the atom could occur simultaneously and in a parallel 
form in these three-state subspaces {|P + fifco)|so), |-P)|e), \P — fiki)\gi)} for 
all the given momentum P within the effective momentum region [P] if the 
adiabatic condition is met in the effective momentum region [P]. It is suffi- 
cient to examine the adiabatic evolution process of the atom in a three-state 
subspace {\P + fiko)\go), |-P)|e), \P — fik-i)\gi)} with a given momentum P of 
the effective momentum region [P] to set up the adiabatic condition for the 
STIRAP experiments of the decelerating and accelerating processes. Obvi- 
ously, any atomic state in the three-state subspace during the adiabatic evo- 
lution process can be expanded in terms of the three basis vectors of the sub- 
space. If now the three basis vectors of the three-state subspace {\P + hko)\go), 
|P)|e), \P — Hki)\gi)} are taken as the three orthonormal adiabatic eigenvec- 
tors {\g^{P,t)), \g^{P,t))} of the Hamiltonian H{P,t), then the atomic three- 
state vector |$(-P,t)) = {Ao{P,t), Ai{P,t), A2{P,t))'^ in the Schrodin ger equa- 
tion (16) at any instant of time t in the adiabatic evolution process may be 
expanded as [4, 5, 10, 12, 15, 17, 18, 19, 20, 25] 



By substituting the adiabatic cigcnstatcs of Eqs. (19a) and (19b) into the state 
|$(P, t)) of Eq. (21) one can find that the coefficients {ao{P, t), a±{P, t)} in Eq. 
(21) are related to those {Ai{P, t), / = 0, 1, 2} in Eq. (16) by 




to 




to 



(21) 



Ao{P,t) = ao(P,i)cos6'(t)exp[-i7(i)] 





(22a) 





(22b) 



A2{P,t) = exp[-i(ap(P,i) - as{P,t))]{-ao{P,t)sm0{t)exp[-i'y{t)] 
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1 /■* 

H — j=a- {P, t) cos 6{t) exp [-i6{t)]exp[-i / dt'n{t')]}. (22c) 

v2 J to 

Thus, the atomic three-state vector \^{P, t)) may be determined if one knows the 
coefficients ao{P,t) and a±{P,t). Assume that at the initial time to the atom 
is in the adiabatic cigcnstatc |g''(P, to))- Then the adiabatic theorem shows 
that the atom is kept in the adiabatic eigenstate \g'^{P,t)) of the Hamiltonian 
H{P, t) at any instant of time t during the adiabatic process if the ideal adiabatic 
condition is met. Then it follows from the atomic state |$(P, t)) of Eq. (21) that 
if the ideal adiabatic condition is met, the time-dependent coefficient \ao{P,t)\ 
should be kept almost unchanged and is very close to unity over the whole 
adiabatic process, while two other coefficients {|a±(P, t)\} should be very close to 
zero. Generally, these coefficients may be evaluated by solving the Schrodinger 
equation (16) that the state \^{P,t)) of Eq. (21) obeys. Inserting the state 
|<I>(P, t)) of Eq. (21) into the Schrodinger equation (16) one obtains a set of the 
three differential equations with the variables ao(P,t) and a±{P,t), 

d /■* 

— ao(f , t) + ioJoiP, t)aoiP, t) + a+{P, t)iJo,+ {P^ *) exp[i J dt'fl{t')] 



+ a_{P,t)(jQ_{P,t)exp[-i [ dt'n{tf)]=0, (23a) 

J to 

a±{P,t) +i(j±{P,t)a±{P,t) +a^{P,t)uj±^^{P,t)e^[T'2i f dt'n{t')] 

Jto 



to 

+ ao{P,t)u!±fi{P,t)exp[Ti [ dt'n{t')]=0, (23b) 

Jto 

where the coefficients {wk{P,t), u}k,i{P,t)} are defined as 

iu,k{P,t) = {g\P,t)\dg\P,t)/dt) {k = 0,±), 

Uk,i{P,t) = {9HP,t)\-^^W{P,t)) {k^l, k,l = Q,±). 

These coefficients can be obtained explicitly from the adiabatic eigenstates of 
Eqs. (19), 

ujQ{P,t) = -^7(i) - siv? e{t) — [a^{P,t) - as{P,t)], (24a) 

d 1 ^ 

w+{P,t)=^.{P,t) = --5{t)---a^{P,t) 

-\cosH{t)-^^[ap{P,t)-a,{P,t)] (24b) 



and 



too,±{P,t) = -w±fl{P,tr = exp[z(7(i) - Sm 
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X {—e{t)+isme{t)cos9{t) — [ap{P,t) - a«(P,i)]}, (24c) 

- i^cosHit)-^^[apiP,t) - asiP,t)]. (24d) 

Though those coefficients contain the global phases 7(t) and S{t), the final re- 
sults of the adiabatic process obtained from these coefficients will not be affected 
by these global phase factors [18b]. as can be seen later. In order to solve con- 
veniently the equations (23) it could be better to make variable transformations 
[25, 30]: 



ak{P,t) = &fe(P,t)exp[-z f dt'ujk{P,t')] {k = 0,±). (25) 

J to 

Then with the new variables {bk{P,t)} these equations (23) are changed to 



^boiP,t) = ^^p[^(7fa)-^fa))] e(P,t)* 

C/t Y 2 

X {b+{P,t)eyip[i f dt'n+{P,t')]+b-iP,t)e^j>[-i f dt'n_{P,t')]}, (26a) 

J to J to 

§-^b±{P,t) = -i^b^{P,t)T{P,t)eMTi J^' dt'2n{t')] 



exp[-i(7(to) - 6{to))] , 



bo{P,t)e{P,t)exp[Ti [ dt'n±{P,if)]. (26b) 
v2 J to 

The coefficients in Eqs. (26) are obtained from those in Eq. (24): 
n±{P,t) = n{t) ± {[as{P,t) - ^ap{P,t)] sin" 0{t) 

+ [ap(P,t) - ^a,{P,t)]cosH{t)}, (27a) 

Q{P,t) = -e{t)+i^sm2e{t)[ap{P,t) - d«(P,i)], (27b) 

r(P, t) = sin^ 9{t)ap{P, t) + cos^ 9{t)as{P, t), (27c) 

where e{t) = ^0{t), a^iP^t) = §-^as{P,t), and so on. The equations (26) 
are the basic differential equations to describe completely the STIRAP-based 
decelerating and accelerating processes. The set of basic equations (26) is a 
generalization of the basic equations to describe the conventional three-state 
STIRAP experiments [17, 16, 18] when the effect of a momentum distribution 
is taken into account on the STIRAP experiments. The basic equations (26) 
may be used to set up a general adiabatic condition for the three-state STIRAP- 
based decelerating and accelerating processes. 
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The basic differential equations (26) may be rewritten in the matrix form 

i^B{P,t) = MiP,t)BiP,t). (28) 

Here the normahzation three-state vector B{P, t) is defined as {bo{P, t), b+{P, t), 
b-{P,t))'^ and the 3 x 3— dimensional hermitian Hamiltonian M{P,t) is given 

by 

Mi2 Mi3 
Mi*2 M23 



Mi*3 M2*3 



M{P,t): 

in which the matrix elements {My} are defined by 

M12 = exp[i(7(to) - 6{to))MP, ty exp[i / dt'Q+{P, t% 
v2 J to 

M,s = ^cxp[iij{to) - S{to))]eiP,ty exp[-i f dt'fl_{P,t% 
v2 J to 

M23 = ^r(P,i)exp[-j / dt'2n{t')]. 
2 Jto 

Notice that the momentum P in the hermitian Hamiltonian M{P, t) is a pa- 
rameter instead of an operator. The three-state Schrodinger equation (28) may 
have the formal solution: 

B{P,t) = Te^{-i [ dt'M{P,t')}B{P,to) 

Jto 

= {! + {-)[ dtiM{P,ti) + {\f f r dtidt2M{P,ti)M{P,t2) + ...}B{P,to). 

* Jto * Jto Jto 

(29) 

The Dyson series solution (29) may be useful to set up a general adiabatic 
condition for the STIRAP-based decelerating and accelerating processes. The 
detailed discussion will appear in the section seven of the paper. 

4. The STIRAP state-transfer process in the ideal adiabatic condi- 
tion 

Consider first the special case: the ideal adiabatic condition. The ideal 
adiabatic condition usually means that the time interval T of the adiabatic 
process is infinitely large or the adiabatic eigenstates of the Hamiltonian (17) 
rotate infinitely slowly. Here the ideal adiabatic condition means that for any 
instant of time of the adiabatic process the integrations of the basic differential 
equations (26) approach zero for any given momentum P within the effective 
momentum region [P] of the atomic wave-packet motional state. The ideal 
adiabatic condition may be expressed as 



/ 

Jtr 



* d 

dt'[—bi{P,t')\ ^0, l = 0,±; to<t<to+T; Pe [P]. (30) 

to 
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The ideal adiabatic condition (30) shows that the basic equations (26) have 
the solution bi{P,t) bi{P,to) {I = 0,±) for any instant of time t of the 
adiabatic process and any given momentum P within the effective momentum 
region [P]. If at the initial time the atom is prepared in the adiabatic cigen- 
state |(?°(P, to))) which is the initial atomic state |$(P, to)) of Eq. (21) with 
ao(P, to) = 1 and a±(P, to) = 0, then according to the adiabatic theorem [25, 
30] one should find that at any instant of time t of the adiabatic process the 
atom is in the adiabatic eigenstate \g'^{P, t)) of the Hamiltonian H{P, t). In fact, 
by the equations (25) and (26) one can obtain |a;(P, t)| — > |ai(P, to)| (/ = 0,±) 
in the ideal adiabatic condition (30). Hence |ao(P, t)| |ao(P, to)| = 1 and 
•^±(-^1^) (i±{P,to) = for the atomic state |$(P, t)) at the instant of time 
t of the ideal adiabatic process, while this atomic state |$(P, t)) is just the 
adiabatic eigenstate \g'^{P,t)) of the Hamiltonian H{P,t), as can be seen from 
Eq. (21). The basic equations (26) show that the ideal adiabatic condition (30) 
could be achieved if the coefficients ©(P, t) and r(P, t) approach zero sufficiently 
and the Rabi frequencies ri(t) and il±{P,t) are sufficiently large. The Rabi fre- 
quencies r2(t) and n±(P, t) can have a great effect on the adiabatic condition 
of the STIRAP experiments. This point is very important for the quantum 
control process, since a sufficiently long time period T of the adiabatic process 
is not accepted for the quantum control process. Then one may likely make the 
adiabatic condition to be met by setting suitably the Rabi frequencies of the 
Raman laser light beams, although the adiabatic condition is met visually by 
making the time interval T sufficiently large. A general adiabatic condition will 
be discussed in detail later. 

In order to optimize the adiabatic condition the carrier frequencies uo and 
(jji of the Raman laser light beams could be chosen suitably such that 

MP^t) = [ ^^^° + ^^° - (wo2 - u;o)]t + Mt) = ^kot + Mt) (31) 

and 

MP, t) = [ ~^^2M ~ ^'^'2 ~ + "^i^*^ = t'''* + '^i^*^ ^^^^ 

where the momentum difference AP = P — Po and the momentum Po is just the 
central position of the effective momcntiim region [P] of a general momentum 
wave-packet state. When the carrier frequencies are chosen according to Eq. 
(31) and (32), the maximum momentum difference value within the effective 

momentum region [P] is minimum. These equations (31) and (32) can be sat- 
isfied when the carrier frequencies {oJi} {I = 0, 1) are determined from the two 
equations: 

^ - {1V02 - Wo) + Co = -^ko (33a) 

and 

_i^(^,2-a;i)+ci = ^fci (33b) 
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where the wave number ki = uJi/c, cq = Aij(Pq), ci = As{Pq), and ■jj,(f)i{t) = 
Ci + ■^(fi{t). On the other hand, if the carrier frequencies are given in advance, 
then the detunings cq = Ap(Po) and Ci = As(Po) may also be determined from 
these two equations (33) , respectively. In the decelerating and accelerating 
processes the momentum value Pq is generally different for each basic STIRAP- 
based decelerating or accelerating process. Then one may adjust the carrier 
frequencies {wi} or the detunings {c(} so that the equations (33) can be satisfied. 
Now using the phase ap{P, t) of Eq. (31) and as{P, t) of Eq. (32) one can rewrite 
the phase difference of Eq. (20) in the simple form 

AP 

ap{P,t) - a,{P,t) = —{ko + h)t + Mt)-'fii{t)- (34) 

If in the STIRAP experiments at the initial time to the Rabi frequency ils{to) 
of the Stokes pulse is much greater than the one ^lp{to) of the pumping pulse, 
then at the initial time the mixing angle 6{t) in the adiabatic eigenstates of Eqs. 
(19) satisfies [15, 16, 17, 18] 

cos e{to) = - 1 or e{to) - 0. (35) 

In practice both the Rabi frequencies ^s{to) and ^p{to) usually could be rela- 
tively small at the initial time, but the Rabi frequency O,p{to) of the pumping 
pulse is much less than r2,s(to) of the Stokes pulse. The initial mixing angle 
9{to) leads the initial adiabatic eigenstates of the Hamiltonian H{P,to) of 
Eq. (17) to the asymptotic forms 

|5°(P,io))^exp[-n(to)]|l), (36a) 

\g^{P,to)) ^ ^exp[-iS{to)]{Tcxp[-iap{P,to)]\2) 

+ cxp[-i{ap{P,h) - a,{P,to))]m. (36b) 

Notice that the three-state basis vectors |1), |2), and |3) stand for the basis vec- 
tors \P + tiko)\go), \P)\e), and \P — tiki)\gi), respectively. The adiabatic eigen- 
state \g°{P,to)) of the Hamiltonian H{P,tQ) is much simpler than \g^{P,to)), 
the latter is more complex in that their expansion coefficients are dependent 
on the momentum. In general, at the initial time an atomic system may be 
prepared more easily in the adiabatic eigenstate |g*'(P, to)) up to a global phase 
factor. An important example is that at the initial time the atom is completely 
in the ground internal state l^o) and in the Gaussian wave-packet motional state 
or in a superposition of the momentum states. Now consider that the initial 
state of the atom is prepared in the superposition state |4'(x, r, to)} given by Eq. 
(12a). By comparing Eq. (12) with Eq. (12a) with P' = P + hko one sees that 
the coefficients of the state |4'(x,r, to)) of Eq. (12) are given by Ao(-P, to) = 1, 
^i(P, to) = 0, and A2{P,to) = at the initial time to. Then it follows from 
Eqs. (15) that the coefficients {Ai{P,to)} can be obtained from {Ai{P,to)} ■ 

Ao{P,to) = exp[^{^-^±^ + Eo)to], A,{P,to) = A^iP,^) = 0. (37) 
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These coefScients associated with the initial mixing angle 0{to) = (here 0{to) 
is so small that it can be taken as zero without losing generality) are inserted 
into Eqs. (22) one obtains, by solving the equations (22), 

ao{P,to) = exp[ij{to)]exp[^{^^^^^^ + Eo)to], a+{P,to) = a-{P,to) = 0. 

(38) 

Indeed, at the initial time the atom is completely in the adiabatic eigenstate 

|.g''(P, to)) of the Hamiltonian H{P,tci) of (17) because both the coefficients 
a+{P,to) and a_(P, to) are zero in the initial atomic state |<i>(P, to))- It follows 
from Eq. (25) that at the initial time to the coefEcient ai{P,to) = bi{P,to) 
{I = 0,±). Then at the initial time to the variables {6;(P, to)} of the basic 
equations (26) can be obtained from those coefficients of Eq. (38), 

bo{P, to) = exp[n(to)] exp[^( + ^o)to], 6+(P,to) = 6_(P,to) = 0. 

(39) 

These coefficients {bi{P,to)} of Eq. (39) provide the basic equations (26) with 
the initial values. 

Now the atom with the initial wave-packet state |5'(x,r, to)) of Eq. (12a) 
undergoes the STIRAP-based decelerating process (11) in the ideal adiabatic 
condition (30). Then one can calculate the atomic state at the end time tf = 
to + T of the ideal adiabatic process (11). This atomic state |4'(a;, r,t)) is 
still given by Eq. (12), but the coefficients in Eq. (12) need to be calculated 
explicitly. By integrating the basic equations (26) and using the ideal adiabatic 
condition (30) the solution to the basic equations (26) is given by 

bi{P,t) = bi{P,to), l = 0,±; to<t<to+T; P e [P]. (40) 

Here the initial values {6;(P, to)} are given by Eq. (39). Once the coefficients 
{bi{P,t)} are obtained from Eqs. (40) one may use equation (25) to calculate 
the coefficients {ai{P,t)}: 

ao{P, t) = bo{P, to) exp[i(7(t) - 7(to))] 

X exp{i J dt'sin^ e{t')-^[ap{P,t') - a.(P,t')]}, (41a) 

a+(P,t) = a_(P,t) = 0. (41b) 

Here the frequency parameter uJo{P,t) of Eq. (24a) has been used. The fact 
that the coefficients a+(P, t) = a_(P, t) = shows that the atom in the adi- 
abatic eigenstate \g°{P,to)) of the Hamiltonian H{P,to) at the initial time to 
is transferred adiabatically to the adiabatic eigenstate |g"(P, t)) of the Hamil- 
tonian H{P,t) and finally to the desired adiabatic eigenstate |(?°(P, t/)) of the 
Hamiltonian H{P,tf) at the end of the adiabatic process. This is just consistent 
with the prediction of the adiabatic theorem in quantum mechanics [25, 30]. By 
substituting the coefficients {a;(P, t), I = 0,±} of Eqs. (41) into Eqs. (22) one 
obtains 

Ao{P, t) = exp[-i7(t)]ao(P, t) cos e{t), (42a) 
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Ai{P,t)=0, (42b) 

A2{P,t) = -ao(P,t)expH7(i)]expH(aj,(P,t) - a,(P, i))] sin6i(t). (42c) 

The coefficient Ai (P, t) = shows clearly that in the ideal adiabatic condition 
the atomic excited internal state |e) indeed does not appear during the STIRAP 
state-transfer process. Though one has the coefficients {ai{P,t), I = 0,±} of 
Eqs. (41) at hand, it is not sufficient from Eqs. (42) to determine imiqucly the 
coefficients {Ak{P,t)} if one does not know the mixing angle 9{t) at the end 
time tf = to +T of the ideal adiabatic process. Suppose that at the end of the 
ideal adiabatic process the Rabi frequency Qp{t) for the pumping pulse is much 
greater than the one f2s(f) of the Stokes pulse. Then the mixing angle 9{t) at 
the end of the ideal adiabatic process takes the asymptotic form [15, 17] 

Now inserting the mixing angle 0{tf ) = 7r/2 into Eqs. (42) one obtains uniquely 

AoiP,tf) = A,iP,tf) = 0, (43a) 

A2(P,t/) = -ao{P,tf)e^p[-i-f{tf)]exp[-i{ap{P,tf) - a.iPtf))]. (43b) 

Thus, with the aid of Eqs. (15) and (43) the atomic three-state vector {Ao{P, t), 
Ai{P,t), A2{P,t)}'^ at the end time tf of the ideal adiabatic process is deter- 
mined by 

AoiP,tf)=Ai{P,tf) = 0, (44a) 
A2{P, tf) = - exp[-( ^ + Eo)to] exp[--( ^ + Er)tf] 

xexpH(ap(P,t/)-a,(P,f/))]exp{i J ' dt' sin^ 0{t')-^[ap{P,t')-a,iPt')]}. 

(44b) 

The coefficients Ao{P,tf) = Ai{P,tf) = show that at the end of the ideal 
adiabatic process the atom is transferred completely to the ground internal 
state \gi) from the initial internal state \go) by the basic STIRAP decelerating 
process (11). This is just the desired result of the ideal adiabatic process (11). 
By using the phase difference in Eq. (34) one can further express the coefficient 
A2{P,tf) as 

A^iP, tf) = expmtf,to)] cxp[-(L^^ + Eo)to] 

, iAP-hki)^ , 
xexp[--(^^^+Er)i/] 

AP f*^ 
X exp{-i — (fco + kx)[to + / dt' cos^ 0{t')]} (44c) 
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where the global phase factor is given by 



exp[i/3(t/,io)] = - exp[-i((^o(*/) - (*/))] 
X exp{i / ^ dt' sin"^ 0{t')[<po{t') - ifiiit')]}. (45) 

J to 

Once the atomic three-state vector {Ao{P,t),Ai{P,t),A2{P,t)}^ is obtained at 
the end time tj of the basic STIRAP decelerating process (11), through the 
equation (12) one can calculate the motional state of the atom at the end of 
the ideal adiabatic process. In order to calculate conveniently the wave-packet 
motional state of the atom one may use the momentum P' = P+Hko as variable 
to express the three-state vector {Ao{P' , t), Ai{P' , t), A2{P' , t)}'^, and then the 
vector at the end time tf can be determined by 

Ao{P',tf) = A,iP',tf) = 0, (46a) 
A2iP',tf) = eMiP'{tf,to)] ^M^iQ + Eo)to] 

xexp[--( — +E,)tf] 

xexp{-i— — (fco + fci)[to+ / dt' cos^e{t')]}. (46b) 

^ J to 

Note that here the momentum difference AP' still represents the deviation of 
the momentum point P' from the central point of the effective momentum region 
[P] in the initial wave-packet motional state. 

For the STIRAP-based accelerating process (11a) in the ideal adiabatic con- 
dition the three-state vector {Ao{P' , t), Ai{P' , t), A^iP', i)}^ at the end of the 
ideal adiabatic process (11a) should be determined by 

AoiP',tf) = Ai{P',tf)=0, (47a) 
MP',tf) = eM^fi'{tfM]eM^{Q+Eo)to] 

. i .{p' + nko + nkif , „ 1 

2M 

Af" f^^ 
X exp{i— — (fco + fci)[io+ / dt' cos^ e{t')]}. (47b) 

^ Jto 
This is because the propagating directions of the Raman laser light beams in 
the accelerating process are just opposite to those in the decelerating process, 
respectively. 

5. The decelerating and accelerating processes of a Gaussian wave- 
packet state in the ideal adiabatic condition 
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As a typical example, consider that an atom with a Gaussian wave-packet 
state is decelerated by the basic STIRAP-based decelerating sequence (11). For 
simplicity, here consider the simple situation that the STIRAP decelerating 
process (11) satisfies the ideal adiabatic condition (30). The theory developed 
in previous sections can determine the wave-packet motional state of the atom at 
the end of the basic decelerating process (11). Suppose that the initial motional 
state of the atom is a standard Gaussian wave-packet state [25, 31] in one- 
dimensional coordinate space: 



*o(a;,io) = exp(iyo)[- 



(Aa:)^i/4 1 



27r y(Ax)2+^(f^) 



2 



Here the complex linewidth of the Gaussian wave-packet state \l/o(a;,<o) is de- 
fined as 

W{To) = {Axr+i{^). 
The probability density of the state ^'o(x, to) is given by 



I -r / m2 1 /l 1 r 1 (X — Zo)^ 

By comparing |^'o(x,to)P with the standard Gaussian function G{x) = [ey/Tr\~^ 
X exp[— (a; — a;o)^/e^] one sees that the center-of-mass position and the wave- 
packet spreading of the state ^o(a;, to) are zq and Eq = y^2[(Aa;)^ + (217^^)^]' 
respectively. If the atom is a free particle, the Gaussian wave-packet state 
^o{x,to) has an explicit physical meaning that the atom with the Gaussian 
wave-packet state moves along the direction +x with the velocity po/M. One 
may expand the coordinate-space Gaussian wave-packet state 4*0 (x, to) in terms 
of the momentum eigenstates {\p)} of Eq. (7), 

*o(x,io) =^Po(P,^o)b). (49) 
p 

The expansion coefficient or the amplitude Po{p,to) is determined by 

Po{P,to) = j dx'^o{x,to)^p{x)* (49a) 

where ipp{x) = exp{ipx/ti) is just the momentum eigenstate |V'p(a;)) or \p) 
of Eq. (7). By substituting the wave-packet state ^'0(2;, to) of Eq. (48) and 
the momentum eigenstate \tpp{x)) of Eq. (7) into the amplitude Po{p, to) of Eq. 
(49a) one obtains, by a complex calculation, 

Po{p,to) = exp(i^o)[^^]^/'exp{-(Ax)2(^)2} 
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xexp{-z(— )(— ^) }exp{-t(— ^)2;o}. (50) 

The amplitude Po{p,to) is really the Fourier transform state of the coordinate- 
space Gaussian wave-packet state ^'o(a;,to)- It also has a Gaussian shape and 
this can be seen more clearly from its absolute square: 

IPoipMf = [^^]^/^exp{-2(Aa=)2(^)n. (50a) 

This is a standard Gaussian function with the propagation-vector variable k = 
p/h. Therefore, it satisfies the normalization, 

/+00 POO 
dk\po{k,to)? = dx\<l!o{x,to)? = 1- (51) 
-CX) J —oo 

The center-of-mass position of the Gaussian function is po for the momen- 
tum variable or po/H for the propagation- vector variable. Tlic^ wavc^-packet 
spreading of the Gaussian function for the momentum p is determined by 
{Ap) = fi/[\/2(Ax)] and for the propagation- vector variable k is given by Ak = 
[\/2{Ax)]^^ . The state PQ{p,to) is called the momentum-space Gaussian wave- 
packet state of the atom. There is a Gaussian decay factor exp[— (Aa;)^( ^"^^° )^] 
in the momentum-space wave-packet state Po{p,to), which decides the effec- 
tive momentum region [P] of the momentum wave-packet state. Obviously, the 
probability density |Po(f^>^o)P approaches zero rapidly (exponentially) when the 
momentum p takes a value such that the deviation |p — po I is greater than the 
wave-packet spreading (Ap). Equation (49) really shows that the coordinate- 
space Gaussian wave-packet state ^'o(a;,io) can be expanded in terms of the 
momentum-space Gaussian wave-packet states {po{p,to)}- 

Consider the superposition of the momentum wave-packet states {Po(P> *o)} : 

Mx,to)= J2 Po{P,to)\p), (52) 

|p-Po|<APm/2 

where APm is just the bandwidth of the effective momentum region [P] of the 
momentum wave-packet state p^{p,tf)). When the bandwidth APm oo, the 
superposition state $0(2^,^0) is just the Gaussian wave-packet state ^'o(a;,to), 
as can be seen from Eqs. (49) and (52). The deviation of the state $0(2;, io) 
from the state '^q{x, t^) may be measured by the probability 

P{^o{x,to)-Mx,to)} = \ Yl Po{P,to)\p)\^ 

\p-Po\>APm/2 

poo 

|2 



/•oo 

! / dk\po{k,to)f 



2 

)[APM/2+vo\/n 

where the momentum eigenstates \p) of Eq. (7) and their orthonormal relations 
and the momentum wave-packet state Pq{p, to) of Eq. (50) have been used. Now 
using the probability density \Po{p, to)|^ of Eq. (50a) one has 

2 

P{^o{x,to) - Mx,to)} = ^ dyeM-y^) (53) 
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where the lower integral limit is 

{APm){Ax) 

The probability P{\E'o(a;, io) - ^o{x,to)} is bounded on by [32] 



The important thing is that the probability P{'ifo{x,to) — ^o{x,to)} decays ex- 
ponentially with the number yf^ . If the bandwidth APm of effective momentum 
region [P] of the momentum wave-packet state Po(P'^o) is chosen such that the 
number j/m >> 1, then the probability P{5'o(a;,to) — ^o{x,to)} is almost zero. 
As a result, the superposition state <I>o(a;,to) is almost equal to the Gaussian 
wave-packet state \l/o(x, to)- 

If an atomic system is in a momentum superposition state which spreads 
from — oo to +00 in momentum space, then it is generally hard to achieve a 
complete STIRAP state transfer in the atomic system by the basic decelerat- 
ing process (11), since the adiabatic condition can not be met over the whole 
momentum space (—00, -l-oo). On the other hand, an almost complete STIRAP 
state transfer could be achieved for an atomic wave-packet motional state with 
a finite wave-packet spreading by the basic decelerating process (11). This can 
be illustrated through the momentum wave-packet state $o(a^,to) of Eq. (52). 
The state $0(2;, io) is also a superposition of the nionicntum eigenstates of Eq. 
(7). All the momentum components of the state $0(2;, io) are within the ef- 
fective momentum region [P] = [po — APm/2, Po + APm/ 2], as can be seen 
in Eq. (52). When the bandwidth APm of the effective momentum region 
[P] is finite and satisfies po — APm /2 > Hko + tiki , a complete STIRAP state 
transfer could be achieved for the state ^o{x, to) within the effective momentum 
region [P] by the STIRAP decelerating process (11) if the ideal adiabatic condi- 
tion (30) is met within the effective momentum region [P] for the decelerating 
process (11). Now using the same STIRAP pulse sequence (11) one can make 
an almost complete STIRAP state transfer for the Gaussian wavc-packcit state 
^o{x,to) as the state 'i!o{x,to) is almost equal to the state ^o{x,to) when the 
number yM » 1- Hereafter suppose that the number ?/m >> 1 so that the 
Gaussian wave-packet state ^"0(2;, to) can be replaced with the state <I'o(a;,to) 
and vice versa without generating a significant error in evaluating the unitary 
decelerating and accelerating processes. 

Now suppose that at the initial time to the atom is in the Gaussian wave- 
packet motional state \E'o(a;,to) of (48) and in the internal state \go). Then the 
total product state of the atom at the initial time is given by 

*o(a;,r,to) = *o(a;, to)|c/o> =Y1po{p,^o)\p)\9o)- (54) 

p 

By comparing the product state '^o{x,r,to) with that state of Eq. (12a) one 
can see that the amplitude Po(P>*o) in the state ^'o(a;,r, to) just corresponds to 
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the amplitude p{P') in the state of Eq. (12a). This means that at the initial 
time to the probability to find the atom in the three-state subspace {|p)|(?o); 
\p — fiko)\e), \p — Hko — fik-i)\g-i)} is given by |po(P)io)|^- Note that during the 
STIRAP decelerating process (11) this probability is not changed with time. 
Obviously, the coefficients Ao{p,to) = 1, Ai{p,to) = 0, and A2{p,to) = at the 
initial time, as can be seen in Eq. (54). Now the initial atomic product state 
^q{x, r, to) of (54) undergoes the basic STIRAP decelerating process (11). Then 
at the end time tf = to + T oi the decelerating process (11) the total product 
state of the atom is given by Eq. (12), 

*o(a;,r,i/) = ^ pQ{p,to){Ao{p,tf)\p)\go) 



+ Ai{p,tf)\p - nko)\e) + A2ip,tf)\p - hko - }iki)\gi)} (55) 

where in the ideal adiabatic condition (30) the coefficients {A2{p, tf), I = 0,1,2} 
are given by Eqs. (46) with the parameter settings P' = p and AP' = p — po. 
Though in Eq. (55) the sum for the momentum p runs over only the effective 
momentum region [P], it will not generate a significant error if the sum really 
runs over the whole momentum region (— oc, +oo), as pointed out before. Since 
in the ideal adiabatic condition (30) the coefficients Ao{p,tf) = A-i{p,tf) = 0, 
the total product state (55) is reduced to the simple form 

'^oix, r,tf) =^Po{p,to)A2{p,tf)\p- hko - hki)\gi). (56) 
p 

The product state '^o(x,r,tf) shows that in the ideal adiabatic condition (30) 
at the end of the decelerating process (11) the atom is completely in the ground 
internal state j^i) and also in the wave-packet motional state: 

^o{x,tf) = Y^Po{p,to)A2{p,tf)\p-hko - hki). (57) 



The initial product state of Eq. (54) and the final product state of Eq. (56) 
show that indeed, the atom is transferred completely from the initial internal 
state l^fo) and the Gaussian wave-packet motional state ^'o(a;, to) of (48) to the 
final internal state \gi) and the motional state 'I'o(a;, tj) of (57), respectively, by 
the STIRAP decelerating sequence (11) in the ideal adiabatic condition (30). It 
can turn out that the motional state ^'o(a;, t/) of Eq. (57) is still a Gaussian 
wave-packet state. By the new variable q = {p — po)/h the coefhcient A2{p,tf) 
of Eq. (46b) with P' = p and AP' = p — po and the amplitude Po{p, to) of Eq. 
(50) are respectively rewritten as 

2 

A2{p, tf) = e^pmtf,to)] exp[^(|^ + Eo)to] 

i {po - hko - hki)^ , ^^ 1 r ■2M*/-*o)i 
2M + E^)tf]eM-^<l ^^^^l 
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X exp[-tq — [tf - to)] exp{-tq — J dt cos^ e{t )} 

and 

Po{p,to) = exp(i<po)[^^^]i/4exp{-g2[(Ax)2 + j(^)]}exp(-ig^o)- 

TT ZIVl 

Inserting the two coefficients and the momentum eigenstate \p — Hko — Hki) of 
Eq. (7) into Eq. (57) one obtains 

M^,tf) = exp(i<^o)[^^^]'/'exp[i/?(tj,io)]exp[|(|| + £o)to] 

X exp[-^(^^^— ^^^^ + E,)tf] exp[i(po - hko - nk{)x/n] 
1 f°° 

X / dqcyi'p{—aq^ + bq) 



where the sum has been replaced with the integral as the momentum p is 
continuous for a free atom, and the coefficients a and b are given by 



a = {^xY + i 



2 , . KTq + tf - to) 



2M 



{po-hko-hki) h{ko + ki) f*f . 2 or^'M 
b = i[x-zo — {tf - to) — / dt cos e{t )]. 



to 



The Gaussian integral in the state 'ifo{x,tf) can be calculated by 



f 

J — ( 



/ TT b 

dqexjp{—aq^ +bq) = — oxp(— ). (58) 



Now by a simple calculation one obtains the final state; 
*o(a;,i/) = exp[iip^{tf,to)][ 



(Ax)^li/4 

2^ ^ \/(Ax)2+i^Mi^ 



2M 

^2 



xexp{-i J^}} — -^}exp{ipix/h}. (59) 

Indeed, the motional state \E'o(a;,f/) is also a Gaussian wave- packet state just 
like the initial motional state ^o{x,to) of Eq. (48). Here the Gaussian wave- 
packet state ^o{x,tf) has the center-of-mass position 

{po - hko - Hki) ,x , fi{ko + ki) 2 an.i\ ccn\ 

zi=zo + — {tf - to) + — / dt cos 6{t ), (60) 



to 



the wave-packet spreading 



35 



the mean motional momentum 



Pi=Po- fiko - Hki, 

and the global phase factor 

It is interesting to compare the final motional state ^"0(2;,*/) with the initial 
state ^o{x, to) of the decelerating process (11). It needs only three parameters to 
characterize completely a Gaussian wave-packet motional state of a free particle: 
the center-of-mass position, the mean motional momentum (or velocity), and the 
complex lincwidth. Hero the wave-packet spreading is determined completely by 
the complex linewidth. For the first point, the atom is decelerated by hko/M + 
Kki/M by the STIRAP pulse sequence (11) as expected, because the average 
motional momentum pi = po ^ Hkft — Hki of the final state '9o{x,tf) is smaller 
than po of the initial state '^o{x,to) and their difference is (fifco + tiki). Here 
suppose that the initial velocity for the moving atom po/M » ti{ko + ki)/M. 
Note that fifco /M and hki /M are the atomic recoil velocities in the two Raman 
laser light beams with the wave numbers fco and ki , respectively. For the second 
point, the atom moves a distance (zi — zq) along the direction +x during the 
decelerating process. If a free atom moved along the direction +x with the 
velocities po/M and Pi/M, respectively, then in the time interval T = tf — to 
the atom would move distances equal to T x po/M and T x pi/M, respectively. 
Here the velocities po/M and pi/M are the atomic moving velocities before and 
after the decelerating process, respectively. One can expect that the distance 
{zi — zo) should lie in between the distances T x p\/M and T x po/M, that 
is, T X pi/M < zi — zq < T X Po/M. Indeed, the equation (60) shows this 
point. For the third point, the wave-packet spreading of the atom at the end 
of the decelerating process is larger than the initial one. If one compares the 
wave-packet spreading of Eq. (61) with the free- atom wave-packet spreading 
(see section 6 below), one can see that though the atom is irradiated by the 
Raman laser light beams, the wave-packet spreading of the atom during the 
decelerating process is not really affected by the Raman laser light beams and 
is just the same as that one of the atom in the absence of the Raman laser 
light beams. This point is important as the wave-packet spreading of the atom 
after the decelerating (or accelerating) process can be calculated easily. Note 
that the wave-packet spreading of a free atom becomes larger and larger as time 
increases, as can be seen in section 6 below. 

From the experimental viewpoint one does not expect that after the de- 
celerating process (11) the wave-packet spreading of the atomic momentum 
wave-packet state could become larger, because this may make a trouble for the 
design of the STIRAP pulse sequence (11). Fortunately, it can turn out that in 
the ideal adiabatic condition (30) the wave-packet spreading of the momentum 
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wave-packet state is not really affected by the STIRAP pulse sequence (11). 
One can expand the wave-packet motional state of Eq. (59) in terms of the 
momentum eigenstates {\p)} of Eq. (7): ^o{x,tf) = X^pPi(p, t/)b) just like 
the expansion (49), and just like the momentum wave-packet state Pq{p, to) the 
momentum wave- packet state Pi{p, tf) can be calculated from Eq. (49a), 

TT n 

X exp{-i ^(^° +^ - il^f} eM-i^^z.]. (62) 

Indeed, the probability density \Pi(p,tf)\'^ is also a Gaussian function and it 
is really equal to \po{p,to)\^ of the initial state \l/o(a;,to) if the initial momen- 
tum po is replaced with pi. Thus, the wave-packet spreading of the Gaussian 
function \pi{p,tf)\'^ is equal to that one of |/9q(p, io)P and is also given by 
{Ap) = h/[V^{Ax)]. This shows that the effective momentum region [P] of 
the initial momentum wave-packet state Po{p,to) is not changed after the atom 
is decelerated by the STIRAP pulse sequence (11), although the center-of-mass 
position of the momentum wave-packet state is changed to pi = po — hko — hki 
from the initial one po after the atom is decelerated. 

In the quantum control process the halting-qubit atom usually needs to be 
decelerated continuously [22] because each STIRAP pulse sequence (11) usually 
can decelerate the atom only by a small velocity value. As shown above, after 
the STIRAP decelerating process (11) the decelerated atom is in the product 
state \E'o(x, r, i/) of Eq. (56), that is, the atom is in the internal state j^i) and 
the Gaussian wave-packet state ^o{x,tf) of Eq. (59). Now the atom needs to 
be decelerated further by another STIRAP pulse sequence. This basic STIRAP- 
based decelerating process may be expressed as 

|P + iilo)\9i) ^ |P)|e) ^ |P - hh)\go). (63) 

In this decelerating process the atom is changed from the internal state |gi) to 
\go). This is opposite to the previous decelerating process (11). Therefore, the 
experimental parameters for the STIRAP pulse sequence (63) needs to be set 
suitably. Now the two Raman laser light beams in Eq. (9) for the STIRAP 
pulse sequence (63) should have the parameter settings: {ELo{t),kLo,^^Lo) = 
(£'Qi(t), ?0) Wq) and {ELi{t),kLi,coLi) = {Ei2{t),h,coi), where the first Raman 
laser light beam (i?Qj^ (t), ?o, t^o) (the pumping pulse) couples the two internal 
states \gi) and |e) and it propagates along the opposite motional direction to 
the atom, while the second beam {£[2(1) , h , (the Stokes pulse) connects the 
two internal states |go) and |e) and it travels along the motional direction of the 
atom. The atomic three-state subspace {|V'/c(^))} in the STIRAP decelerating 
process (63) should be set by |V'o('')) = Ifl'i), \i>i{r)) = \go), and \tp2{r)) = |e), 
and the transition frequencies between the atomic internal energy levels in the 
decelerating process (63) should be defined by ivq2 = '^12 = {E2 — Ei)/h and 
w\2 = ^02 = {E2 — Eo)/h, respectively. If one makes an exchange Eq <-> Ei in 
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all those results obtained in the previous decelerating process (11), then these 
results can be adopted in the present decelerating process (63). 

Now suppose that at the initial time to the atom is in the product state 
\E'o(.T, r, to) of Eq. (54). The atom first imdergoes the STIRAP decelerating 
process (11) and hence the product state \E'o(a;, r, fo) is completely transferred 
to the product state ^o{x,r,tf) of Eq. (56) at the end time tf = to + T of 
the decelerating process (11). Then the atom undergoes the second STIRAP 
decelerating process (63). Now one wants to calculate the atomic wave- packet 
product state at the end of the second decelerating process (63). At the initial 
time ti = tf of the second decelerating process (63) the atomic product state 
is given by '^i{x,r,ti) = '^i{x,ti)\gi) . Obviously, this product state is just 
the product state of the atom at the end of the first decelerating process (11). 
Thus, the initial motional state '^i{x,ti) is the Gaussian wave-packet state of 
Eq. (59): '^i{x,ti) = \E'o(a;,f/). Then in the ideal adiabatic condition (30) at 
the end time = ti +T of the second decelerating process (63) the atom is 
completely in the product state: 

^i{x,r,t'f) = ^i{x,t'f)\go) (64) 

where the wave-packet motional state ^i{x,t'f) can be calculated by 

^i{x,t'f) = Y,Pi{P,ti)A2{p,t'f)\p - nio - hh). (65) 
p 

Here the amplitude p^{p,ti) is given by Eq. (62) with the time tj = ti, the 
momentum eigenstate \p— Hlo — fill) is still given by Eq. (7), and the coefficient 
A2{p, t'f) with the center-of-mass momentum Pq = Pi and Ap = p—pi is written 
as 

2 

A2{p,t'f) = eMWf,ti)]eM^i^+Ei)ti] 
, i ,(p — Hlo — Hli)^ „ N / 1 

xexp{~i^{lo + h)[ti+ J^' dt'cosHiit')]}, (66) 

where the global phase /3i{t'j,ti) is still calculated by Eq. (45) with the related 
parameter settings such as the mixing angle 9i{t) and the phase modulation 
functions ipio{t) and tpnit) of the present STIRAP pulse sequence (63). Then 
by a complex calculation one can obtain from Eq. (65) the Gaussian wave-packet 
motional state at the end of the decelerating process (63): 



^i{x,t'f) = exp[i^2(t/,*i)][^^]'/' 



y(Ax)2+i^^I^±^^i^g±^^iz!l)) 

^^P^- ^^ ^ .i(To+(t-tc)+(t^-tO) > exp[»P2x/fi] (67) 
(Aa;)2 + 1 ^ — i 
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where the center-of-mass position 02 is given by 



(pi - ///() - hli) 
M 



Hip + h) 

M 




f 



Z2 = Zi + 



iff -tl) + 



dt' cos^ ei{t') 



(68) 



the wave-packet spreading by 



£ = W2(Aa;)2 + 2[ 



fi{To + (ti - to) + jt'^ - h j) 
2M{Ax) 



(69) 



the mean momentum by 



P2=Pi- filo - fill, 



and the global phase factor by 



e-xp[iip2{t'f,ti)] = exp[i(^i(ti,io)]exp[z/3,(t|y,ii)] 




+ Eo)t'f]. 



Here both the basic decelerating sequences (11) and (63) are studied in detail 

as they arc the basic STIRAP-bascd decelerating processes. Any unitary decel- 
erating process in the quantum control process [22] may be constructed with a 
train of these two basic decelerating sequences. 

When the atom is in the Gaussian wave-packet motional state 4'o(.x,io) of 
Eq. (48) at the initial time to, the complex linewidth of the motional state 
is W{To) = (Ax)'^ + inTo/{2M). After the first basic STIRAP decelerating 
process (11) the atom is in the Gaussian wave-packet motional state 5'o(.x,ti) 
{ti = to + T) of Eq. (59) and the state has the complex linewidth W{ti — 
to + To) = (Ax)^ + ih{To + ii - to)/{2M). Then after the second basic STIRAP 
decelerating process (63) the atom is in the Gaussian wave-packet motional state 
\l/i(a;,tj) (<j = ti +T) of Eq. (67) and the complex linewidth of the state is 
W{1^f -to + To) = (Ax)'^ + in[To + (ti - to) + (t'f - ti)]/{2M). Thus, one can 
sec that the real part of the complex linewidth of the Gaussian wave-packet 
motional state of the atom keeps unchanged during these decelerating processes 
(11) and (63), while the imaginary part increases linearly with the time periods 
of these decelerating processes. This result is found not only in the decelerating 
processes but also in the accelerating processes and the free-particle motional 
process. 

In a general case a unitary decelerating process may consist of a train of the 
two basic STIRAP decelerating processes (11) and (63). For convenience, here 
each basic STIRAP decelerating process is set to have the same time period 
T = tfi and suppose that at the initial time of the imitary decelerating process 
the atom is in the internal state \go) and has a large motional momentum such 
that the atom still moves along the initial direction +x even after the unitary 
decelerating process. The basic decelerating sequences (11) and (63) may form 
a basic STIRAP decelerating cycle in such a way that first the decelerating 
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sequence (11) and then the sequence (63) is apphed to the decelerated atom. 
The unitary decelerating process may consist of many these basic STIRAP de- 
celerating cycles. Denote C/'*(ll) and U'^{63) as the unitary propagators of the 
basic decelerating processes (11) and (63), respectively. Then a general unitary 
decelerating process may be expressed as 

Un{2N) = f/2^^(63)C/2Vi(ll)-C^2'„(63)f/2Vi(ll)-f/2(63)C/{'(ll) (70a) 

or 

Ud{2N - 1) = U^^_,{n)...Ut{6S)Ut-i{n)...U^{63)Ut{n) (70b) 

where C/2„_i(ll) is the propagator of the (2n — 1)— th basic decelerating unit 

for n = 1, 2, N in the imitary decelerating process Ud{2N) or Ud{2N — 1), 
while U2n{Q3) is the propagator of the 2n— th basic decelerating unit. Each basic 
decelerating unit with an even index 2n in the unitary decelerating process is 
taken as the basic decelerating process (63), while that with an odd index (2n— 1) 
is taken as the basic decelerating process (11). Thus, t^2n-i (-'--'-) = U'^{11) 
and Ui„{63) = ?7<^(63) for n = 1,2, ...,iV. In particular, Ud{0) = E (the 
imit operator), Ud{1) = U'^{^^), and Ud{2) = ?7''(63)?7''(11). The unitary 
decelerating processes J7d(1) and Ud{2) have been investigated in detail in 
the preceding paragraphs. Obviously, the unitary decelerating process Ud{2N) 
consists of 2A'' basic decelerating processes (11) and (63) alternately or N basic 
decelerating cycles, while Ud{2N — 1) consists of N basic decelerating process 
(11) and N — 1 basic decelerating process (63) alternately. 

The time evolution process of the atom in the presence of the unitary decel- 
erating sequence Ud{2N) or Ud{2N — 1) can be calculated exactly in the ideal 
adiabatic condition (30). For the simplest cases N = 1 and 2 the time evolution 
processes of the atom have already calculated in the ideal adiabatic condition 
(30) in the previous paragraphs. In order to calculate the time evolution pro- 
cess of the atom in a general unitary decelerating process one may first set up 
the recursive relation between the two atomic product states at the end of the 
unitary decelerating processes C/£i(2n — 1) and UD{2n) for n = 1,2, N. Since 
the initial internal state of the atom is \go) in both the unitary decelerating pro- 
cesses Ud{2N) and Ud{2N — 1), after the unitary decelerating process UD{2n) 
(or UD{2n — 1)) (1 < n < N) the final internal state of the atom is clearly \go) 
(or l^i)). Then the initial internal states of the atom in the basic decelerating 
processes ?/2n+i(ll) and U2n{Q3) should be \go) and \gi), respectively. It is 
known that at the initial time to the atom is in the Gaussian wave-packet mo- 
tional state '^o{x,to) of Eq. (48) and the product state ^'o(a;,r, to) of Eq. (54). 
It is also known that after the imitary decelerating sequences ?7_d(1) and U]j(2) 
act on the initial product state "^oi^j to) of Eq. (54) the initial motional state 
\l/o(x, ^o) of Eq. (48) is converted into the Gaussian wave-packet motional states 
^o{x,tf) of Eq. (59) and 'i'i{x,t'j:) of Eq. (67), respectively. This means that 
the unitary decelerating sequences Ud{1) and Ud{2) do not change the Gaus- 
sian shape of the atomic motional state. Therefore, it is reasonable to deduce 
that after the unitary decelerating process UD{2n) for n = 0, 1, 2, N the atom 
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is in the Gaussian wave-packet motional state: 



2M 



and also in the atomic wave-packet product state: 

*l(a;,r,tL-i) = "^LMn-M = $^p(p,t2n-i)b)bo). (72) 

p 

where ^^-i is the end time of the unitary decelerating process Uoi^n). In an 
analogous way to calculating the amplitude Po{p,to) via the equation (49a) one 
can calculate the amplitude i^n-i) of Eq. (72) from the motional state 
^2n{x,t^2n-i) of Eq. (71). The resuh is 

p(p,tL-i) = exp(i^t)[^i^]V4exp{-(Aa.)^(^^)n 

X eM-ii ^^^'^^'^'^h i^-^neM-ii^^^L}- (73) 

It will prove below that the states ^'2ri(''^i ^2ri-i) of (71) and 'I'2?i(-^i "^i *2ri-i) of 
(72) are indeed the wave-packet motional state and product state of the atom 
at the end of the unitary decelerating process UD{^n), respectively. 

First of all, the product state '^o{x, r, to) of Eq. (54) and the motional state 
\E'o(a;, to) of Eq. (48) are the initial product state and motional state of the atom 
in the presence of the unitary decelerating process Uoi^n) (or Un{^n — 1)), 
respectively. Of course, these two states may also be formally thought of as 
the final wave-packet product state and motional state of the atom after the 
unitary 'decelerating' process Un{0) = E (the unity operator), respectively. 
This means that the atomic wave-packet product state '^q{x, r, fLi) of Eq. (72) 
and the motional state '^^{Xjtti) of Eq. (71) with n = should be equal to 
*o(a;,r,to) of Eq. (54) and ^Q{x,to) of Eq. (48), respectively, 

*^(a;,r,ili) = *o(a;,r,to), ^oix,tU) = ^o{x,to), 

while the momentum wave-packet state p{p, tti) of Eq. (73) thus is just Pq{p, to) 
of Eq. (50). Indeed, these equations (71), (72), and (73) show this point, if in 
Eqs. (71), (72), and (73) one sets the parameters: n = 0, fL^ = to = to, 
(pg — LfQ, Td = To, Zq = zo, Pq = Po, and td = T, where T is the time period 
of the basic decelerating sequence (11). Next, the product state ^i(x,r, t^) of 
Eq. (64) and the motional state '^i{x, t'^) of Eq. (67) with the time t'^ = t\ are 
just the product state ^"2(2;, r, ) of Eq. (72) and the motional state ^'2(a;,tf) 
of Eq. (71) with n = 1 at the end time t\ of the unitary decelerating process 
J7£)(2), respectively, 

^i{x,r,tl) = ^i{,x,r,t'f), *^(x,tf) = *i(x,i'^). 
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This can be confirmed by setting the following parameters in Eqs. (71) and (72): 
ti = tf = h= to+td,tl= ff =ti + td = to + 2td, = <^2(ii, ), 4 = Z2, 
and P2 = P2, One therefore shows that the motional state 'i'2n{x, t^n-i) ^q. 
(71) and the product state ^'fnl^^, iin-i) of Eq. (72) are correct for both the 
unitary decelerating processes Ud{0) (n = 0) and C/£)(2) {n = 1). It will prove 
below that both the motional state "^ini^^^'in-i) of Eq. (71) and the product 
state ^'fnl^^j '^2n-i) of Eq. (72) are also correct for the unitary decelerating 
process Uoi'^n) with n = 0, 1, A*". 

Suppose that the states '^ini^^^n-i) o^ Eq. (71) and \E'2„(a;, r, t|„_;^) of 
Eq. (72) are correct for the imitary decelerating process Uoi^n). It is known 
that the internal state of the atom is l^o) at the end of the unitary decelerating 
process UDC^in + !))■ Then one needs only to prove that the motional state of 
Eq. (71) is also correct for the unitary decelerating process UD('2{n + 1)). The 
propagator of the unitary decelerating process U£){2{n +1)) can be written as 

f/z,(2(n + 1)) = Ut+2ie5)UD{2n+l) = U^^+2iQmL+i{^WDi2n). 



According to the assumption the motional state *2n(^'^2n-i) J^^^ the final 
motional state of the atom when the atom is acted on by the unitary propagator 
Uoi^n). Obviously, the motional state ^'2n(^'^2n) = ^2n(2^' *2n-i) ^^^o the 
initial motional state of the (2n + l)— th basic decelerating process (11) with the 
propagator ?72n+i(ll) ™ the unitary decelerating process /7_D(2(n + 1)), where 
the initial time is denoted as t2n = ^2n-i and there are the recursive relations: 
tU =4 = to, tl=ti = 4 + td, tl=4 = 4 + 2td, and = 4n-i+td {n > 0). 
Then the initial product state of the basic decelerating process [^2n+i(ll) is 
given by *^„(x,r,t^„) = *2n(2^' *2n-i) of Eq. (72). Now the initial product 
state ^2n{^'^'4n) is acted on by the unitary propagator f/2n+i(ll)- Then it 
can turn out that at the end time t^^ = 4^ + td of the basic decelerating process 
^2n+i(ll) the initial motional state ^'2„(-i', 4n) and product state *2n(^' 4n) 
are respectively transferred into the motional state: 
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X '-^^-hjK^^f^^^^SlA^ exp{^P,t,,x/fi} (74) 
and the product state: 

*2„+l(^,^,i2j = *2„+l(^,*2„)l5l) (75) 

where 

Pin+i = PL-fiko-fih, (76) 

zU, = zL + ^td + ^^h+bl dt' cos^ 9{t% {77) 
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xexp{--(^-gi^ + i;i)iL}- (78) 

The computational process from the initial state ^2n(^>^2n-i) the final state 
^2n+i(^!^2n) ^ ^he Same as the previous one from the initial state \l/o(a;,to) of 
(48) to the final state ^'o(x,t/) of (59). There are the relations: 

*2n+i(a^, r, tin) = <+i(ll)*2„(a^. r, 4j = U ^{^n + l)vE'o(x, r, to). 

These relations show that both the states *2n+i(^)'')^2n) ^"^^ *2n+i(2^' ^in) 
are also the product state and the motional state of the atom at the end time 
^2n — ^2n + of the Unitary decelerating process UdC^ti + 1), respectively. 

The atomic product state ^'2n+i(^''^'^2n) of Eq. (75) at the end of the 
unitary decelerating process UdC^ti + 1) will be used to calculate the atomic 
product state at the end of the unitary decelerating process [/£)(2n + 2). This 
computational process is the same as the previous one from the initial state 
*i(a;,ii) = ^o{x,tf) of (59) to the final state ^i(x,tf) of (67). There are the 
relations: 

^in+2{x,r,tl^+r) = Ut+2(Q3)'^in+iix,r,tlJ = UDi2n + 2)^o{x,r,to).. 

These relations show that the atomic product state '^2n+2{x, r, t2n+i) the end 
of the unitary decelerating process UD{2n + 2) can be obtained by applying the 
propagator ?/2ri+2(63) to the atomic product state ^'2^+1(2;, ^in) of Eq. (75). 
For convenient calculation, the atomic motional state '^in+A^^^'^n) of Eq. (74) 
is rewritten as (n' = n + 1) 



*L-i(^,iL-i) = exp(i<pL_i)[^li/4 / ^ 
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1 {^-4n'-i? 

4(Ax)2+zMIkt^^^i)M) 
Then the atomic product state of Eq. (75) can be rewritten as 



^L'-l{x,r,4n'-l) = *2n'-l(a;,iL'-l)lffl) = E'^(^''*2n'-l)|P)lffl) (80) 



where the amplitude p(p, t2n'-i) calculated from the motional state of 

Eq. (79) and is given by 

p(p,4'-i) = exp(»<,_,)[^^]^/^exp{-(Aa:)^( ^ " f "'"^ f} 
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Now the atomic product state *2n'-i(2^' ^' ^2n'-i) of Eq. (80) is applied by the 
unitary propagator C/2„,(63) (n' = n+ 1). Then it can turn out that at the end 
time i2n'-i — ^2n'~i + of the unitary decelerating process Uni'^n + 2) the 
atomic wave-packet motional state takes the form 



2^ ' Y(Ax)2 + iMI}<+|iM 

xexp{-i — }exp{iP2''„,a;/fi}. (82) 

and the atomic product state is 

*t,(a;,r,i^„,_i) = ^tA^An'-M, (83) 

where 

PL'=PL'-i-filo-nh, (84) 

^2n' = Z2n'-1 + -Jf^d H j dt COS ^[{t ), (85) 

exp(«</52n') = exp(i(p^„,_i)exp[i/3;(i^„,_i,t^„,_i)] 

"^P[i^ ^^'2M +^i)4'-i]exp[-^(^^ +^o)tL'-i]- (86) 

Now by comparing the motional state ^2n'(^' ^2n'-i) of Eq. (82) with the 
motional state ^2n{^^ ^2n-i) of Eq. (71) and the product state ^2n'i^' *2n'-i) 
of Eq. (83) with the product state 5'2„(a;, r, t|„_i) of Eq. (72) one can conclude 
by the mathematical principle of induction that the motional state \I'2„(a;, i2n-i) 
of Eq. (71) and the product state ^'2n(^)'')*2n-i) of Eq. (72) are indeed the 
states of the atom at the end of the unitary decelerating process Uoi'^n) for 
n = 0,1,2, N . In an analogous way, one can prove that the product state 
*2n+i(a;,r-,if„+i) of Eq. (75) and the motional state *2n+i(a^' *2n+i) of Eq. 
(74) are the states of the atom at the end of the unitary decelerating process 
J7£,(2n+1) forn = 0,l,2,...,iV-l. 

Now one can prove that the atomic motional momentum ^2^+1 of Eq. (76) 
and P2„/ of Eq. (84) are given by, respectively, 

Pin+i = Po + l)(fifco + hki) - n{nio + Hh), 

P^^, = - n'ihko + nki)- n'{nio + nh). 

It is known that the recursive relations for the atomic motional momentum 

are given by P^n+i = Pin ^ ^^0 — fifci for n = 0, 1, ...,N — 1, and Pl'^, = 
Pin'-i ^ ^ ^'1 '^or n' ~ 1,2, A^, which are obtained from Eqs. (76) and 
(84), respectively. The two recursive relations together can lead directly to the 
two formulae for the atomic motional momentum P2„+i and P^^^, . 
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One also can calculate exactly the time evolution process of the atom in 
the unitary STIRAP-based accelerating process in the ideal adiabatic condition 
(30). There are also two basic STIRAP-based accelerating sequences which 
correspond to the basic decelerating sequences (11) and (63), respectively. One 
of which is already expressed as (11a). The basic accelerating sequence (11a) 
corresponds to the basic decelerating sequence (11). Another may be expressed 
in an intuitive form 



This basic accelerating sequence corresponds to the basic decelerating sequence 
(63). In an analogous way to constructing the unitary decelerating processes 
UdC^N) and Ud{2N — 1) one may build up the unitary accelerating processes 
UaC^N) and UaC^N— 1) out of the basic accelerating sequences (11a) and (63a), 



[/A(2iV) = [/2V(63a)[/2Vi(ll«)-t^2n(63a)[/2Vi(lla)-C/2(63a)[/f(lla) 



C/a(27V - 1) = C/2V-i(lla)-C/2„(63a)[/2Vi(ll«)-C^2 (63a)C/f (11a) (87b) 

where C/2„_i(lla) and C/2n(63a) for n = 1, 2, A'' are the unitary propagators of 
the (2n— 1)— th basic accelerating sequence (11a) and (2n)— th basic accelerating 
sequence (63a), respectively. Here also suppose that the atom is in the internal 
state \go) at the initial time in both the unitary accelerating processes Ua{'2N) 
and Ua{2N - 1). 

The time evolution process of the atom in the unitary accelerating process 
UaC^N) (and UaC^N — 1)) can be calculated exactly in the ideal adiabatic 
condition (30) in a similar way to that one in the imitary decelerating process 
Ud{2N) (and Ud{2N- 1)). Actually, the recursive relations (71)-(78) and (79)- 
(86) of the unitary decelerating process Un{2N) or Un{2N — 1) can be used as 
well for the unitary accelerating process Ua{2N) or Ua{2N — 1) if one makes 
transformations: ko —k^, fci — fc", and 9{t) da{t) in those recursive 
equations (71)-(78) and Iq — > —Iq, Zi — > — Z", and 6i{t) — > 9ia{t) in those recur- 
sive equations (79) (86). As an example, suppose that the initial wave-packet 
motional state for the atom in the unitary accelerating process UA{2n) is given 



P - nio)\gi) ^ \P)\e) ^ \P + fih)\go). 



(63a) 



(87a) 



or 



by 



and the atomic wave-packet product state by 




(88) 



n{x, r, ig) = nix, ig)|5o) = E PiP^ to)\p)\9o). 



(89) 



p 
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Then it can turn out that the momentum wave-packet state p{p, fg) of the 
motional state \E'g(a;,tg) can be written as 

pipXo) = exp(i^S)[^^]V4exp{-(Aa=)2(^^)2} 
IT n 

Now the initial wave-packet product state r, tg) of the atom undergoes the 

unitary accelerating process Ua (2n) . Then it can be proved that at the end of 
the unitary accelerating process [/^(Zn) the atomic wave-packet motional state 
is given by 



*2„(^,iL-i)=exp(i^L)[^]'/' 



27r ' y (Aa;)2+iMZV^ 

4 (Ax)^ + l^^^g '- 

and the atomic wave-packet product state by 

nn{^,r,e2n) = nn{xAn-M (92) 

where the end time of the imitary accelerating process UAi^n) is t|„_i = — 
+ 2nta for n = 0, 1, 2, ...N, the atomic motional momentum is given by 

PL = PS + n{fikS + tik1)+ n{m^ + ft/? ) , (93) 

and the center-of-mass position is determined from these recursive relations: 

%^ia - / dt' cos^ 0„(O, (94a) 

pa f:(ja I ja\ ptg+2kta 

p^^^_n^o+hl dt' cos' eun (94b) 

-P2/S-I = -P2fe-2 + H^O + k^), Pgfe = P^-l + ^(^0 + ^1 )i 

where the initial center-of-mass position and momentum are Zq and P^, respec- 
tively, the index 1 < fc < n, both the basic accelerating sequence (11a) and 
(63a) have the same duration ta, and the global phase factor exp(z(^2n) can also 
be calculated through the recursive relations similar to Eq. (78) and (86). In 
an analogous way, one also can calculate exactly the time evolution process of 
the atom in the unitary accelerating process ?7^(2n — 1) in the ideal adiabatic 
condition (30). 

6. The space- and time-compressing processes based on the uni- 
tary decelerating and accelerating processes 



^2fc-l — ^2fe-2 



^2k — ^2k-l 
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Suppose that in the quantum control process [22] the first unitary decelerat- 
ing process consists of 2nd basic STIRAP decelerating sequences (11) and (63) 
alternately, which is given by Uoi'^nd) of Eq. (70a), and the ideal adiabatic 
condition (30) is met for all these basic decelerating sequences. According to 
the quantum control process the unitary decelerating sequence is applied selec- 
tively in the given spatial region [Dl, D^] in the right-hand potential well of the 
double- well potential field, where Dl and Dij are the left- and right-boundary 
positions of the spatial region in the coordinate axis, respectively. The halting- 
qubit atom can be decelerated by the unitary decelerating sequence J7D(2nd) 
only when the atom enters into the spatial region [Dl, D^]- Thus, the spatial 
region [D^, Dr\ may be called the decelerating spatial region. The decelerating 
spatial region must cover sufficiently the whole wave-packet motional state of 
the atom during the whole imitary decelerating process when the atom is de- 
celerated in the decelerating region. The decelerating region is so wide that for 
the wave-packet motional state of the atom the Raman laser light beams of the 
unitary decelerating sequence Uoi'ind) can be thought of as infinite plane-wave 
electromagnetic fields. Suppose that the halting-qubit atom is in the product 
state ^o{x, r, to) of Eq. (54) (or in the motional state ^o{x, to) of Eq. (48) and 
the internal state Igo)) and in the decelerating region [D^, Dj^] when the unitary 
decelerating sequence Uu{2na) is turned on at the initial time to. Here the spa- 
tial position of an atom is defined as the center-of-mass position of the atomic 
wave-packet motional state. Now the center-of-mass position and wave-packet 
spread of the initial motional state \E'o(a;,fo) are zo and Eq, respectively. Since 
the halting-qubit atom is in the decelerating region [Dl, Dr] at the time to, that 
is, Zo G [Dl, Dr], both the distances (zo — Dl) and {Dr — zq) must be much 
greater than the wave-packet spreading eo, that is, (Dr — zo) > (zq — Dl) >> £o, 
meaning that the decelerating region [Dl, Dr] covers sufficiently the whole ini- 
tial wave-packet motional state ^o{x,to). The halting-qubit atom starts to be 
decelerated by the unitary decelerating process U d (2nd) at the initial time to 
and in the decelerating region [Dl, Dr]. With the help of the recursive relations 
(71) (78) and (79) (86) one can prove that at the end time = + "^n^t^ of 
the unitary decelerating process Uoi'^nd) the motional state of the halting-qubit 
atom is given by 





(95) 



,2 I „• tiiTd+2ndtd) 
' 2M 



and the atomic product state by 



(96) 



where the atomic motional momentum 



is given by 



Pin^ =Po- nd{hko + hk\) - nd{hlo + Hh), 



(97) 
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and the ccnter-of-mass position Z2n^ can be calculated by the recursive relations 
(77) and (85), 



^271-1 = ^271-2 H U-^d -I U / "-^ ^v* )> (98a) 




'2n-2 



zL = 4n-i + + ^^^^ T"" dt' cos' em, (98b) 
pL-1 = pL-2 - {fiko + nh), = Pin-i - {f^k + nh), 



where 1 < n < n^; 4 = zq, = po, Td = U; 4 = h, t^+i =4+ ^d, and 
fj. = t'l_^_-^ for < fc < 2nd — li and the global phase factor exp(z(p2„^) can 
be calculated by Eq. (78) and (86) with the initial phase ip^ = ip^. Both the 
initial atomic product state 5'o(a;,r, io) of Eq. (54) and the final product state 
^2rad(^! '^^ ^2n<i) ^iq. (96) show that before and after the unitary decelerating 
process Uoi^nd) the halting-qubit atom is in the same internal state \go), while 
its initial motional state ^"0(2;, ^o) of Eq. (48) is changed to the motional state 
*2nd(^' *2n<j) (^5) after the unitary decelerating process. The final 

motional state ^'2na(a;, i2nd) ^'i- (^5) has the center-of-mass position ^jn^ 
and the wave-packet spreading 



Though the atom moves a distance {z2nj — ^o) during the unitary decelerating 
process, the center-of-mass position £ [Dl, Dh] as the atom is still in the 
decelerating region [Dl, Dji] at the end of the unitary decelerating process. 
Then both the distances {z2nj — Dl) and {Dn — 22ra<j) must be much greater 
than the wave-packet spreading e, that is, {z2nd~DL) > {DR — ^ind) >> ^- This 
means that the decelerating region [Dl, Db] also covers sufficiently the whole 
final motional state *2n<j('^' ^2nd)- There are two extra constraint conditions on 
the decelerating region [Dl, Dfi\. If the atom has not yet entered into the 
decelerating region when the unitary decelerating sequence is switched on or it 
leaves the decelerating region after the unitary decelerating process is switched 
ofi^, then it will not be aflPected by the unitary decelerating sequence or by next 
unitary decelerating sequences. The two constraint conditions are stated below. 

In the quantum control process [22] the halting-qubit atom may enter into 
the right-hand potential well from the left-hand one in any i— th cycle of the 
quantum program for i = 1,2, ...,TOr. The i— th possible wave-packet motional 
state of the halting-qubit atom is just defined as the atomic motional state 
when the atom enters into the right-hand potential well in the i— th cycle of 
the quantum program. Thus, there is a different time for any possible atomic 
motional state such as the i— th wave-packet motional state to enter into the 
right-hand potential well. If the time period of each cycle of the quantum 
program is AT, then the time difference between the i— th and j— th {i < j) 
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wave-packet motional states to enter into the right-hand potential well is given 
by AT{i,j) = (j — i)AT for i < j and i,j = 1,2, ...,mr- This time difference 
results in a center-of-mass distance in space between the two wave-packet mo- 
tional states. If the halting-qnbit atom moves along the direction +x with the 
velocity po/-^ J then the distance is given by AL{i,j) = (j — i)AT{po/M). Now 
examine two consecutive possible wave-packet motional states: the i—th and 
(i + 1)— th wave-packet motional states. Here for convenience the i—th wave- 
packet motional state is set to the motional state ^'o(a;,fo) of Eq. (48). Then 
at the time to the i—th wave-packet motional state is in the decelerating region 
[Di, Dj^] and its center-of-mass position is zq, while the center-of-mass position 
of the {i + 1)— th wave-packet motional state is clearly [zq — AT{po/M)]. It is 
known that the total duration for the unitary decelerating sequence UDC^nd) 
is 2n(it(i- Then at the end time to + '^natd of the imitary decelerating pro- 
cess the center-of-mass position of the [i + 1)— th wave-packet motional state 
becomes [zo — (AT— 2ndtd) {po /M)] . Obviously, the distance between this center- 
of-mass position and the left-end position of the decelerating region [Dl, Dr\ 
is Dl — [zq — (AT — 2ndtd){po/M)]. Denote ej+i(io + 2ndtd) as the wave- 
packet spreading of the {i + 1)— th wave-packet motional state at the time 
to + 2ndtd- The wave-packet spreading ei+i(to + 2ndtd) may be calculated with 
the help of the z— th wave-packet state ^'o(a;,<o) and the free-particle propaga- 
tor. Then this distance must be much greater than Si_^.i{to + 2ndtd), that is, 
Dl - [zo - (AT - 2ndtd){po/M)] » e,+i{to + 2ndtd), so that the (i + l)-th 
wave-packet motional state is not affected by the unitary decelerating sequence 
UD{2nd) during the whole unitary decelerating process. This is a constraint 
condition on the decelerating region [Dl, Dr]. 

It is known that at the end time to + 2ndtd of the unitary decelerating pro- 
cess Un{2nd) the i—th wave-packet motional state is the state ^2rad(^'*2nd) 
of Eq. (95), which has the center-of-mass position and the motional mo- 
mentum Pina- After the unitary decelerating process the i—th wave- packet 
motional state (i.e., the halting-qubit atom) moves along the direction +x with 
the velocity P2n^/M. Since the atom is usually decelerated greatly by the 
decelerating sequence Uoi^nd) the atomic velocity P2n^/M is much less than 
the original velocity po/M. Obviously, the i—th wave-packet motional state 
moves to the position Z2n^ + (AT — 2ndtd)P2ndl^ when next unitary decel- 
erating sequence starts to work at the time to + AT. Then the distance be- 
tween this position and the right-end position of the decelerating region [Dl, 
Dr] is given by z^„^ + (AT - 2ndtd)Pi„jM - Dr. Denote e,(to + AT) as 
the wave-packet spreading of the i—th wave-packet motional state at the time 
to + AT. The wave-packet spreading e,if.o + AT) can be calculated with the 
help of the motional state '^tudi^^^'in,,) of Eq. (95) and the free-particle prop- 
agator. Then this distance must be much greater than Siito + AT), that is, 
z^^^ + {AT - 2ndtd)P2nJM - Dr » ei{to + AT), so that, from the time 
to + AT on, the i—th wave-packet motional state is no longer affcicted by the 
unitary decelerating sequences. This is another constraint condition on the 
decelerating region [Dl, Dr]. 
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The i—th wave-packet motional state is decelerated from the initial time to 
to the end time to + 2ndtd by the unitary decelerating process Uni'^nd) in the 
decelerating region [Dl, Dr]. It moves a distance Z2n^ — zq along the direction 
+x and spends the time 2ndtd and it is decelerated down to P^,^ /M from the 
initial velocity po/M during the unitary decelerating process. According to the 
quantum control process [22] the {i + 1)— th wave-packet motional state arrives 
at the position zq in the decelerating region [Dl, Dji] at the time to + AT. Then 
the (i + 1)— th wave-packet motional state at the time to + AT is really equal to 
the i—th wave-packet motional state at the time to up to a global phase factor, 
indicating that the {i + 1)— th wave-packet motional state at the time to + AT is 
also equal to the motional state ^'o(a;, to) of Eq. (48) up to a global phase factor. 
Generally, according to the quantum control process each of these rur possible 
wave packet motional states is really equal to the motional state ^0(2;, to) of 
Eq. (48) up to a global phase factor when the wave-packet motional state 
arrives at the same position zq in the decelerating region [D/,, Dn]. Just like 
the I— th wave-packet motional state at the time to the (i + 1)— th wave-packet 
motional state at the time to + AT is decelerated by the unitary decelerating 
process UoC^nd)- It also moves the distance Z2n^ — zo along the direction -\-x and 
spends the time 2ndtd and it is also decelerated down to P2n,,l-^^ from the initial 
velocity po/M during the unitary decelerating process. Generally, each of these 
rUr possible wave-packet motional states moves the same distance Zjnd ~ ^0 along 
the direction +.x and also spends the same time 2ndtd and it is also decelerated 
down to the same velocity P^„^/M from the same initial velocity Po/M during 
the unitary decelerating process. The difference among these rrir possible wave- 
packet motional states is that the starting time is different to decelerate each 
one of these wave-packet motional states by the unitary decelerating process 
UD{2nd). 

In the quantum control process [22] the unitary decelerating sequence is 
used to decelerate the halting-qubit atom so that the center-of-mass distances 
between these rUr possible wave-packet motional states of the atom can be 
narrowed greatly. Thus, the unitary decelerating process Uoi^nd) is really a 
space-compressing process for these possible wave-packet motional states. Since 
each one of these possible wave-packet motional states spends the same 
time 2ndtd in the imitary decelerating process Uiy{2nd), the time difference 
AT{i,j) = [j — i)AT between the i—th and j— th wave-packet motional states 
{i < j; i,j = 1,2, ...,mr) does not change before and after the unitary decel- 
erating process. It is known that each possible wave-packet motional state has 
the initial moving velocity po /M before the unitary decelerating process and the 
final moving velocity P^nJM > after the unitary decelerating process. Here 
the atomic moving velocity P^n^jM can be obtained from Eq. (97), 

Kj^i = [po - udinko + nky) - ndinio + nh)]/M. (99) 

If the number Ud of the unitary decelerating process Uni^nd) is chosen suitably, 
then the velocity P2nj^/M can be much less than the initial one Po/M. Before 
the unitary decelerating process the distance between the i—th and j— th wave- 
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packet motional states is ALo{i, j) = {j — i)AT{po/M), since the velocity is 
Po/M and the time difference is AT(i, j) = {j — i)AT before the unitary de- 
celerating process. After the unitary decelerating process the atomic moving 
velocity is P^,^ /M and the time difference is still AT(i, j) = (j — i)AT. Then 
after the unitary decelerating process the distance between the i—th and j— th 
{i < j) wave-packet motional states is equal to 



where i < j; i,j = l,2,...,mr- Since the velocity (P2„^/M) << po/M, the 
distance AL{i,j) « ALo{i,j), indicating that the spatial region to cover all 
these nir possible wave-packet motional states is greatly compressed after the 
unitary decelerating process. The distance AL{i,j) of Eq. (100) has been 
obtained in the previous paper [22], where the atomic velocity {P2nj/M) is 
denoted as vq after the unitary decelerating process. Then the ratio of the two 
distances AL(i,j) and ALo(«, j) is the space-compressing factor for these wave- 
packet motional states after and before the unitary decelerating process, which 
can be calculated by 



The space-compressing factor is not dependent upon the indices i and j, since 
the time difference AT{i,j) does not change before and after the unitary decel- 
erating process and since all these possible wave-packet motional states have the 
same initial motional momentum pq and also the same final motional momen- 
tum P2nj after each of these possible wave-packet motional states undergoes 
the same unitary decelerating process Uoi'^nd) in the same decelerating region 



Before the unitary accelerating process comes to making a real action on the 
halting-qubit atom, the atom needs to stay in the right-hand potential well for 
a time period to wait for the quantum program running to the end according 
to the quantum control process [22] . The time period during which the halting- 
qubit atom stays in the right-hand potential well is different and dependent 
upon how early the halting-qubit atom enters into the right-hand potential well 
from the left-hand one. When the atom enters into the right-hand potential 
well at an earlier time, it will stay in the right-hand potential well for a longer 
time. Denote that Ts{i) = Tg — [i — 1)AT with the index i — l,2,...,mr is 
the time period during which the atom moves freely along the direction +x 
in the right-hand potential well after the atom is decelerated by the unitary 
decelerating sequence Uoi^nd) and before the atom starts to be accelerated at 
the end time of the quantum program. The index i indicates that the halting- 
qubit atom enters into the right-hand potential well from the left-hand one in 
the i—th cycle of the quantum program. Here suppose that the last unitary 
decelerating process Uoi^nd) is turned off before the quantum program comes 
to the end. The calculation for the time evolution process of the halting-qubit 
atom moving freely during the time period Ts{i) needs to use the free-particle 



AL{i,j) = {j-i)AT{PtjM), 



(100) 



Rg — 



AL{i,j) _ [po - nd{nko + Hki) - nd.{hlo + Hh)] 
ALo{i,j) Po 



(101) 



[Dl, Dr]. 
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unitary propagator. Now the unitary propagator of a free particle is written as 
[25] 

Then the time evolution process of an atom in a free-particle motion with the 
time period T = t' — t can be calculated by 

^{x',t')= J dxG{x',t';x,t)^{x,t). (103) 

It is known that the i— th wave-packet motional state of the atom is given 
by the motional state '^'in^i^^'^^na) '^^ (^5) at the end time ijnd = + 
'^ndtd of the unitary decelerating process. When the wave-packet motional state 
^tridi^T^irid) moves freely along the direction +x for the time period Ts{i) 
from the time f^n-^ to the time ifn^ + Tsii), it will change to another Gaussian 
wave-packet motional state. This Gaussian wave-packet motional state can 
be calculated from the equation (103) by taking the initial state ^'(.x, t) as 
^2nd(^' ^iud) (^^)' using the free-particle propagator G{x' , t'; x, t) of Eq. 

(102), and denoting t! = tg^^ +^s(0 ^^'^ * = ^2nd- ^ complex calculation, in 
which the Gaussian integral (58) has been used, the final Gaussian wave-packet 
state \l/(x',t') can be obtained explicitly, which now is renamed 'iff{x,t2nj + 

^n^, tLd + T.d)) = eMi^Ld) ^-P{-^^^^y2^} 



X 
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xexpj--- — , _i_T^ I exp{^P2r^.,a;/?^}. (104) 

^ [{Ax)"^ I i f^\Td+2ndtd+T,(i)) -^ J fi zud / J V / 

On the other hand, the atomic internal state \go) and the motional momentum 
Piud ^^^P unchanged during the free-particle motion of the atom. Therefore, be- 
fore the unitary accelerating process starts at the end of the quantum program, 
these rrir possible wave-packet motional states are given by 'i'f {x,t2n^ + ^s(i)) 
of Eq. (104) for i = 1, 2, ...,mr and each of them has a different center-of-mass 
position: Z2n^ + {P2n^/M)Ts{i), a different global phase factor: 

eMivLd)^M-i{PLdfTs{i)/{2nM)}, 
mlex linewidth: 



and a different complex linewidth: 



\2 , -rKTci + ^ndtd + Tsd)) 



WiiTd + 2ndtd + T,{i)) = {AxY + i[ ' " ' " ]. (105) 

An important fact is that the imaginary part of the complex linewidth of the 
motional state (x, fj^^ + Tgii)) increases linearly with the time period T's(i), 
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while the real part keeps unchanged. Though each one of these rrir possible 
wave-packet motional states has the same wave-packet spreading and the same 
complex linewidth W{Td+2ndtd) = {Ax)"^ before the free-particle 
motion, as can seen from the state "i'^n^ (-^i ^an^) ^f Eq- (95), each possible wave- 
packet motional state has a larger wave-packet spreading and a different com- 
plex linewidth Wi{Td + 2ndtd + Ts{i)) when the unitary accelerating sequence 
starts to act on the halting-qubit atom at the end of the quantum program. 
Obviously, the first wave-packet motional state *f (a;,t2nd +Ts) has the largest 
wave-packet spreading, while the last motional state {x, f^^^ + Ts{mr)) has 
the least one. These show that the free-particle motion of the halting-qubit 
atom leads to the difference among the wave-packet spreads of these possi- 
ble wave-packet motional states and makes these wave-packet motional states 
broader. This difference may have a significant impact on the quantum control 
process [22]. On the other hand, the free-particle motion of the halting-qubit 
atom does not change the time differences and the distances in space between 
those m,. possible wave-packet motional states. This is because the motional 
momentum Pjnd same for all these rrir possible wave-packet motional 

states and keeps unchanged during the free-particle motion. Thus, the distance 
between the z— th and th [i < j) wavc^-packct motional states is still given by 
AL{i,j) of Eq. (100) and their time difference by AT{i,j) = (j - i)AT. Par- 
ticularly, the distance between two nearest wave-packet motional states is given 
by AT(P^„^/M). Obviously, the halting-qubit atom moves a distance equal 
to {P2n^/M)Ts{i) along the direction +x in the time period Ts{i) of the free- 
particle motion. This distance is dependent upon the index i. The first wave- 
packet motional state ^'f (a;,t2„^ + T,) moves the largest distance {Pina/^'^)'^^ 
which decides mainly the dimensional size of the right-hand potential well, while 
the last wave-packet motional state '^m^i^j'^ina + Ts{inr)) moves the shortest 
distance {P^„JM)[T, - (m^ - 1)AT]. 

The atomic wave-packet states {^'f(a;,i2nd + '^s(*))} of Eq. (104) show 
that just before the unitary accelerating sequence is switched on, all these 
possible wave-packet states of Eq. (104) are in the spatial region [xo{mr) — 
^difTT-r), xq{1) +£^(1)], whcrc XQ{j) and Sd{j) {j = l,2,...,m,.) arc the center- 
of-mass position and the wave-packet spreading of the j— th wave-packet state 
^j'{x,t2n^ + Ts{j)), respectively. Suppose that all these possible wave- 
packet states are accelerated uniformly by the unitary accelerating sequence 
and each possible wave-packet state moves the same distance La during the 
unitary accelerating process. The distance La will be obtained later. Obviously, 
after the unitary accelerating process all these to^ possible wave-packet motional 
states are in the spatial region [xo{mr) + LA — Sa{mr), xo{l) + LA + Sa{l)], where 
£a(j) {j = 1,2, ...,mr) is the wave-packet spreading of the j— th wave-packet 
motional state of the atom after the unitary accelerating process. Therefore, 
during the unitary accelerating process any possible wave-packet motional state 
of the halting-qubit atom is within the effective spatial region: 

[Al, Ar] = [xo{mr) - Ed, Xo{l) + La+ Sa] 
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where Ea » £a(l) and £d » £d{mr). The effective spatial region \Al, A^] 
covers all these possible wave-packet motional states of the atom during 
the whole unitary accelerating process. Now the spatial region of the unitary 
accelerating sequence must encompass sufficiently the whole effective spatial re- 
gion [Al, Afl], so that for all these possible wave- packet motional states the 
Raman laser light beams of the unitary accelerating sequence can be thought 
of as infinite plane-wave electromagnetic fields, and the most important is that 
the unitary accelerating sequence can act on all these possible wave-packet 
motional states simultaneously and uniformly during the whole unitary accel- 
erating process. The spatial region [Al, Afi\ may be called the accelerating 
spatial region. 

According to the quantum control process [22] the halting-qubit atom is ac- 
celerated by a unitary accelerating sequence at the end time of the quantum 
program. Here the unitary accelerating sequence may be given by [/A(2na) of 
Eq. (87a), which consists of tIq pairs of the basic STIRAP accelerating sequences 
(11a) and (63a) in an alternate form and each basic accelerating sequence has 
the same time period ta- The unitary accelerating process UAi'^ria) has a total 
time period 2nata- The ideal adiabatic condition (30) is also met in the uni- 
tary accelerating process. Now one may use the recursive relations (88) (94) 
to obtain the final wave-packet motional state of the halting-qubit atom after 
the atom is accelerated by the unitary accelerating sequence J7^(2nQ). Here 
the starting time of the imitary accelerating process is the end time of the 
quantum program. At the initial time each possible wave-packet motional 
state of the halting-qubit atom is given by "^f{x,t2nj + Ts{j)) of Eq. (104) 
for j = 1, 2, TO,-. All these to,- possible wave-packet motional states start 
to undergo the same unitary accelerating process C/A(2na) at the initial time 
trur simultaneously. In order to use the recursive relations (88)-(94) the initial 
motional state ^"0(2;,^^) of Eq. (88) needs first to be obtained from the state 
*f (a;, t2n^+Ts{j)) of Eq. (104). By comparing the initial state *g(a;, ig) of Eq. 
(88) with {x,tj„^ + Ts{j)) of Eq. (104) one can see that at the initial time 

= trrir- thc ccntcr-of-mass position, momentum, and global phase factor of the 
initial state M/g(a;,ig) are given by = z^ij) = + {PLJM)Ts{j), = 
Pt,, and exp[i(^g] = exp[i^g(j)] = eMivLj ^M-i{PLfT,{j) / {2nM)], re- 
spectively, and in the complex linewidth W{Ta) of the initial state 5'g(a;,tg) 
the time interval Ta = Ta{j) =Td + 2ndtd + Ts{j). It is known that the initial 
internal state is \go)- The initial atomic wave-packet product state then is given 
by ^'Q(a;,r, tg) = "^fi^Ttina + ^s(j))|<7o)- After the unitary accelerating process 
C/A(2na) the wave- packet motional state of the halting-qubit atom will take the 
form, according to the recursive relations (88)-(94), 




X exp{-- 



4 (Ax)2 + i 




(106) 
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and the atomic wave-packet product state is given by 

*2„.j(a;,r-,tLJ = *2„„j(^,*2nJl50), (107) 

where the end time of the unitary accelerating process is f^^ = tm,. + '^nata, 
and the atomic motional momentum is given by 

An„ = pL, + riaiHK + Hkl) + Uaini^ + Ml), (108) 

and the center-of-mass position (j) can be determined from the recursive 
relations: 

^2fe-l(j) = 4.-2(j) + -^t. - ^ / dt' COsHait'), 

^ Jtg + (2fe-2)t„ 

(109a) 

pa fcZ/a I la\ /■to+2kta 

iW ytg+(2fe-l)t<. 
PSk-l — P2k-2 + ^(^0 + ^l)' P2k = -^2^-1 + ^(^0 + ^l)' 

where 1 < fc < ria- The global phase factor exp[i(^2na0)] (^O^) ^1^° 

bo obtained from the recursive relations similar to Eq. (78) and (86). Now 
one can find from the final motional state *2n„j(^! *2ra„) of Eq. (106) that the 
moving distance La of the halting-qubit atom is La = -^fn^O) ~ ^o(i) during 
the imitary accelerating process, which appears in the accelerating region [Al, 
Ar] above. Note that the distance La is the same for each one of these rrir 
possible wave-packet motional states. 

The unitary accelerating process tells ones some facts. For the first point, 
the halting-qubit atom indeed is accelerated by na{fikQ+hkf)+na{hlQ+Hlf) and 
this accelerating process is uniform, that is, the accelerating process is the same 
for each one of these m,r possible wave-packet motional states {'^f{x,t2n^ + 
Ts{j))}. Thus, after the unitary accelerating process the atom is accelerated to 
the velocity {P2n^/M). For the second point, it can be seen from the motional 
states {^2na ^2na)} of Eq. (106) that in the complex linewidth the imaginary 
part increases linearly with the time period of the unitary accelerating process 
and is increased by h{2nata) / {'2M) , which is also independent of any index value 
J, while the real part keeps unchanged in the unitary accelerating process. For 
the third point, the distances between these nir possible wave-packet motional 
states keep unchanged during the unitary accelerating process. This fact can 
be deduced from the recursive relations (109a) and (109b) because the motional 
momentum P" {I = 0, 1,2, ...,2na), the mixing angles 9a{t) and 9ia{t), and the 
wave numbers (fcg + k^) and {Iq + l") all are independent of the index value 
j. This means that each one of these possible wave-packet motional states 
moves the same spatial distance during the unitary accelerating process. Since 
the distance AL{i,j) between the i— th and j—th {i < j) wave-packet motional 
states is still given by Eq. (100) and the atomic moving velocity is (P2„^/M) 
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after the unitary accelerating process, the time difference between the two wave- 
packet states *2n„.i(^7 *2n„) ^^'^ *2na j' (-^i '^2n„ ) ^hcn is givcn by 

AT{i,j) = AL{i,j)/{PS^JM) = (j - i)AT{PljPS^J, (110) 

where i < j; i, j = 1,2, ...,m,r. It is known that the time difference ATa{i, j) = 
{j — i)AT before the unitary accelerating process. Since the atomic veloc- 
ity {P2n^/M) after the accelerating process is much greater than the veloc- 
ity {Piud/-^'^) before the accelerating process, the time difference AT{i,j) << 
ATo{i,j), indicating that the time differences are compressed greatly for these 
rrir possible wave-packet motional states after the unitary accelerating process. 
Then the time-compressing factor for these possible wave-packet motional states 
after and before the unitary accelerating process i7^(2na) can be calculated by 

^ _ AT{i,j) 

* ATo{i,j) Pi^^ + na{nk^ + hki) + na{ni^Q + niiy 

The time-compressing factor Rt is independent of the indices i and j. Thus, 
the time-compressing process is uniform. The time-compressing factor Rt has 

been obtained in the previous paper [22], where Rt = (vq/v) and vq and v are 
denoted as the atomic moving velocities P^„^/M and P^„^/M before and after 
the unitary accelerating process, respectively. 

7. General adiabatic conditions and the error estimation for the 
decelerating and accelerating processes 

The starting point to set up a general adiabatic condition for a basic STI- 
RAP decelerating or accelerating process is to solve the basic equations (23) to 
find the coefficients {ak{P, t)} or to solve the basic equations (26) to obtain the 
coefficients {bk{P,t)}. Then it is to seek under what experimental conditions 
a real adiabatic condition for the basic STIRAP decelerating or accelerating 
process can be sufficiently close to the ideal adiabatic condition (30). This is a 
routine procedure in quantum mechanics [25] . There are three basic parameters 
to affect the real adiabatic condition of a STIRAP experiment: the time period 
of the STIRAP experiment, the Rabi frequencies and the phase-modulation 
functions of the Raman laser light beams. From the point of view of quan- 
tum computation one usually does not expect the quantum control process to 
consume a long time. However, a long time period of the STIRAP experiment 
usually can lead to that the adiabatic condition for the STIRAP experiment 
is met better [30]. If the time period of each basic STIRAP pulse sequence in 
the STIRAP-based unitary decelerating and accelerating processes is not long 
enough, then the adiabatic condition could not be met well. Then in this situ- 
ation one may use jointly the time period, the Rabi frequencies, and even the 
phase-modulation functions to achieve a better adiabatic condition for these 
decelerating and accelerating processes. Actually, the Rabi frequencies of the 
Raman laser light beams are very important to achieve a better adiabatic con- 
dition for the STIRAP experiment [15, 18b]. Without losing generality here 
take the basic STIRAP decelerating sequence (11) as an example to discuss a 
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general adiabatic condition. The obtained results ean be used as well for other 
basic STIRAP decelerating and accelerating processes. The STIRAP adiabatic 
conditions have been discussed in detail in many references [15, 16, 17, 18] in the 
conventional STIRAP experiments without considering explicitly the atomic or 
molecular momentum distribution. The conventional adiabatic conditions [4, 
15, 17, 18] usually are based on the first-order approximation solution to the 
basic equations similar to the present basic differential equations (26). These 
adiabatic conditions are usually a qualitative and approximate description to 
the adiabatic theorem. In the following two strict and different general adia- 
batic conditions are derived analytically. They are a quantitative description 
to the adiabatic theorem. The first general adiabatic condition is based on the 
Dyson series solution (29) of the basic differential equations (26). The second 
is based on a new method to solve the basic differential equations (26). This 
new method uses the equivalent transformations to solve the basic differential 
equations (26). That is, by making repeatedly the equivalent transformations 
the three basic differential equations (26) are transformed to the three equiva- 
lent linear algebra equations. Though the final solution to the basic differential 
equations (26) obtained with the new method is approximate, the truncation er- 
ror of the approximation solution can be controlled as desired. The two general 
adiabatic conditions may be used to set up the conventional STIRAP experi- 
ments. Thus, they may be used to design the STIRAP pulse sequence to realize 
the perfect state (or population) transfer for a quantum ensemble of the atoms 
or molecules. But their more important application is that they may be used 
to set up the basic STIRAP unitary decelerating and accelerating processes for 
a free atom and an atomic or molecular ensemble. 

The basic differential equations (26) or their matrix form (28) can be in- 
tegrated formally. The formal solution to the basic equations (28) may be 
expressed as the Dyson series (29). Here one needs to use the initial condi- 
tion of the basic STIRAP decelerating sequence (11). At the initial time to of 
the basic STIRAP decelerating sequence (11) the three-state vector B{P,to) = 
{ba{P,to),b+{P,to),b_{P,ta))'^ is given by Eq. (39). The initial condition (39) 
has been used to set up the ideal adiabatic condition (30). For a real adiabatic 
condition the initial condition may be generally given in (130) below. At first 
the formal solution (29) may be rewritten as 



where to < t < to+T and T is the time period of the basic STIRAP decelerating 

process, and the error term Er{P,t) measures the deviation of a real adiabatic 
condition from the ideal adiabatic condition and it may be expressed as 



B{P,t) = B{P,to) + Er{P,t) 



(112) 




+ i-f I r r dtidt2dt3M{P,ti)M{P,t2)M{P,t3) + ...}B{P,to). (113) 
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The upper bound of the error term is evaluated accurately below. Denote the 
maximum norm of the hermitian matrix M(P, tn) which is given in (28) in the 
time region [to,tn-i] as 

||M(P,f„)l|max = ^ max {\\MiP,tn)\\}, 

where to <•••<*« < ^n-i < ••• < ti <to + T. Obviously, there are the following 
relations for the maximum norms {||M(P,t„)||max} '■ 

||M(P,io)|Ua.< •••< ||M(P, 

) 1 1 max < \\M{P, 

^n— l) 1 1 max 

< ||M(P,t2)|Uax < \\M{P,h)\U^ = \\M{P,t)\U^. (114) 
Here the maximum norm ||M(P, f)||max is defined as 

\\MiP,t)\U. = ^ max {||M(P,i)||}. (115) 

to<t<.to + -l 

Then with the help of (113) and (114) it can turn out that the upper bound of 
the deviation Er (P, t) may be determined from 

\\EriP,m < exp[(||M(P,t)|Uax)r] X ||P(l)(P,t)|kax, (116) 

where the first-order approximation solution B^^\P,t) to the basic differential 
equations (26) is given by 

1 ft 

BW(P,i) = (-) / dhM{P,h)B{P,to), {to<t<to + T), (117) 
« Jto 

while ||P(^^(P, f)||max is the maximum norm of the solution B^^\P,t) in the 
time region to <t <to+T. This norm ||i?(^)(P, t)||i„ax is written as 

||B('H^,i)l|max = V(l''o'H^,i)P + l&+^(^.*)P + l&-^(^.*)PW. (118) 

On the other hand, it follows from the matrix M(P, t) in (28) that the maximum 
norm | |M(P, f)| |max is bounded by 

||M(P,i)|U,.<^X]|My(P,t)|2< ^{2\Q{P,t)\ + \T{P,t)\U^. (119) 

The adiabatic condition (116) is strict because it is required that at any instant 
of time in the whole STIRAP decelerating or accelerating process the deviation 
from the ideal adiabatic condition (30) be limited within a given small value, 
that is, the upper bound of the error term (P, t) is less than a given small 
value at any instant of time. Notice that in theory at the initial time <o the 
atom is prepared to be in the trapped state \g^{P,to)) of (19a) completely. If 
the error term Er{P,t) is large, then this will mean that during the STIRAP 
decelerating or accelerating process there is a large probability for the atom to 
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be excited to the two cigcnstatcs \g^{P,t)) of the instantaneous Hamiltonian 
H{P,t) of (17). It is known from (19b) that any one of the two eigenstates 
\g'^{P, t)) contains the excited internal state of the atom. Then the atom could 
be easily affected due to the atomic spontaneous emission if it is in any one 
of the eigenstates \g'^{P,t)). On the other hand, a high probability for the 
atom to stay in the trapped state |5'''(P, t)) may lead to that the atom is not 
easily affected by environment and may avoid the spontaneous emission. The 
adiabatic condition (116) indicates that the probability for the atom to leave 
the trapped state \g°{P,t)) may be limited to a small value as desired during 
the STIRAP decelerating or accelerating process. Therefore, it ensures that the 
atom is almost completely in the trapped state \g°{P,t)) during the STIRAP 
decelerating or accelerating process. The adiabatic condition (116) is more 
severe than those in the conventional STIRAP experiments [15, 17, 18]. The 
latter usually require that the probability for the atoms or molecules under 
investigation in the two eigenstates \g^{P,t)) be much smaller than one. This 
is a qualitative description for the adiabatic theorem. The present adiabatic 
condition (116) is closely related to the requirement that Gaussian shape of the 
Gaussian wave-packet motional state of the decelerated or accelerated atom keep 
unchanged before and after the basic STIRAP decelerating and accelerating 
processes. It measures the deviation of a real adiabatic condition from the 
ideal adiabatic condition (30), while the deviation may occur not only in the 
two cigcnstatcs \g^{P,t)) but also in the trapped state \g^{P,t)). The present 
adiabatic condition (116) limits the upper bound of the deviation to a given 
small value. This is a quantitative description for the adiabatic theorem. This 
results in that the present adiabatic condition (116) is more severe. 

When the adiabatic condition (116) is met, the error term \ \Er{P,to + T)\\ 
of the final state at the time t = to +T is clearly not more than the upper 
bound (116) and the real error term \\Er(P,to + r)|| could be much less than 
the upper bound (116). It may be required in theory that the real error term 
\\Er{P,to + T)\\ of the final state be less than some given value which is much 
less than the upper boimd (116). This requirement is not severe in theory with 
respect to the adiabatic condition (116). It may be met by setting the suitable 
experimental parameters at the final time t = to+T for the STIRAP decelerating 
or accelerating process. However, in practice the lower bound of the error term 
Er{P,tQ + T) of the final state is generally affected by the adiabatic condition 
(116). If the upper bound (116) is large, then the lower bound of the error term 
Er{P, to + T) usually is large too. 

According to the superposition principle in quantum mechanics in a real 
adiabatic condition the atomic product state at any instant of time t {to < t < 
to + T) in the basic STIRAP decelerating process (11) may be calculated from 
Eq. (12) {P = P'- hko), 

\^r{x,r,t)) =Y,p{P){[A{P,t)+5o{PM\P + r^ko)\go) 
p 

+ [A[{P,t)+5i{P,t)\\P)\e) + [A^^{P,t)+5^{P,t)]\P - (120) 
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where the coefficients {Al,{P, i)} for fc = 0, 1, 2 are obtained in the ideal adiabatic 
condition (30), while the coefficients {6^{P,t)} measure the deviation of the 
real adiabatic condition from the ideal one. The product state (120) may be 
rewritten as 

\^'r{x,r,t)) = \^,{x,r,t)) + Er{x,r,t) 

where the wave-packet state |\E'i(x,r, t)) is the atomic state at the time t in the 
basic STIRAP decelerating process (11) in the ideal adiabatic condition and it 
may be written as 

\^i{x, r, t)) = J2 PiP){Ai{P, t)\P + Hko) \9o) 
p 

+ Ai(P,t)|P)|e) +A^(P,i)|P-fifci)|5i)}, (121) 
and the error term Er{x, r, t) is given by 

Er{x, r,t) = J2 PiPK^oiP, t)\P + fiko)\go) 
p 

+ 6t{P,t)\P)\e)+di{Rt)\P-Hk,)\g,)}. (122) 

It turns out in the preceding section 5 that the final state |\E',(a;,r, t)) with 
t = to+T in the ideal adiabatic condition is a perfect Gaussian wave-packet state 
if the initial state of the basic STIRAP decelerating process (11) is a Gaussian 
wave-packet state. Obviously, it follows from (122) that the probability for the 
error term Er{x, r, t) at any time t may be calculated by 

\\E^{x,r,t)\\'' = E |p(^')r{l^o + \StiP,t)\' + \6UP,m- (123) 

p 

(Notice that the error probability (121) in the previous versions of this paper 

which is denoted as Er{P,t) is just equal to 2\\Er{x,r,t)\\'^ of (123)). In or- 
der to use directly the solution to the basic differential equations (26) or their 
matrix form (28) to calculate the upper bound of the error term Er{x, r, t) one 
may use the coefficients bf){P,t) and b±{P,t) to express the error probability 
||£'r-(a^)?')i)|P of (123). Notice that there is the unitary transformation between 
the two three-state vectors {bo{P,t),b+{P,t),b-{P,t))'^ and {Ao{P,t), Ai{P,t), 
A2{P,t)Y . The three-state vector _(Ao(P,t)^Ai(P,t)^A2(P,t))^ is ffi-st con- 
verted into the three-state vector {AQ{P,t), Ai[P,t), A2{P,t))^ by the unitary 
transformation of (15a)-(15c), then into the three-state vector (ao(P, t), a+{P, t), 
a-{P,t))'^ by the unitary transformation of (22a)-(22c), and finally into the 
three-state vector {bo{P,t), b+{P,t),b-{P,t))'^ by the unitary transformation 
(25). Thus, under these unitary transformations there is the relation: 

iAoiP,t),A,{P,t),A2{Pt)f = UAb{boiAt),b+{P,t),b_{P,t)f (124) 

where UAb is the unitary transformation between the two three-state vectors 
{bo{P,t),b+{P,t),b-{P,t))'^ and {AQ{P,t), Ai{P,t), A2{P,t))'^ . If now the so- 
lution to the basic equations (26) in the ideal adiabatic condition is given 
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by (bQ{P,t),b^_^{P,t),b^_{P,t))'^ , then after the unitary transformation UAb one 
obtains the three-state vector {AQ{P,t), A\{P,t), A2{P,t))'^ of the ideal adi- 
abatic condition and then the state \^i{x,r,t)) can be calculated from Eq. 
(121) by using the thrcc-statc vector. If the solution to the basic equations 
(26) in a real adiabatic condition is given by {bo{P,t),b+{P,t),b-{P,t))^ with 
bk{P, t) = b\{P, t) + d\{P, t) for k = 0,+,—, then after the unitary transforma- 
tion UAb one obtains the thrcc-statc vector {Ao{P. t), Ai{P, t), A2{P, t))'^ of the 
real adiabatic condition, where Ak{P, t) = Al.{P, t) + S^{P, t). Thus, there is the 
unitary transformation between the two three-state deviation vectors: 

{S^{P,t),St{P,t),S^{P,t)f = UAb{5l{P.t),6\{P,t),5t{P,t)f. (125) 

It is well known that the unitary transformation UAb does not change the norm 
of the three-state deviation vector {5^Q{P,t),5''j^{P,t),5^_{P,t)Y' . This indicates 
that there is the relation: 

|5^(P,i)|2 + |<5f (P,i)|2 + |<5^(P,i)|^ = \6l{P,t)\'' + \5\{P,t)\'' + \6'_{P,t)\\ 

This relation leads to that the error probability ||£'r-(a;,r, f)|p of (123) may be 
expressed as 

|K(a;,r,t)||2 = ^ \p{Pt{\5l{P,tt + I^U^.*)!' + (126) 
p 

It is convenient to calculate the error upper bound [[-©^(x, r, <)|| by using the 
equation (126), since the deviation vector {5l{P,t),5\{P,t),5_{PJ)Y 

can be 

obtained conveniently by solving the basic differential equations (26). Thus, an 
accurate error upper bound | |£'r(a;, r, t)| | could be obtained directly by comput- 
ing the equation (126) by using the deviation vector {5^q{P, t), 6''_^_{P, t),5^_{P, t))'^ 
for the basic STIRAP decelerating or accelerating process. Obviously, the three- 
state deviation vector ((5q(P, f), (5']_(P, t), 5^ (P, f))^ has the maximum norm or 
length over the effective momentum distribution region [P] and in the time 
period Iq <t <to+T, 

^|,5^(P,t)|2+|J^(P,i)|2 + |5i(P,i)|2 

< \] {\5l{P,t)\^ + \5X{P,t)\^ + |(5^_(P,t)|2)„,„ for P G [P] and to<t<h+T. 
Then the upper bound of the error term Er{x, r, t) may be determined from 

\\Er{x,r,t)\\ < yJ{\5l{PM^ + + l-^-C^. i)l')xnax, (127) 

where the normalization relation X^p|p(P)P = 1 is used and the truncation 
error is neglected for any momentum components outside the effective momen- 
tum region [P]. The inequality (127) is a general adiabatic condition for the 
basic STIRAP decelerating or accelerating process of a free atom in a wave- 
packet motional state. There is also a simpler method to obtain the error 
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upper bound ||£'r(a;, r, i)||, as stated below. It uses the general adiabatic condi- 
tion (116). It is known that the formal solution to the basic differential equa- 
tions (26) or their matrix form (28) may be expressed as (112), where the 
solution in the ideal adiabatic condition (30) is given by B{P,t) = _B(P, to), 
as shown in Eq. (40) in the previous section 4. Then the equation (112) 
shows that the three-state deviation vector is just Er{P,t) and hence one has 
Er{P,t) = (SgiPJ,), S''_^{P,t), 6t{P,t))'^ . Furthermore, the adiabatic condition 
(116) and the equation (126) show that there are the relations: 

\\Erix,r,t)\\ = {T.\p{P)mr{P,t)\\''r^' 

p 

< {Y,\p{P)\'eMmMiP,t)\U.)T] X \\B^'\P,mle.V^' 
p 

<exp[(||M(P,i)|Uax)T] X ||p(i)(P,t)|Uax, (128) 

where the relation J^p |p(-f )P = 1 is used and the truncation error has been ne- 
glected for any momentum components outside the effective momentum region 
[P] , and the maximum norms | |M(P, i) | |max and | \B^^^ (P, t) \ l^^x are respectively 
defined as 

(||M(P,t)|Ux) = max^(||M(P,t)||^ax), 

1 1^(1) (P, t) I Uax = max I |P(1) (P, t)\\n,^. 

The last inequality in (128) is a real adiabatic condition of the basic STIRAP 
decelerating or accelerating process for a free atom in a wave-packet motional 
state. It could be useful to design the basic STIRAP decelerating or accelerating 
process. 

At first the adiabatic condition (128) requires one to calculate the norm 
(||M(P, t)||inax) and the first-order approximation solution B^^\P,t). It is easy 
to calculate the first-order approximation solution to the basic differential equa- 
tions (26). Actually, the first-order approximation solution may be obtained 
from the equation (117). Here for convenience setting the global phases 7(to) = 
(5(to) = in the basic equations (26) and the phase-modulation fimctions of the 
Raman laser light beams to be </'o(0 = in Eq. (31) and ipi{t) = in Eq. (32). 
It should be pointed out that the following methods are available as well for the 
phase-modulation Raman laser light beams. It follows from Eqs. (31) and (32) 
that 

d , AP, d , AP, 

— a„{P,t) = kt), — a^(P,t) = fci . 

dt ' ' M dt ^ ' ' M 

By inserting these two equations into Eqs. (27a)-(27c) one obtains 



n±(P, t) = n{t) ± Ko{t)AP, (129a) 

e(p, t) = -e{t) + iKi{t)AP, r(p, t) = ^ap - K2{t)AP, (129b) 



62 



where 

Ko{t) = — [(fco-fci) + 3(fco + fci)cos20(t)], 

Ki{t) = sm2e{t), K2{t) = — cos 9{t). 

The initial condition for the basic equations (26) is given by (39). If the initial 
mixing angle ^?(io) is very small {9 (to) « 1) but not equal to zero, then the 
initial condition (39) may be changed to the general form 

6o(P,io) = exp[- C ' +Eo)to]cos0{to), (130a) 

b+{P,to)=b.{P,to) = -^exp[-( ^ +Eo)to]sme{to). (130b) 

Then in the initial condition (130) the first-order approximation solution to the 
basic equations (26) for the coefficient bo{P, t) may be written as, by integrating 
by parts the integral (117), 

b^^^ (P, t) = bo{P, t) - bo{P, to) = C^^^ {P, t) + C^^^ (P, t), (131a) 

where the main term C'q^\p, t) that is proportional to 6(P, t)* /fl{t) is written 



as 



C^^\P,t) = V2b+{P,to) ^^^;!'^* exp[zAP f dt'Ko{t')]sin[i f dt'^{t% 

^'■V-) J to J to 

(131b) 



(2) 

and the secondary term Cqj! (P, t) is given by 



rt 



C^,^{P,t) = ^^/26+(P,^o) dt,{[z—i^^) AP] 



xexp[iAP dt'Ko{t')]sm[i dt'fl{t')]}. (131c) 

Jto Jta 

The first-order solution for the coefficients b± (P, t) is given by 

bi\P,t) = b±{P,t)-b±{P,to) = c2\P,t) + C^^l{P,t) + Fi'\p,t)+Fi'UP,t) 

(132a) 

where the main terms C^^ (P, t) are given by 

C^±\P,t)=Ti^bo{P,to)^^exp[-iAP dt'Ko{t')]exp[Ti f dt'n{t')] 
v2 Jto Jto 
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and the secondary terms are 

c^^{p,t) = ±-^bo{p,to) "^^'^^'df.^^^^^y^ + — fx^o — ^ 

xexp[-iAP f dt'Ko{t')]exp[Ti f dt'n{t')]}, (132c) 

J t(j Jto 

(132d) 

F^^{P,t) = T^-b^{P,to) dh{[-^^i^^)]eMTi dt'2Q{t')]}. 

(132e) 

The dominating terms in the first-order approximation solution of (131a) and 
(132a) arc c'^^ (P, t), which are proportional to 6(P, t)/il{t) and bo^P, to). It can 

(2) 

turn out by integrating by parts the integral (132c) that the terms C^j,{P,t) 

are proportional to f](f)^^. Thus, the terms Cj3^(P, t) are secondary in the first- 
order solution for a large Rabi frequency Q{t). Similarly, it can turn out that 
C^y^(P, t) (X b+{P, tQ)n{t)-'^ and F^^{P, t) ot b+{P, ^o)^^(^)"^ by integrating by 
parts the integrals (131c) and (132e), respectively. Thus, these terms CQr^{P,t) 

and (P, t) are secondary with respect to the terms c'^\P,t) and F^^^ {P,t)i 
respectively. On the other hand, the imperfection for the initial conditions 
b+{P,ta) = b^{P,to) ^ could mainly affect ci'^\p,t) and F^\p,t). Its mag- 
nitude is approximately proportional to the factors \b±{P,to)\\Q{P,t)\/rt{t) or 
\b±{P,to)\\T{P,t)\/n{t). Since \b±{P,to)\ « \bo{P,to)\, these terms C^^\P,t) 
and F^^ (P, t) are secondary with respect to the main terms C^"^ (P, t). Thus, the 
error upper bound ||i?'^^^(P, t)||max in (128) may be determined from the main 
terms C^\P,t). Now by inserting 6^^^(P,t) of (131a) and b'i\P,t) of (132a) 
into (118) it can be found that the norm ||B(i)(P,t)|| is bounded by 

iiBWrp/^ii<M^ + M^ 



^Je{t)^ + K,{tr\AP\^ ^Je{tor + K,ito)^\AP\^ 



where those secondary terms of the first-order approximation solution are ne- 
glected and only the main terms c£^'' (P, t) of (132b) arc used and |&o(^: ^o)l ^ 1 
is also used. It is known that the initial mixing angle 9{to) 0, as can be seen 
in (35). Note that Ki{t) = (fco + h) sin 20{t)/{2M) and the time derivative l9(i) 
of the mixing angle is given by 

.... »p(/)co.s^(/)-r),(/)smg(t) 

^^^^ = m ■ 
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If the Rabi frequencies Qp{t) and fls{t) arc chosen suitably in experiment such 
that at the initial and final times the mixing angle satisfies the relations: 

tan6'(to) = np{to)/ns{to) « np{to)/n,{to) 0, 

f2p(io + T) cos 9{to + T) - n,{to + T) sin e{to + T) ^ 0, 

then the time derivative ^(to) ~ and 6{ti^ + T) « 0. On the other hand, the 
momentum distribution satisfies |AP| < /S.Pm/'2- for a momentum wave-packet 
state with an effective momentum bandwidth APm- Therefore, at the initial 
time |ii'i(to)AP| < |ii'i(to)|APA//2 = APM(fco + fci)| sin26'(to)|/(4M) « if 
the momentum wave-packet state has a finite wave-packet spread. Then in these 
conditions the error upper bound | (P, t)| l^ax may be determined from, by 
neglecting the second term on the rightest side of (133), 



iiB<->(p,oii.( ^^"''"'^^"'''\;;'''-'°'^^"'^'^^"^ .„. (134) 

Here the subscript 'max' means that the function on the right-hand side of 
(134) is taken as the maximum value in the time period < t < to -\- T 
of the STIRAP process. The first-order approximation adiabatic condition is 
that the maximum value of the function on the right-hand side of (134) is 
controlled to be smaller than some desired small value. (Notice that this first- 
order adiabatic condition (134) is slightly different from that one (128a) in the 
previous versions of this paper). As shown below, in the initial and final time 
periods of the STIRAP process the adiabatic condition (134) still may be met 
even if the Rabi frequency ^{t) is small in these time periods. This is a global 
adiabatic condition, since it is involved in the whole time period of the STIRAP 
process. Here the global adiabatic condition has a different definition from the 
conventional one in Ref. [4]. In the conventional STIRAP experiments [4, 15, 
17, 18] the (local) adiabatic condition is defined as that at any instant of time 
of the STIRAP process the population or probability in the two eigenstates 
|(7^(P, t)) that contain the excited internal state is much smaller unity. This 
is approximately equivalent to the first inequality of (133) for the first-order 
approximation, where the inequality symbol ' <' is replaced with ' «' . On 
the other hand, according to (119) the maximum matrix norm (||M(P, <;)||inax) 
is determined from 

||M(p,t)||<i={2|e(p,i)| + |r(p,t)|} 



< 



^{2Vl^Wlmax + {ko + fci)2(APM)V(16M2) + max(fco, h)} 

(135) 

where |^(t)|max is the maximum value of \9{t)\ in the whole STIRAP decelerating 
or accelerating process. After the upper bound \\B^^\P,t)\\m_a.^ and maximum 
norm (||M(P, f)||max) are determined from (134) and (135), respectively, one 
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may determine the upper bound of the error term Er{x, r, t) from (128). It can 
be seen from (134), (135), and (128) that the adiabatic condition (128) may 
be better satisfied for a small time derivative 6{t), a short time period T, a 
large Rabi frequency rJ(t), and a narrow momentum wave-packet state (APm is 
small). If the time derivative 9{t) is smaller, then the time period T usually is 
larger. Thus, there is a compromise between the settings of the time derivative 
6{i?i and the time period T in experiment. Obviously, one has the relation for 
any basic STIRAP decelerating or accelerating process: 

6»(io + T) - e{to) = f"^ dt'9{t') = 7r/2. 

J to 

This relation may be used to determine the time period T if one knows the time 
derivative 9{t) of the mixing angle. 

It seems that the adiabatic conditions (116) and (128) could not be better 
satisfied in the initial and final time periods, since the Rabi frequency Q{t) 
takes a smaller value in these time periods. It is well known that any STIRAP 
experiment requires that at the initial and final time the mixing angle 9{t) satisfy 
the constraint conditions: 

9{to) ^ 0, 9{to +T)^ Tr/2. 

The two constraint conditions are compatible with the adiabatic conditions (116) 
and (128). This can be seen from (134). At the initial time period the Rabi 
frequency fl{t) takes generally a smaller value, but at the same time the mixing 
angle 6{to) and its time derivative ^(to) may be set to a value close to 
zero, leading to that the value | (P, f)| | may be kept at a smaller value at 
the initial time period. At the final time the mixing angle 9{to +T) 7r/2 or 
sin^ 29(to +T) and the time derivative ^(^o + T) also may be set to a value 
close to zero. Then the value ||-B*^^)(P, t)|| still may be kept at a smaller value, 
although the Rabi frequency Cl{t) takes a smaller value at the final time period. 
These results show that in theory the adiabatic conditions (116) and (128) still 
may be met in the initial and final time periods as long as the Raman laser 
light beams of the STIRAP experiment are suitably designed. The adiabatic 
conditions (116) and (128) could be better used for a conventional three-state 
STIRAP experiment that uses a pair of copropagating Raman laser light beams 
and those STIRAP-based decelerating and accelerating processes of the atomic 
or molecular systems with a narrow momentum distribution. 

The deviation Er{P,t) of (113) and its upper bound (116) are obtained for 
a single basic STIRAP decelerating sequence (11). If a unitary decelerating 
process consists of rid pairs of the basic STIRAP decelerating sequences (11) 
and (63), then the total deviation generated in the unitary decelerating process 
is bounded by 

\\Er{P,t)\ \ <ndexp[(||Mfe(P,i)||n,ax)T] X ||Pi''(P,i)|Uax 

+ ndexp[(||M;(P,i)|Ua.)r] X ||B|')(P,t)||^ax, (136) 



66 



where the subscript k marks the basic STIRAP decelerating sequence (11) that 
uses a pair of the Raman laser light beams with the Rabi frequency ^{t), the 
mixing angle 6{t), the wave numbers fci and k2, and the carrier frequencies uiqi 
and ujQ2, while the subscript / denotes the basic STIRAP decelerating sequence 
(63) that may use another pair of the Raman laser light beams with the Rabi 
frequency fli{t), the mixing angle Oi{t), the wave numbers k[ and ^2, and the 
carrier frequencies uj^i and cjq2- Both the basic STIRAP decelerating processes 
have the same time period T. The upper bound of the total deviation Er{P, t) 
on the right-hand side of the inequality (136) for the unitary decelerating or 
accelerating process may be controlled by setting suitably the experimental 
parameters of the Raman laser light beams, which include the Rabi frequencies 
fl{t) and fli{t), the mixing angles 6{t) and Oi{t), and the time derivatives (){t) 



The basic STIRAP decelerating or accelerating process is a time-dependent 
unitary quantum dynamical problem from the viewpoint of quantum mechan- 
ics. The three-state basic STIRAP decelerating or accelerating process is a quite 
simple unitary dynamical process, but it can describe completely the complex 
STIRAP-based unitary decelerating and accelerating processes of a free atom. 
It is relatively simple to solve approximately such a unitary dynamical problem 
as the three-state basic STIRAP decelerating or accelerating process in quan- 
tum mechanics, although this problem is time-dependent and it is difficult to 
solve exactly the basic differential equations (26) except for some special cases 
[39]. For example, one may use the conventional methods of successive approx- 
imations [33] to solve these basic differential equations approximately when the 
Rabi frequencies il{t) is large. As mentioned before, the Dyson series solution 
(29) to the basic differential equations (26) or (28) may be used to calculate 
the deviation of a real adiabatic condition from the ideal adiabatic condition. 
The problem to be answered is that one needs to calculate how many leading 
terms in the Dyson series (29) so that the result obtained is enough accurate. 
The leading term number is dependent upon the desired error value and the 
maximum norm of the matrix M{P,t) in (28). If one uses the first n terms 
on the right-hand side of (113) to calculate the error term Er{P,t), then the 
residual term may be given by 



and Oi {t), and so on. 




dtidt2...dtn+lM{P, ti)M{P, t2)...M{P, tn+l) 




dtidt2...dtn+2M{P,ti)M{P,t2)...M{P,tn+2) 



+ }B{P,to). 

Then the residual term R{P, t) is bounded by 



(137) 



CO 



k\ 



\\R{PM< E 



m{p,t)\u^f 



k=n+l 
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< '"^'^■""Ti',-'°''"" exp|(||«(P,,)||„„)(,-MI (138) 



where \\B{P,to)\\ = 1 is used. Now (n+ 1)! w v^7r(n+l)[(n + l)/e]"+i. Then 
the upper bound of the residual term is determined from 

s/2TT(n + 1) 

If (||M(P,i)|U^)(i - to) <{n + 1), then expi^^Mpf^H^^] < e. Thus, 

if (n + l)/e > ||M(P,t)|Uax(i - to) exp[^Mmkp)(£^], then the residual 
term t) is exponentially small. Suppose that the upper bound of the error 
term Er{P^ t) is set to a given value Sr ■ \ \Er{P,t)\ \ < Sr- Then this requires the 
residual term to satisfy << Sr- The condition ||i?(P,f)|| << Sr can 

be easily met by setting a minimum integer n such that 

1 l|M(P,t)|| ..T (||M(P,t)|U..)T 

V27r(n + 1) ^ (n + l)/e n + 1 r v ; 

Once the minimum integer n is determined, one may use the first n terms on 
the right-hand side of (113) to calculate the error term Er{P,t) and its upper 
bound, while the residual term R{P,t) does not affect significantly the final 
result. 

It can be seen from (133) that the first-order approximation solution shows 
that the error upper bound ||B'^^(P, t)|| is mainly dependent upon the parame- 
ter |0(P, f)|/ri(t) and almost independent of the parameters r{P,t) and KQ{t). 
Actually, in the first-order approximation solution the parameter T(P,t) ap- 
pears in the secondary terms F2\p,t) and F^^{P,t) and i^o(ii) in the sec- 
ondary terms c!^^{P,t) (these parameters may appear in the phase factors, 
and if so, they do not make a contribution to the error upper bound). It is 
known that the parameter 0{P,t) = -9{t) + i^{AP/M){ko + h) sm2e{t). If 
the two Raman laser light beams are copropagating, then the wave-number sum 
(fco + ki) will be changed to the wave- number difference (/cq ^ ^i) in the param- 
eter 0{P,t). Then the effect of the momentum distribution on the error upper 
bound ||P(^)(P, t)|| will be greatly weakened and the momentum distribution 
will become a higher-order effect on the STIRAP state transfer. Therefore, in 
this sense the copropagating Raman laser light beams used to construct the 
STIRAP pulse sequence may be better than the counterpropagating ones to re- 
alize the perfect STIRAP state transfer. Unlike the parameter &{P, t) these two 
parameters r(P, t) and Koit) are dependent on both the wave-number sum and 
difference. That the momentum distribution affects the STIRAP state trans- 
fer is mainly through the parameters r(P, t) and Ko(t), which appear in the 
higher-order terms in the Dyson series (29), if the copropagating Raman laser 
light beams are used in the STIRAP state transfer. In fact, the maximum norm 
(||M(P, t)||inax) determined from (135) shows that it is proportional to |r(P,t)|. 
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These parameters make an important effect for tire momentum distribution on 
the STIRAP state transfer. Thus, that the conventional three-state STIRAP 
experiments [4, 15] use the copropagating Raman laser light beams is favorable 
for the perfect state transfer and may minimize the effect of the momentum 
distribution of the atomic and molecular systems under investigation on the 
perfect state transfer. However, in the basic STIRAP decelerating or accel- 
erating process the countcrpropagating Raman laser light beams are generally 
used so that the fast moving atom can be decelerated or accelerated more ef- 
ficiently. On the other hand, the first-order solution (117) could not exactly 
account for the momentum distribution in any case that either copropagating 
or countcrpropagating laser light beams is used in the STIRAP experiments. 
The adiabatic conditions (116) and (128) are accurate, but due to that there is 
an exponential correction factor in (116) and (128) they could not be met for a 
broad momentum distribution. A broad momentum distribution is often met in 
an atomic or molecular quantum ensemble. It is necessary to consider the effect 
of the momentum distribution when these physical ensembles are decelerated 
(or accelerated) by the STIRAP decelerating (or accelerating) pulse sequence. 
Thus, it is necessary to find a more useful adiabatic condition that can account 
for the effect of a broad momentum distribution on the STIRAP state transfer. 

It is still complex to use the Dyson series solution (29) to calculate the error 
term Er{P,t) of (113) and its upper bound. In the following an equivalent 
transformation method based on the integration by parts to solve the basic 
differential equations (26) is proposed so that the error term Er{P,t) of (113) 
and its upper bound can be obtained conveniently in a high accuracy. On the 
other hand, by solving the basic differential equations (26) to obtain an enough 
accurate solution one may further use the solution to calculate conveniently 
the time evolution process for the basic STIRAP decelerating and accelerating 
processes for a free atom. Generally, it is quite inconvenient to calculate the time 
evolution process of an atomic decelerating or accelerating process by directly 
solving the Schrodinger equation. The present scheme is convenient to calculate 
the time evolution process of the basic STIRAP decelerating or accelerating 
process because it does not solve directly the original Schrodinger equation 
but solves the three first-order differential equations (26) that are equivalent 
to and much simpler than the original Schrodinger equation. The equivalent 
transformation method to solve the basic equations (26) is based on the fact that 
the Rabi frequency Cl{t) may be set to a large value in experiment. Though the 
solution to the basic equations (26) obtained by this method is approximate, the 
truncation error of the solution can be controlled as expected. The procedure 
to solve the basic equations (26) with the equivalent transformation method 
may be described below. By integrating the basic differential equations (26) 
one obtains the equivalent integral equations: 

1 /■* 

bo{p,t) - bo{p,to) = ^j^ dti{b+{p,h)e{p,tiy 
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xexp[iAP / dt'Ko{t')]exp[i j dt'Q.{t')]} 

J to J to 

f dti{b-{P,ti)e{P,ti)* exp[iAP f ' dt'Ko{t')]exp[-i [ ' dt'n{t')]}, 

V 2 J to J to J to 



to '^to 

(141a) 

t rti 



b±{P,t) - b±{P,to) = -i\ I dti{6T(P,ii)r(P,ii)exp[Ti / ' dt'2^{t')]} 

^ Jta J to 



^ / dti{bo{P,ti)@{P,ti)exp[-iAP [ ' dt'KQ{t')]exp[Ti j' dt'n{t' 

V2 Jto Jto Jto 



)]}■ 



(141b) 

Hereafter the initial condition (39) is used for convenience. The equivalent trans- 
formation method is that by the integration by parts and another transformation 
(see below) the integral equations (141) may be approximately reduced to the 
three linear algebra equations. These three linear algebra equations are equiva- 
lent to the original integral equations (141) if the initial condition (39) is taken 
into account and when the truncation error can be neglected. At the first step 
of the equivalent transformation method the integrals on the right-hand sides 
of (141) are calculated by the integration by parts. Then the initial condition 
(39) is used to simplify the calculated results If there are the time derivatives 
of the variables 6o {P, t) and b± {P, t) in the integrands after the integration by 
parts, then one may use the basic differential equations (26) to replace these 
time derivatives. As an example, by integrating by parts the equation (141b) 
for the variable 6+ (P, t) and then using the basic differential equations (26) and 
the initial condition (39), one can obtain the following equation: 

b\{P,t) = b+{P,t) = b-{P,t)T-^{P,t)ex-p[-i j dt'2VL{t')] 

Jto 

+ibo{P,t)T\^{P,t)exp[-iAP [ dt' Ko{t')]exp[-i [ dt'n{t')] 

Jto Jto 

-ibo{P,to)Tl,{P,to)+i f dh{b+{P,h)Tt{P,h)} 

Jto 

+i I dii{6_(P,ii)e:^i(P,ti)exp[-i dt'2n{t')]} 

Jto Jt() 

+ I dti{6o(J^,ii)6+i(-P,ii)exp[-iAP / dt' Ko{t')\exp[-i j dt'n{t')]} 

J to J to J to 

(142) 

where the five amplitudes are given by 

r-(Pt)-ilW) (Pi)__J_^ 
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e- (P t) - + MM^i 



ftO .p,^_ 1 ^■d MP,t) i mt)Q{P,t) Ko{t)e{P,t) 

The equation (142) is almost completely equivalent to the original equation 

(141b) for the variable b^{P,t). The unique difference between the two equa- 
tions is that the equation (142) uses the initial condition (39), while the original 
equation does not. The first four terms on the right-hand side of (142) may 
be considered as the main terms, since these terms have a greater contribu- 
tion to the solution b}^{P,t) of (142). On the other hand, each one of the 
last two integrals on the right-hand side of (142) contains the integrands with 
the largely oscillatory phase factor exp[— i J^*^ dt'V,{t')] or exp[— i J^*^ dt'2V,{t')]. 
These largely oscillatory phase factors make the two integrals secondary in the 
solution (142). This can be seen by integrating by parts the two integrals once 
again. If these two integrals are neglected, then one obtains the first-order ap- 
proximation solution to the original equation (141b) for the variable b+{P,t), 
since all the five amplitudes including r^j(P, f), r°^(P, t), etc., appearing in 
the equation (142) are inversely proportional to the Rabi frequency n{t). 

One may further obtain a better approximation solution than the first-order 
one. This can be done by integrating by parts the last two integrals on the 
right-hand side of (142) again. However, the fourth term (or the first integral) 
contains the solution b+{P,t) itself on the right-hand side of (142). While one 
may substitute the first-order approximation solution of b\ {P, t) into the integral 
to obtain a better approximation, the calculation process becomes so complex 
that one can only obtain a lower-order approximation solution. In order to 
avoid this complex one may make a transformation on the solution b+{P,t) 
to cancel the integral before integrating by parts the last two integrals. This 
transformation is given by 

bX{P,t)=bX{P,t)eM-i [ dhT+{P,h)]. (143) 

J to 

This transformation is the key point to the present equivalent transformation 
method to solve the basic differential equations (26). By this transformation 
and the initial condition (39) the transformed solution bX{P,t) may be written 

as ^ 

b].{P,t) = b-{P,t)t-^{P,t)exp[-i [ dt'2n{t')] 

J to 

+ibo{P,t)tl^{P,t)eiip[-iAP [ dt' Ko{t')]e^p[-i [ dt'Q{t')] 

Jto Jto 

-ibo{P,to)tl^{P,to)+i f dh{b-{P,ti)e-^{P,h)e^[-i r dt'2fl{t')]} 

J tn Jtn 
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+ f dti{bo{P,ti)e°_^{P,ti)ex.p[-iAP f' dt'KQ{t')]exp[-i f' dt'Q{t')]}. 

J to v to «' to 

(144) 

Now there is not any integral containing the solution b]_ {P, t) itself on the right- 
hand side of (144). The amplitudes of the solution b\{P, t) are related to those 
amplitudes of the solution h\{P,t) of (142) by the recursive relations: 

f;i(P,t) = r;i(P,t)exp[-i f dtiT+{PM)l (145a) 

J to 

f%(P,t) = r%(P,t)exp[-i f dtiT+{PM)l (145b) 

Jto 

Q+,{P,t) = [e^,{P,t) + T-,{P,t)r+{P,t)]eM-i I dtiT+{P,h)], (145c) 

Jto 

el,{P,t) = [el,{P,t)-Tl,{P,t)T+{P,t)]eM-i f dhT+{PM)]- (145d) 

Jto 

The transformation (143) and the integration by parts may be called the equiv- 
alent transformations as they does not generate any error term in these trans- 
formation processes. The transformation (143) does not improve essentially the 
first-order approximation solution, but it docs simplify greatly the calculation 
process to further obtain a higher-order approximation solution. Now by inte- 
grating by parts the last two integrals on the right-hand side of (144) and using 
the basic equations (26) and the initial condition (39) the solution h\{P,t) may 
be written as 



b\{P, t) = b\{P, t) = b-{P, i)r-2(P, t) exp[-i / dt'2n{t')] 

Jto 

-Fi6o(P,i)r°2(P,t)exp[-iAP / difKo{t')]exp[-i j dt'9.{t')\ 
-ibo{P,to)r%{P,to) + i f dh{bl{P,ti)T+{P,h)} 

Jto 

+i I dti{b-{P,ti)^+2{PM)ex.v[-i I dt'2n{t')\] 

Jto Jto 

dti{&o(i^,ii)6+2(-P,ii)exp[-iAP / dt'ifo(t')] exp[-i / dt'n{t')]} 

I to J to »' to 

(146) 

where the amplitudes satisfy the recursive relations: 

r;.(P,t) = r%(p,i) = f%(p,t)+^±i|M, (147a) 
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(147b) 

n-^p,^ .1 d , e+,{P,t) ^ 1 e(P,t)*9^i(P,^) 



«o ,p,,_ 1 9(fit)e;i(p,t) .a e%(^ ^oWe^i(P^ 

(147d) 

The unique difference between the equation (146) and the original equation 
(141b) for the variable b+{P,t) is that the initial condition (39) has been used 
in (146), while it is not used in (141b). Now the amplitudes r^(P, t), 6^2 *)> 
and 8*]_2(^ji) of the last three integrals on the right-hand side of (146) are 
inversely proportional to 0(4)^. Thus, these three integrals are secondary with 
respect to the first three terms on the right-hand side of (146). Moreover, the 
last two integrals contain the integrands with the largely oscillatory phase factor 
exp[— « J^^ dt'il(t')] or exp[— « J^^ dt'2^l{t')]. Then the integration by parts shows 
that these two integrals are secondary with respect to the other four terms on 
the right-hand side of (146). Thus, by the integration by parts the last two 
integrals becomes less important in the solution than before. 

The above equivalent transformations can be repeated many times that the 
solution b'^{P,t) {k = 1,2,...,) is transformed to the solution b^{P,t) by the 
equivalent transformation similar to (143) and then the solution b^{P,t) is 
changed to the solution b^^{P,t) by integrating by parts the last two inte- 
grals of the solution V^{P,t). The equivalent transformation from the solution 
b\{P,t) to the transformed solution b\{P,t) may be generally given by 

bX{P,t) = bX{P,t)eM-i I dtiT+{P,ti)], k= 1,2,.... (148) 

J to 

The solutions 6^(P, t) and &+(P, t) are called the /c— order exact solutions to the 
equation (141b) for the variable b+{P, t). From the solution b\{P, t) to the trans- 
formed solution b\{P,f) the amplitudes {T'l,^{P,f), Q+kiP^t)} {a = 0, -) of the 
solution b^iP, t) are transformed to the amplitudes {T'^}.{P, t), 6^^^(P, t)} of the 
solution bX{P,t) according to the recursive equations (145) if in the recursive 
relations (145) one makes the following replacements: r"i(P, t) T^i^{P,t), 

e%,{p,t) ^ e^fe(P,t), r«i(p,f) ^ r^jpt), e%,{P,t) ^ e^jpt), for 

a = 0, - , and r+ (P, t) ^ r+ (P, t) . On the other hand, from the solution b^ (P, t) 
to the solution b'^^{P,t) the recursive relations for their amplitudes are still 
given by (147) except the relation (147b), in which one needs to make the follow- 
ing replacements: f^i(P,t) ^ f^fe(P,t), e^i(P,t) ^ e^jP,^, r^2(^,0 ^ 
T%„+^{P,t), &%2{P,t) ^ &lk+iiP^t), for a = 0,-, and T+{P,t) ^ r+(P,0 
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and T^{P,t) ^ r'l_^_^{P,t). The relation (147b) is modified to the form 

r+ (PA A mn^+,ip,t) , 1 e{Rir^%(Rt) 



xexp[i/" dtiT+{P,h)]eMi f dhT+_^{P,h)]...exp[i f dtiT+{P,h)]. 

It can be found that after integrating by parts the last two integrals of the 
solution b^{P,t) many times, the two integrals become less and less important 
in the solution. Actually, it can turn out that the amplitudes ©^^.(P, t) and 

9'^j.(P, t) of the solution b'!^{P, t) is inversely proportional to the fc— th power of 

the Rabi frequency n{t), that is, e~,^{P,t) oc n{t)-'' and e%(P,i) a n{t)-''. 
Thus, by making only a few equivalent transformations of the integration by 
parts one c;an obtain a highly accurate solution to the original c;quation (141b) 
for the variable b+{P,t) even if the last two integrals are neglected. 

Now it is easy to obtain the first-order approximation solution to the orig- 
inal equation (141b) for the variable b+{P,t) from the equation (144) by ne- 
glecting the last two integrals on the right-hand side of (144). The first-order 
approximation solution is a linear algebra equation with the three variables 
{bQ{P,t) - 6o(P,to), b±{P,t)} and is given by 

b+iP,t) = b\{P,t) = b-iP,t)T-^{P,t)eiip[-i [ dt'2n{t')] 

Jto 

+i(5o(P,i)-6o(P,io))r°i(P,f)exp[-iAP / dt' Ko{t')]exp[-i f dt'^it')] 

J to Jto 

+iboiP,to)T°^i{P,t)e^[-iAP f dt' Ko{t')]exp[-i f dt'n{t')] 

J t(\ •/to 



ibo{P,to)Tl^{P,to)eMi [ dhT+{P,h)], (149) 

Jto 



'to 

while the truncation error is just the last two integrals: 



E+{P,t)=exp[tf dhr+{P,ti)] 

Jto 

f dti{b-{P,ti)&-i{P,ti)exp[-i r dt'2il{t')]} 

Jto Jto 



(iti{6o(P,ti)©+i(P,ii)exp[-iAP f' dt'Koit')]exp[-i f' dt'n{t')]}}. 

'to ^ to ^to 

By the integration by parts it can turn out that the error term E^-^{P,t) is 
bounded by 



i |e;i(P.^)l |e%(p,^)| |e%(p/o)i 
2 n{t) n{t) n{to) 
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+ f dti{\T+{P,h)\ + \Q^^{P,h)\ + \e%{P,h)\}. (150) 

Jto 

It is clear that this error upper bound is inversely proportional to f2(t)^. The 
second-order approximation solution is given by 

bl{P,t) =b-{P,t)T-2{P,t)exp[-i [ dt'2n{t')] 

Jto 

+i{bo{P,t)-bQ{P,to))T%{P,t)exp[-iAP f dt'Ko{t')]exp[-i f dt'n{t')] 

Jto Jto 

+i6o(^',io)r%(P,t)exp[-iAP / dt' KQ{t')]ex^[-i f dt/nit')] 

Jto Jto 

-ibo{P,to)Tl2iP,to)exp[i f dhT+{P,h)] (151) 

Jtn 



and the truncation error is bounded by 

ip+rp.^i< , \'dUP^t)\ , |e%(p,MI 



2 n{t) 0(i) 0(io) 

+ f dh{\r+{p,h)\ + \Q-^{p,t,)\ + \eUp,t,)\}. (152) 

Jto 

Obviously, this error upper bound is inversely proportional to fl{t)^ . The first- 
and second-order approximation solutions may be used to set up the adiabatic 
condition. 

By using the similar equivalent transformations mentioned above one may 
obtain the fc— order exact solution b^{P,t) {k = 1,2,...) from the equation 
(141b) for the variable b-{P,t), 

btiP, t) = b+{P, t)Ttk{P, t) exp[i / dt'2Q.{t')] 

Jto 

+ibo{P,t)T°_^{P,t)exp[-iAP f dt' Ko{t')]exp[i f dt'n{t')] 

Jto Jto 

-ibo{P,to)r\{p,to) + i I dh{b'i{P,h)r-{p,h)} 

Jto 

+t f dti{b+{P,ti)etk{P,ti)exp[i f dt'2il{t')]} 

Jto Jto 



+ I dti{bo{P,ti)e°_k{P,ti)exp[-iAP f' dt'Ko{t')]exp[i f' dt'n{t')]}. 

J to J to J to 

(153) 
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The equation (153) is equivalent to the original equation (141b) for the variable 
b-{P,t) if the initial condition (39) is taken into account. In particular, the 
first-order exact solution (P, t) = b- {P, t) and its five amplitudes are given 

The equivalent transformation from the solution 6^ (P, t) of (153) to the solution 
6^ (P, t) is given by 

lt{P,t) =bt{P,t)e-K^[-i [ dhT^{P,ti)]. (155) 

J to 

By the transformation (155) and the initial condition (39) the fourth term on 
the right-hand side of (153) is cancelled and the solution 6^(P, f) is changed to 
the transformed solution b'!L{P,t): 

rt 



VL (P, t) = b+{P, i)f +fe(P, t) exp[i / dt'2n{t')] 

J to 

+ibo{P,t)t°_k{P,t)exp[-iAP j rft'ii'o(t')] exp[i / dt'9.{t')\ 

J to J to 

-MP,to)t\{P,to)+i j dti{6+(P,ti)e+fe(P,ti)exp[i dt'2n{t')]} 
+ I dti{bo{P,ti)&°_k{P,ti)eyip[-iAP dt'Ko{t')]e^p[i f' dt'n{t')]}. 

J to "to "to 

(156) 

The recursive relations for the amplitudes of both the solutions b'^ (P, t) and 
6^ (P, t) are given by 

fifc(P,i) = Tt,{P,t)eM-i fdhT^{P,h)], (157a) 

Jto 

f°fe(P,t) = r°fe(P,t)expH / dhT-{P,h)], (157b) 

Jto 

eU{P,t) = [e+,(p,t) + r+,(p,t)r-(p,t)]expH f dhT-{p,h)], (i57c) 

Jto 
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e\{P,t) = [e\iP,t) - T°_,{P,t)T^iP,t)]eM-i fdUTliPM)]. (157d) 

After making the integration by parts on the last two integrals on the right- 
hand side of (156) the solution b'^{P,t) is changed to the {k + 1)— order exact 
solution b'^^{P, t) which is also given by (153). The recursive relations between 
the amplitudes of both the solutions b'^{P, t) and h^^{P, t) are generally given 

by 

Tt^,^,){P,t)=tt,{P,t) + \^^^, (158a) 
r°(,+i)(P,t) = t\{P,t) - (158b) 

,„,,_,i r(P,t)9^,(p,t) ^ 1 e(p,t)*e°,(p,t) ^ 

X I dhT+{PM)]eMi I dtiT+_^iP,h)]...eMi [ dhT+{P,h)], 

J to -J Iq J tQ 

(158c) 

n+ (P,^ -^d fiU{P,t) 1 Q{P,tYQ\{P,t) 
(p,. 1 9(p,^)el,(p,^) d Q%{p,t) Mt)e'_,{p,t) 

0_(fe+i)(P,i) - ^ + z-{—^^) + — AP 

(158c) 

The recursive relations (157c), (157d), (158d), and (158c) show that in the last 
two integrals of the fc— order exact solution b^{P,t) of (156) the amplitudes 
0^j.(P, f) and 9°j.(P, t) are inversely proportional to the A;— th power of the 
Rabi frequency n{t), that is, etfe(P,f) a n{t)~'' and Q\{P,i) a Vl{t)-^ . 
Thus, by making a few equivalent transformations one may obtain a high-order 
approximation solution from the exact solution 6^(P, t) or b'^{P,t). The first- 
order approximation solution may be obtained from the exact solution b^ (P, t) 
by neglecting the last two integrals on the right-hand side of (156) {k = 1), 

b-{P,t) = bl{P,t) = 6+(P,t)r+i(P,t)exp[i / dt'2n{t')] 

J to 

-hi(6o(P,i) - 6o(P,io))r°i(P,t)exp[-iAP / dt' Ko{t')]exp[i [ dt'n{t')] 

Jto Jto 

-Fi6o(P,io)r°i(P,t)exp[-iAP / dt' Ko{if)]exp[i [ dt'n{if)] 

Jto Jto 

-ibo{P,to)T°_^{P,to)exp[i [ dhT^{P,ti)], (159) 

Jto 
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and the truncation error is just given by these last two integrals neglected, and 
by the integration by parts it can turn out that the truncation error is bounded 

by 

l^.il^,fJ|S2 + j^(^) + ^^^^^ 

+ / dti{|r2-(P,ti)| + |et2(P,ti)l + |e°2(P,ii)|}. (leo) 

Jta 

The error upper bound is proportional to Q.{t)~^ . The second-order approxima- 
tion solution is obtained from the exact solution (P, t) by neglecting the last 
two integrals on the right-hand side of (156) (A; = 2), 



6?. (P, t) = b+{P, t)Tt2iP, t) exp[i j dt'2n{t')] 

J to 

+ibo{P,t)r%{P,t)exp[-iAP f dt'Ko{t')]exp[i / dt'n{t')] 

Jto Jto 

-ibo{P,t„)r'l^{P,to)cxp[i f dhT^iP,h)] (161) 

Jto 

and the truncation error is bounded by 

i |e±2(P,^)l , \e'-2{P,t)\ , |e°,(p,to)| 
\^r2{P^*>\ ^ 2 n{t) + Q(i) + n{to) 

+ f dh{\T^{p,h)\ + \ets{P,h)\ + \e%{p,tr)\}. (162) 

Jto 

The error upper bound is proportional to ^{t)~^. The first- and second-order 
approximation solutions may be used to set up the adiabatic condition below. 

With the help of the equivalent transformations similar to those used above 
one may obtain the fc— order exact solution for the variable 6o(P-i) from the 
equation (141a). At first the first-order exact solution may be given by 

bo{P, t) - bo{P, to) = ib+{P, i)rji (P, t) exp[iAP / difKo{t')] exp[i f diffl{t')] 

Jto Jto 

+ib-{P,t)To^{P,t)exp[iAP f dt' Ko{t')]exp[-i ( dt'9.{t')\ 

J tn J to 



h / dti{b+{P,ti)@^^{P,ti)exp[iAP r dt'Ko{t')]exp[i dt'n{t')]} 

J to J to J to 

/ dti{&-(P,ti)eoi(P,ti)exp[iAP / dt'Ko{t')]exp[-i / dt'n{t')]} 

J to J to J to 

(163) 
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where the four ampUtudes are given by 

1 Q{P,t)* 



T+,{p,t) = -r^,{p,t) 



V2 n{t) 



Q+fPi)- Q-rpi^-.9 MP,tr i T{p,t)e{p,tr Ko{mp,tr 

Uoi(P,f)- i)o,[P,t)-tg^( ^^^^ ^^^^ ^^^^ AP. 

A special point in the first-order exact solution (163) is that there is not any 
term containing the solution 6o(^) t) itself on the right-hand side of (163). Thus, 
one may make directly the integration by parts on the last two integrals on the 
right-hand side of (163). Here still denote that Sq{P, t) = 6o(P, t) — 6o(-P, to) and 
5^(P, t) = b±{P,t) — b±{P,tQ). Then the second-order exact solution d^{P,t) 
may be expressed as 



6f{P, t) = 6oiP, t) = ib+{P, i)r+ (P, t) exp[iAP / dt'Ko{t')] exp[i / dt'n{t' 

Jto Jto 

-|-i6_(P,i)r- (P,f)exp[iAP f dt' Ko(t')]exp[-i [ dt'n(t')] 

J la Jto 
rt 



+ibo{P,to)T9,°{P,t) + i [ dt^{dfiP,ti)T''^{P,ti)} 

Jto 

+ I dii{6+(P,ti)ej2(P,ii)exp[iAP f' dt'Koit')]exp[i f' dt'n{t')]} 

Jto J to J to 

+ j dti{b-{P,ti)eo2{P,ti)exp[iAP f' dt'Ko{t')]exp[-i f' dt'n{t')]} 

J to J to J to 

(164) 

where the amplitudes are given by 

r+ (p,t) = r+ (p,t) - T^,{p,t) = ToKP.i) + 

Tf{P,t) = j'dt.TliPM), T%{P,t) = -V2^i^L^^iM, 



dt^ n{t) 



9+, (P,i) = -eo-,(P,t) = ^|:(%S^) 



j^o(t)e+ (p, t) 1 r(p, t)Q^^ (p, t) 
n{t) 2 n{t) 

Now the first-order approximation solution may be obtained from the first-order 
exact solution (163) by neglecting the last two integrals on the right-hand side 
of (163). It is given by 

5o{P,t) =ib+{P,t)r+^{P,t)exp[iAP f dt'Ko{t')]exp[i f dt'n{t')] 

Jto Jto 
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+ ib-{P,t)Toi{P,t)exp[iAP [ dt'Ko{t')]eyip[-i [ dt'n{t')] (165) 

and the truncation error is just the last two integrals on the right-hand side of 
(163), and it can turn out by the integration by parts that the truncation error 
is bounded by 

K°i(p,i)| < ^^|^ + j^''(iti{|r^(P,ti)| + 2|e+(P,ti)|}. (166) 

This error upper bound is clearly proportional to Q.{t)~'^ . A higher-order ex- 
act solution than the second-order one (164) may be obtained by integrating 
by parts the last two integrals in (164), but it needs first to eliminate the 
fourth term that contains the solution 5^{P,t) itself on the right-hand side 
of (164). This can be done by making an equivalent transformation on the so- 
lution 5^q{P, t). In general, for the fc— order exact solution 6q''{P, t) {k = 2, 3, ...) 
the equivalent transformation is written as 

'C{P,t) = 5^'=(P,t)exp[-z f dHTl{P,H)]. (167) 



Then the fc— order transformed solution 5q (P, t) may be expressed as 

-6fe ' /■* /■* 

6o {P,t)=ib+{P,t)T+^{P,t)exp[iAP dt' Ko{t')]e^p[i dt'n{t')] 

Jto Jto 

+i6_(P,i)fo^(P,i)exp[iAP / dt'KQ{t')]exp[-i [ dt'n{t')]+ibQ{P,tQ)tl°{P,t) 

Jto Jto 

+ f dti{6+(P,ti)e+ (P,ii)exp[iAP r dt'Ko{t')]exp[i dt'n{t')]} 

J to J to J to 

+ f dh{b-{P,ti)%^{P,h)exp[iAP r dt'Ko{t')]exp[-i dt'n{t')]} 

J to to to 

(168) 

where the amplitudes satisfy the recursive relations: 

^okiP, t) = ^okiP^ t) exp[-i f dhTliP, ii)], (169a) 

Jto 

f ofc(^. *) = ^okiP^ t) exp[-i /* dhVliP, h)], (169b) 

Jto 

tl\P,t) = f dh{[-^rl°{P,ti)]exp[-i r dt'rl{P,t')]}, (169c) 

Jto Oti 

t) = [Q^kiP, t) - r+ (P, t)rl{P, t)] exp[-» / dhFliP, ii)], (169d) 

Jto 
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0ofe(^. t) = [QokiP^ t) - T^^iP, t)TUP, t)] exp[-i r dt^TliP, h)]. (169e) 

J to 

From the fc— order transformed solution 6q {P, t) to the (fc + 1)— order sohitfon 
5f+\P,t) there is an equivalent transformation of the integration by parts. The 
{k + 1)— order exact solution 5^f^^^{P,t) is still given by (164) as long as one 
makes the replacement: 5^Q{P,t) 5^Q^^{P,t) and the following replacements 
for the amplitudes: 

n%p,t) - n%,{p,t), rO(P,t) - Tl^,{P,t). 



The recursive relations for the amplitudes between the two solutions Sq (P, t) 
and S^o''+\P,t) {k > 2) are given by 

Y+ (pt)-r+(Pt) ^okjP'^) p- (p ,\ _ f- (p .\ I Qofc(-P^^) 

(170a) 
(170b) 

po .p,^_r 1 e(p,^)9+(p,t) 1 e(p,t)e„,(p,^) 

xe^i f dhTl{P,h)]eM^ j dhTl_^iP,h)]...e^p[t f dhT9.{P,h)], 

J to J tQ v to 

(170c) 

«+ f p ,^ - . ^ ^ Qo/c(^' ^) ^ -^oft)^^il^(^- A p 1 r(p /)e,-,^(P, t) 

•^ofc+il^^iJ ^t) > n{t) 2 0(i) 

(170d) 
(170e) 

It can turn out that the amplitudes 6qj,(P, t) and Ogj,(P, t) of the fc— order exact 
solution (5q'^(P, t) are inversely proportional to Q{t)~''. Thus, the last two inte- 
grals of the fc— order exact solution S^Q^P^t) have a negligible contribution to 
the solution if the Rabi frequency ^l{t) is large. Now the second-order approxi- 
mation solution may be obtained from the second-order exact solution 5q (P, t) 
of (168), 

6'^^{P,t) =ib+{P,t)T+2{P,t)exp[iAP [ dt'Ko{t')]exp[i [ dt'n{t')] 

Jto Jto 

+ib-{P,t)To2{P,t)exp[iAP f dt' Ko{t')]eyip[-i f dt'n{t')] 

Jto Jtn 
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-6o(P,io){l-exp[i / rftirO(P,ti)]}, (171) 

Jta 

while the truncation error is given by the last two integrals on the right-hand 
side of (168), and the truncation error upper bound is determined from 

^^r2{P,t)\< ^^^^ + ^^^^ 

+ f dh{\Tl{p,t)\ + \e+,{p,t)\ + \e^,{P,t)\}. (172) 

J to 

This error upper bound is clearly proportional to fl{t)~^. 

Now by solving the three first-order approximation solutions (149), (159), 
and (165), which are all linear algebra equations, one may obtain the first- 
order approximation solution to the basic equations (141). At first according 
to the first-order approximation solutions (149), (159), and (165) and their 
truncation errors the first-order exact solution to the equations (141) may be 
formally written as, (this is really the first-order approximation solutions plus 
their truncation errors), 

b+{P,t) = atb-{P,t) + a+dl{P,t)+f3++E+{P,t), (173a) 

6_(P,t) = a+b+iP,t) + aoSliP,t) + f3o + E;,iP,t), (173b) 

S''„{P,t) = alb+{P,t) + a'Lb-{P,t) + E°i{P,t), (173c) 

where the truncation errors E^-^ {P, t) and E^-^ [P, t) may be considered as small 
parameters and any other parameters such as ali, a^, a5}_, etc., can be obtained 
directly from the first-order approximation solutions (149), (159), and (165). 
These three equations are linear algebra equations as E^-^ {P, t) and E^-^ (P, t) 
are considered as the small parameters. The exact solution to the three linear 
algebra equations (173) may be written as 

5l{P, t) = S±m{P, t) + Ef^{P, t), S'oiP, t) = Som{P, t) + E^,{P, t), (174) 

where the main terms are written as 

S±m{P,t) = Ti^bo{P,to)F{P,t) 

x{^^^exp[-iAP f dt'Ko{t')]eyip[Ti f dt'n{t')] 

^H*) J to J to 

Q{P,to) .1 /■* . r(p,ti)^ + 4|9(P,ti)p ^^ 

domiP,t)=0. (175b) 
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Here the factors F{P,t) and G±{P,t) are defined by 



F{P,t) = 



2 



, G±{P,t) = 



1 r(p,t) , 1 |e(p,t)|^ 
4 a(t) 2 ~n(tp~ 



1 + 



1 r(p,ty 
16 n(t)^ 



+ 



1 _L 1 r(P,t)^ , |e(p,t) 




The factors G±{P,t) appear in the error terms E^^{P,t) and E^^[P,t), as can 
be seen below. Then the upper bound of the error term that is generated by 
the main terms 5±rn{P-,i) and 5om(-P, i) is determined from 



It is clear that this upper bound is proportional to Q.{t) ^ . Notice that the factor 
F{P,t) is unity approximately if |e(_P,f) I Vi^(t)^ << I and T{P,tf /D.{tf « 1. 
By comparing the inequality (176) with the inequality (133) one can find that 
if the factor F{P, t) is unity, then the inequality (176) is really the inequality 
(133) that leads to the first-order approximation adiabatic condition (134). 

Now investigate the secondary error terms E^^{P,t) and E^^{P,t) in the 
exact solution (174). These error terms contain the first-order truncation errors 
£^i(P,i) and £'°i(P,t) (See: (150), (160), and (166)). They may be written as 



£0,(P, t) = E°UP, t) + E°UP, t) + E°t{P, t) + E°UP, t), (177b) 

where (P,t), E^UP,*) cx ^(t)-^; (P,i), E^fiP,*) cx n{t)-^; Ef^{P,t), 

El!f{P,t) a n{ty^; E';*.f{P,t) oc n{t)-^ . It is not difficuh to obtain the upper 



bounds for all these error terms, but one needs only to consider the dominating 
terms Ef'^{P,t) and E^f{P,t) that are proportional to Q{t)^^, since the trun- 
cation errors Ef^{P,t) and E^i{P,t) are proportional to n(t)~^ and the other 
error terms E^^''{P,t) and E°l{P,t) (fc = 3,4; / = 3,4,5) are higher-order and 
can be neglected with respect to the dominating terms E^^{P, t) and E^f{P, i). 
The dominating error terms are given by 



\\El^\PM = ^/\S+m{P,tr + \S-miP,tr + \S0m{P,tr 




(176) 



£±(P, t) = EtfiP, t) + EtfiP, t) + Et,\P, t), 



(177a) 






(178a) 



E^-^iP, t) = E°,{P, t) - bo{P, to)F{P, t) 



|e(p,t)p 
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xcos(/;...[ofa)+ ^'^-'-';;;|^'^-'-"' ]}. 

The last term on the right-hand sides of (178) for each one of the error terms 
Ef^^{P, t) and E°'^{P, t) is proportional to the value e(P, tQ)/Q.{to) at the initial 
time <o of the basic STIRAP decelerating or accelerating process. As discussed 
before (See: (133) and (134)), the initial value |0(P, io)|/^^(^o) may be controlled 
to be so small that it can be neglected. Then the norm for the error vector 
{E^j{P,t),E+\P,t),E-^\P,t)Y is bounded by 



\\Ei^\PM = \\Ett\P,t)\^ + \E;,\P,t)\^ + |S02(P,f)|2 



< 



^|E0,(P,t)|2 + \F{P,t)n\E+,{P,t)\'^ + |£-i(P,i)|2) 



This error upper bound is proportional to Vl{t)~'^ . Now it follows from (174) 
and (177) and then the inequalities (176) and (179) that the total deviation of 
a real STIRAP adiabatic process from the ideal one at any instant of time is 
bounded by 



\\Er{PM = ^J\5liP,t)\^ + \6''_{P,t)\^ + \6UP,tW 
<\\E('\P,t)\\ + \\E?\P,m 
- + \/\E?iiP,tW + \F{P,m\E+{P,t)\^ + \E;,{P,t)\^) 



^mp,t)\ 



\^ yP'*)\ —PTTI^ + 1 r(Pf.\2 |S77?7TT2— , (li^U) 



where the upper boimds for the truncation errors \E^j^{P,t)\, \E^j^{P,t)\, and 
\E^-^^{P,t)\ are obtained from (150), (160), and (166), respectively. In order to 
obtain the global adiabatic condition one needs to limit the maximum value on 
the rightcst side of (180) not to be more than a desired small value over the 
effective momentum region [P] and in the whole time period to < t < to + T 
of the STIRAP decelerating or accelerating process. Denote Ad{P,t) as the 
function on the rightest side of (180). Then the global adiabatic condition may 
be expressed as 

\\Er{P,t)\\ < MP,t) < max {Ad{P,t)} < Er, (181) 
-Pe[P], to<t<to+T 
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where Sr is the desired vahie and Sr << 1- Unhke the adiabatic condition 
(128) there is not an exponential correction factor in the adiabatic condition 
(180) and (181). The adiabatic condition (180) and (181) consists of the first- 
order term that is proportional to f7(f)~^ and the second-order terms that are 
proportional to f2(t)~^, while the adiabatic condition (128) is the first-order 
term with the factor F{P, t) = 1 times the exponential correction factor. Thus, 
the adiabatic condition (180) and (181) is much less severe than the adiabatic 
condition (128), the latter is most severe for a quantum ensemble with a broad 
momentum distribution. However, just like the adiabatic condition (128) the 
adiabatic condition (180) and (181) is also strict and accurate (It is not difficult 
to include the omitted secondary error terms of (177) in the adiabatic condition 
(180) and (181)). Therefore, if any one of the two general adiabatic conditions 
is met, then a better STIRAP pulse sequence is obtained and the perfect state 
(or population) transfer may be realized by the STIRAP pulse sequence. The 
adiabatic condition (180) and (181) may be more useful in practice to realize 
the perfect STIRAP state (or population) transfer in an atomic or molecular 
ensemble with a broad momentum distribution. It may also be used to realize 
the STIRAP decelerating and accelerating processes in the laser cooling and 
the quantum coherence interference experiments of a cold atomic or molecular 
ensemble. 

The adiabatic condition (180) and (181) is still slightly severe. This is be- 
cause it is based on the first-order approximation solution to the basic equations 
(141), leading to that the upper bounds for the truncation errors \E^-^{P,t)\, 
\E~^{P,t)\, and \E^i{P,t)\ are not the lowest ones. As known before, these trun- 
cation errors are proportional to Q.{t)^'^ . A better adiabatic condition may be 
set up on the basis of the second-order approximation solution to the basic equa- 
tions (141), which are given by (151), (161), and (171), and the corresponding 
truncation error upper bounds |_E^(P, i)|, [^^^"^(P, t)|, and |i?J?2(P, t)| arc given 
by (152), (162), and (172), respectively. These upper bounds are proportional 
to Q.{t)~^. Thus, they are the lower ones with respect to those truncation error 
upper bounds of the first-order solution, leading to that the adiabatic condition 
based on the second-order approximation solution is better one. 

It is known that in the adiabatic condition (180) and (181) the truncation 
errors \E\(P,t)\, \E~^{P,f)\, and \E^.^{PJ)\ are proportional to n{t)-'^ and 
the factor F{P,t) is a function of the ratio T{P,t)^ /U,(t)^ . Then the parame- 
ter T{P,t) appears only in those terms that are proportional to Q.{t)~'^ on the 
rightest side of (180). This is different from the case that the parameter 6(P, t) 
may appear in the first term on the rightest side of (180) that is proportional 
to Q.{t)~^. Thus, the parameter T{P,t) has a smaller contribution to the adi- 
abatic condition (180) and (181) than the parameter ld{P,t). As pointed out 
before, the effect of the momentum distribution is dependent upon whether the 
Raman laser light beams of a STIRAP experiment are copropagating or coun- 
terpropagating. Consider that the Raman laser light beams are copropagating. 
Then the momentum distribution could have a relatively small effect, because 
in this case the parameter Q{P,t) has a smaller value. This is correct for the 
first-order approximation adiabatic condition (134). However, it could not be 
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so simple for the inequality (180) in a quantum ensemble with a broad momen- 
tum distribution. When the first term on the rightest side of (180) is smaller 
due to a small parameter Q{P,t), the second and third terms could become 
more important and have a dominating contribution to the adiabatic condition 
(180) and (181). Then in this case the parameters T{P,t) and Ko{t) become 
more important in the adiabatic condition (180) and (181), resulting in that the 
momentum distribution has a large effect on the adiabatic condition. There- 
fore, the momentum distribution needs to be considered explicitly even for a 
conventional STIRAP state (or population) transfer experiment in a quantum 
ensemble with a broad momentum distribution, which uses the copropagating 
Raman laser light beams. Obviously, for the STIRAP decelerating and acceler- 
ating processes in a quantum ensemble with a broad momentum distribution, 
which use the coimterpropagating Raman laser light beams, one needs to con- 
sider generally the effect of the momentum distribution on the STIRAP state 
(or population) transfer. 

The present adiabatic theoretical methods including the equivalent trans- 
formation method to solve the basic differential equations (26) not only can be 
used to set up a general adiabatic condition for the basic STIRAP decelerating 
and accelerating processes of a single atom or molecule or a quantum ensemble 
of the atoms or molecules, but also they will have an extensive application in 
other research fields such as the NMR spectroscopy (See: for example, Ref. [41]) 
and the magnetic resonance image (MRI). 

8. Discussion 

In the paper the standard three-state STIRAP population transfer theory 
in the laser spectroscopy [4, 15, 16, 17, 18] has been developed to describe theo- 
retically the STIRAP-based unitary decelerating and accelerating processes of a 
single freely moving atom by combining the superposition principle in quantum 
mechanics [25] and the energy, momentum, and angular momentum conserva- 
tion laws for the atomic photon absorption and emission processes [5, 19]. There 
are similar theoretical works or developments to describe the atomic laser cool- 
ing process [5, 19, 20, 21] in a neutral atom ensemble and the atomic quantum 
interference experiments [10, 12] in a cold atomic ensemble. There are also a 
number of works to investigate the atomic decelerating and accelerating pro- 
cesses by the laser light techniques [21, 23, 24]. However, the present work 
is focused on the analytical and quantitative investigation how the momentum 
distribution of a superposition of the momentum states of a pure-state quantum 
system such as a single freely moving atom affects the state-transfer efficiency 
in these STIRAP unitary decelerating and accelerating processes. It empha- 
sizes the complete STIRAP state transfer and the unitarity of these processes. 
This means that in the present work any decoherence efi'ect of the atomic sys- 
tem under study is not considered. A main purpose to investigate the effect of 
the momentum distribution on the STIRAP state transfer is to build up bet- 
ter STIRAP unitary decelerating and accelerating sequences, so that the time- 
and space-compressing processes of the quantum control process to simulate 
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the reversible and unitary state-insensitive halting protocol [22] ean be realized 
through these decelerating and accelerating processes. Thus, this is involved 
in setting up a general adiabatic condition for the STIRAP unitary decelerat- 
ing and accelerating processes. In the paper two general adiabatic conditions 
have been obtained analytically, one based on the Dyson series solution to the 
basic differential equations to govern the STIRAP processes, another based on 
the equivalent transformation method to solve the basic differential equations. 
Both the general adiabatic conditions may be used to set up a conventional 
STIRAP experiment and also the STIRAP-based decelerating and accelerat- 
ing processes. A complete STIRAP state transfer could be achieved only in 
the ideal adiabatic condition. Generally, it is hard to achieve a complete state 
transfer for the STIRAP processes when the atomic momentum superposition 
state has a broad momentum distribution. However, in the ideal or nearly ideal 
adiabatic condition an almost complete STIRAP state transfer may be realized 
if the superposition of the momentum states has a small effective wave-packet 
spreading or a narrow momentum distribution. When the initial motional state 
of a freely moving atom is a Gaussian wave-packet state, the final motional state 
of the atom is still a perfect or almost perfect Gaussian wave-packet state after 
the atom undergoes the STIRAP unitary decelerating (or accelerating) process 
in the ideal or nearly ideal adiabatic condition. Therefore, in the paper it is 
shown that the time- and space-compressing processes of the quantum control 
process [22] can be realized almost perfectly through the STIRAP decelerating 
and accelerating processes in the ideal or nearly ideal adiabatic condition. This 
is one of the important results in the paper. 

The standard STIRAP population transfer theory is generally based on the 
semiclassical theory of electromagnetic radiation. In the semiclassical theory the 
externally applied electromagnetic fields such as the Raman laser light beams 
are considered as the classical electromagnetic fields, while the atomic system 
itself and the interaction between the atomic system and the external electro- 
magnetic fields are treated quantum mechanically. It has been shown that the 
semiclassical theory can describe almost perfectly the three-state STIRAP pop- 
ulation transfer experiments of the atomic and molecular beams in the laser 
spectroscopy [4, 15, 16, 17, 18]. On the other hand, the semiclassical theory 
is also successful to describe the STIRAP-based laser cooling processes in a 
neutral atomic ensemble [20, 21] and especially the atomic quantum interfer- 
ence experiments in a cold neutral atom ensemble [10, 12, 13, 14]. In the paper 
the semiclassical theory also is directly employed to describe the STIRAP-based 
unitary decelerating and accelerating processes of a single atom. The semiclassi- 
cal theory of electromagnetic radiation generally can not explain reasonably the 
atomic spontaneous emission [1, 25. 26]. However, it is suited to describe these 
STIRAP unitary decelerating and accelerating processes due to that these STI- 
RAP processes can avoid the atomic spontaneous emission by setting suitably 
the experimental parameters. 

As far as the Hamiltonian of Eq. (4) to describe the three-state STIRAP 
experiment of an atom system is concerned, there are three requirements: (i) 
the three-state subspace for the atomic internal states is closed under the Hamil- 



87 



tonian; (ii) the clcctric-dipolc approximation is satisfied; and (Hi) the rotating 
wave approximation (RWA) is reasonable. The first requirement can be sat- 
isfied if one chooses suitably the atom and its three internal states and the 
experimental parameters of the Raman laser light beams. Since the size of an 
atom is generally much less than the wave lengths of the Raman laser light 
beams at the optical frequencies (~ uiqi and W02), the second requirement may 
be met generally. If the Raman laser light beams are strong, it can turn out 
that in the first approximation the strong laser light field may generate a Bloch- 
Siegert shift to the transition frequency of the atomic internal energy levels [1]. 
When the Rabi frequencies {flp{t) and ils{t)) of the Raman laser light beams 
are much less than the resonance frequencies (wqi and 0^02) of the atomic inter- 
nal energy levels and the detunings for the Raman laser light beams are small, 
the magnitude of the Bloch-Siegert shifts generated by the Raman laser light 
beams at the optical frequencies is very small and may be negligible and hence 
the rotating wave approximation is reasonable [1]. However, it is also conve- 
nient to correct the Bloeh-Siegert shifts in the STIRAP experiment because one 
needs only to add the Bloch-Siegert shifts to the resonance frequencies (wqi 
and 0)02 )• A simple evaluation for the Bloch-Siegert shifts can be seen in Ref. 
[1] and a general treatment may use the average Hamiltonian theory [34]. On 
the other hand, it is also possible to apply an extra laser light field for each 
one of the Raman laser light beams in the STIRAP experiment to compensate 
the rotating-wavc approximation. In fact, if each one of the two Raman laser 
light beams in the STIRAP experiment is replaced with a pair of the laser light 
beams with the orthogonal electric field vectors and the suitable phases [38], 
one can eliminate the rotating-wave approximation. Similarly, one may also use 
the circularly polarized lights to prepare the dipole interaction Hamiltonian (10) 
without the rotating wave approximation [19, 38]. 

There is also another condition to be met that the electromagnetic field of 
any Raman laser light beam is considered as an infinite plane- wave electromag- 
netic field when calculating the time evolution process of a freely moving atom 
under the STIRAP unitary decelerating and accelerating processes. Here the 
infinite plane-wave electromagnetic field has spatially uniform amplitude and 
phase. This condition can be met only when the electromagnetic field can en- 
compass sufficiently the whole wave-packet motional state of the moving atom. 
Note that the electromagnetic field of the Raman laser light beam propagates 
along a direction parallel to the atomic moving direction in one-dimensional 
space. Since the atom is in a Gaussian wave-packet motional state which has a 
finite wave-packet spreading, then one can set suitably the experimental param- 
eters for the Raman laser light beam such that the electromagnetic field in space 
is much wider than the effective wave-packet spreading of the atomic motional 
state during the whole decelerating or accelerating process. Then in this case 
the electromagnetic field in space can be reasonably considered as an infinite 
and uniform plane-wave electromagnetic field for the atomic wave-packet mo- 
tional state. The condition may be satisfied more easily for a heavy atom as 
the wave-packet motional state for such atom has a more narrow wave-packet 
spreading. If the electromagnetic field in space has a finite bandwidth less than 
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or comparable to the wave-packet spreading of an atomic motional state, then 
it can not be considered as an infinite and uniform plane- wave electromagnetic 
field and the electric dipole interaction of Eq. (10) with space-independent Rabi 
frequencies is not suited to describe the STIRAP process, since in this case the 
Rabi frequencies are dependent upon the spatial coordinate [35]. 

As an important result in the paper, it is shown that if a free atom is in 
a Gaussian wave-packet motional state at the initial time, then it is still in a 
Gaussian wave-packet motional state after it is decelerated (or accelerated) by 
the STIRAP-based unitary decelerating (or accelerating) sequence in the ideal 
or nearly ideal adiabatic condition. As far as a Gaussian wave-packet state is 
concerned, there are two types of time- and space-dependent unitary evolution 
processes that do not change the Gaussian wave-packet shape of the atomic 
motional state. The first type is that the unitary evolution processes may ma- 
nipulate and control the center-of-mass position and/or momentum of a Gaus- 
sian wave-packet state but can not manipulate at will the complex linewidth 
of a Gaussian wave-packet state. The second type is that the unitary evolu- 
tion processes may manipulate and control the complex linewidth of a Gaussian 
wave-packet state. The STIRAP-based unitary decelerating and accelerating 
processes belong to the first type. This type of the unitary evolution processes 
tend to have the property that in the unitary evolution process the imaginary 
part of the complex linewidth of a Gaussian wave-packet state increases linearly 
with the time period of the unitary evolution process, while the real part usually 
keeps unchanged, i.e., the wave-packet spreading of the Gaussian wave-packet 
state becomes larger and larger early or late as the time period increases. These 
unitary evolution processes which have the property also include the free-particle 
motion and the atomic bouncing process off a hard potential wall in the special 
case. Thus, the free-particle motion and the atomic bouncing process may be 
assigned to the first type. The Hamiltonians of the quantum systems to create 
the first type of the time- and space-dependent unitary propagators usually can 
not be singly used to generate their inverse unitary propagators without any 
help of the interactions from outside the quantum systems. Therefore, this type 
of unitary evolution processes do not have their own inverse unitary propaga- 
tors in these quantum systems separated from the outside or their environment. 
Obviously, these separate quantum systems also include the isolated quantum 
systems in the quantum statistical physics and they are described completely 
by the Hamiltonians of the quantum systems. These unitary evolution pro- 
cesses could be considered to be self-irreversible in the sense that there do not 
exist their own inverse unitary propagators in the same separate quantum sys- 
tems, although these processes obey the unitary quantum dynamics and their 
inverse unitary propagators could be generated with the help of the specific 
interactions from outside the quantum systems. Take a free-particle motion as 
a typical example. A free-particle motion can be described completely by the 
unitary propagator U{t) = ex-p[—ip'^t/{2mH)]. Of course, one may also choose 
the unitary propagator U{t)^ = exp[ip^t/(2mfi)] to describe the free-particle 
motion. However, once one chooses one of the two unitary propagators to de- 
scribe the free-particle motion, another is the inverse propagator of the chosen 
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unitary propagator. The free particle may move toward the left or the right in 
the coordinate axis and both the motions can be described by the same unitary 
propagator U{t). Now the inverse free-particle motion is described by the uni- 
tary propagator U(t)^. The Hamiltonian H ~ p'^/2m of the free particle can 
really generate only the unitary propagator U{t) (or U{t)~^) but can not really 
generate both the unitary propagator U{t) and its inverse propagator U{t)~^ 
simultaneously by the Schrodinger equation ihdU{t)/dt = H{t)U{t). One may 
take a hermite conjugate on the Schrodinger equation and then could obtain 
f7(t)+, but taking a Hermite conjugate is not a real physical process and hence 
there does not exist the inverse propagator U{t)~^ in the free-particle quantum 
system. Then the inverse free-particle motion will never really take place if the 
free particle is not acted on by any external interaction. In other words, one 
may argue that the time reversal process (i.e., the inverse frec-particlc motion) 
may be a real physical process [25, 30], but this process can not really take 
place for the free particle without any help of the specific external interactions. 
Then in this sense the free-particle motion is really self-irreversible, although 
it is governed by the unitary propagator U{t). Of course, it is possible to con- 
struct the inverse unitary propagator of the free particle if the specific external 
interactions are applied to the free particle. How to construct these inverse 
propagators of the type of time- and space-dependent unitary propagators men- 
tioned above will be reported in next paper. As pointed out in the previous 
papers [22, 36], such a situation that a quantum system that obeys the uni- 
tary quantum dynamics can not really have both the unitary propagator and 
its inverse propagator is often met in the quantum systems which are used to 
implement a quantum computation. The spontaneously irreversible processes 
of isolated quantum systems in the quantum statistical physics [40] could be 
related to the situation. It has been stressed in the previous papers [22, 36, 37] 
that these irreversible processes should be understood through the unitary quan- 
tum dynamics instead of the stochastic process and probability statistics [40], 
and they could be related to the difference between the unitary evolution pro- 
cess and its inverse process. On the other hand, there also exist other quantum 
systems that have both the unitary propagators and their inverse propagators, 
where both the types of the unitary propagators can be generated by the same 
Hamiltonians of the quantum systems. A conventional harmonic oscillator is 
typically one of such quantum systems. In the quantum statistical physics such 
quantum systems obey completely the Poincare's recurrence theorem. 

There arc the second type of time- and space-dependent unitary propagators 
that can manipulate and control the complex linewidth of a Gaussian wave- 
packet motional state. A general quadratic Hamiltonian can generate a time- 
and space-dependent unitary propagator that can keep the Gaussian shape un- 
changed for a Gaussian wave-packet motional state under the action of the 
unitary propagator. Obviously, the Hamiltonian is different from those of the 
STIRAP unitary decelerating and accelerating processes. One needs this type 
of unitary propagators to manipulate and control the complex linewidth of a 
Gaussian wave-packet motional state to build up the quantum circuit for the 
reversible and unitary state-insensitive halting protocol and also needs their in- 
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verse unitary propagators to realize the efficient quantum searcli process [22, 
36] . How to construct these unitary propagators will be reported in next paper 
(Arxiv: quant-ph/0708.2129). 
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